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ABSTRACT 


This  thesis  considers  the  problem  of  synthesizing  cyclic  (or 
periodic)  sequences,  with  binary  and  nonbinary  (N-ary)  symbols,  for 
special  cammunication  applications.  For  the  applications  considered, 
the  sequences  are  transmitted  over  channels  with  additive  white  Gaussian 
noise.  For  correlation  detection  the  autocorrelation  function  of  a  se¬ 
quence  or  the  cross -correlation  function  between  sequences  are  required 
for  the  evaluation  of  performance . 

The  first  part  of  the  thesis  is  concerned  with  the  derivation 
of  new  classes  of  sequences.  Some  new  theoretical  developments  are 
presented  on  the  cyclic  correlation  properties  of  sequences  containing 
the  Nth  complex  roots  of  unity  as  symbols.  These  sequences  are  related 
to  real  sequences  containing  phase  modulated  sinusoids  as  symbols. 

A  class  of  sequences  derived  by  the  interleaving  of  two  level 
binary  sequences  is  presented.  These  interleaved  sequences  are  shown 
to  exhibit  autocorrelation  and  cross -correlation  functions  with  inter¬ 
mediate  peaks.  The  locations  of  the  peaks  can  be  controlled  to  synthe¬ 
size  autocorrelation  i unctions  which  are  "almost"  two  level  or  autocor¬ 
relation  functions  with  peaks  of  different  magnitudes.  The  application 
of  these  sequences  to  synchronization  is  considered. 

Classes  of  N-ary  sequences  ( called  cyclically  orthogonal  se¬ 
quences),  which  exhibit  cross -correlation  functions  which  are  zero  for 
all  cyclic  shifts  are  derived. 

The  second  part  of  the  thesis  is  concerned  with  the  application 
of  the  sequences  derived  in  the  first  part,  as  well  as  other  known 
classes  of  sequences,  to  two  special  communication  problems. 
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The  first  application  is  a  binary  asynchronous  linear  multiplex 
system.  There  are  k  transmitter-receiver  pairs,  each  using  sequences 
(which  are  sunned  linearly)  as  carriers  for  binary  information.  Although 
a  particular  transmitter  is  synchronized  with  the  corresponding  receiver, 
the  other  transmitter-receiver  pairs  are  asynchronous  with  this  pair. 

The  average  probability  of  error  with  additive  white  Gaussian  noise  is 
evaluated  when  the  carriers  are  the  binary  cyclically  orthogonal  se¬ 
quences.  The  performance  is  also  determined  when  the  carriers  are  sinus¬ 
oids  of  different  periods, random  binary  sequences,  and  sequences  asso¬ 
ciated  with  Bose -Chaudhur i  Codes. 

The  second  application  is  an  N-ary  synchronous  "hard-limiting" 
multiplex  system.  For  this  application  there  are  also  k  transmitter- 
receiver  pairs,  although  the  respective  carrier  sequences  are  not  summed 
linearly,  but  "hard-limited"  prior  to  transmission.  It  is  shown  that 
the  optimum  set  of  carriers  are  the  cyclically  orthogonal  sequences. 

The  system  performance  is  determined  for  these  sequences.  This  system 
might  be  useful  for  satellite  repeaters  which  usually  employ  "hard- 
limiters  . " 
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CORRELATION  PROPERTIES  OF  MULTI-LEVEL  CYCLIC  SEQUENCES 


1.  INTRODUCTION 


Cyclic  or  periodic  sequences  are  gaining  increasing  Importance 
in  such  canmunications  applications  as  synchronization,  tracking  and 
ranging,  multiplex  carrier  systems,  and  signalling  over  a  continuous 
channel.1  This  thesis  is  concerned  with  the  derivation  of  sane  new 
classes  of  sequences,  and  the  application  of  these  sequences  to  two 
important  multiplex  problems. 

An  N-ary  sequence  S,  of  length  L,  is  defined  as  an  Ir  dimensional 
vector,  S  =  (Bo,a1;  ...,  aL«x),  with  symbols  a^,  i  *=  0,1,  ...,  L-l,  taken 
from  a  finite  alphabet,  usually  the  ring  of  integers  modulo  N,  where  N 
is  any  integer  greater  than  1.  For  practical  applications  the  symbols 
are  mapped  onto  quantities  a^ ,  for  which  the  ordinary  addition  and 
multiplication  operations  are  defined.  A  waveform  or  carrier,  f(t)  of 
time  duration  Ltx,  is  then  associated  with  the  mapped  sequence,  of  the 

0  <  t  <  tx 
tx  <  t  <  2tx 

(L-l)  t-i  <  t  <  LtA 


form 


(m) 


f(t) 


,(*) 


a (®) 

“L-l 


If  f(t+nLti)  =  0,  for  0  <  t  <  Lt*,  n  «  ±1,±2,  ...,  then  the  original  se¬ 
quence  is  referred  to  as  an  aperiodic  sequence.  If  f(t+Lti)  =  f(t),  for 
all  t,  then  S  is  referred  to  as  a  periodic  or  cyclic  sequence.  Only 
cyclic  sequences  are  considered  in  the  thesis. 
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Generally  the  terms  sequence,  mapped  sequence  and  carriers  have 
been  used  interchangeably  in  the  literature.  This  custom  will  be  con¬ 
tinued  in  this  thesis  when  there  is  no  possibility  of  ambiguity,  as  is 
the  case  with  binary  sequences.  However,  when  several  types  of  mappings 
are  possible  the  correct  term  will  be  used. 

1.1  Review  of  Pertinent  Prior  Work 

Several  types  of  periodic  sequences  have  been  studied  previously, 
with  the  research  divided  into  four  phases:  (a)  theoretical  derivation 
of  classes  of  sequences  exhibiting  particular  autocorrelation  and  cross - 
correlation  functions,  (b)  proofs  on  the  nonexistence  of  classes  of  se¬ 
quences  exhibiting  certain  correlation  properties,  (c)  methods  of 
generating  particular  sequences,  and  (d)  applications  of  sequences  to 
special  comnrunication  problems.  A  sampling  of  prior  work  in  each  of 
these  areas  is  now  briefly  described. 

a.  Theoretical  Derivation  of  Classes  of  Sequences 

Considerable  effort  has  been  devoted  to  the  derivation  of 

sequences  exhibiting  autocorrelation  functions  which  are  two-level,  i.e. 

the  autocorrelation  function  is  constant  (usually  -l/L)  for  all  nonzero 

integral  shifts.  The  best  known  of  the  two-level  sequences,  the  maximal 

length  sequences  with  binary  or  noribinary  components,  have  been  studied 
2  3  1* 

by  Singer,  Zierler,  and  Golcanb.  Other  two-level  binary  sequences 

which  have  been  synthesized  are  the  quadratic  residue  sequence, ^  twin 

6  7 

prime  sequence,  and  Hall  sequence.1 
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Perfect  N-ary  sequences  with  elements  mapped  onto  the  roots  of 
unity  (the  autocorrelation  function  is  zero  for  all  nonzero  shifts)  have 
been  studied  by  Heimiller8  and  Franck,  et  al. ^  Perfect  ternary  sequences 
with  elements  mapped  onto  +1,  -1,  0  have  been  found  by  Tompkins10  by  a 
computer  search.  Titsworth11  has  derived  a  class  of  two-level  sequences 
with  the  elements  contained  in  the  field  of  irrational  numbers. 
Titsworth^  has  also  derived  a  class  of  "almost"  two-level  binary  se¬ 
quences  for  some  lengths  not  covered  by  the  known  two-level  sequences. 

Extensive  searches  have  been  conducted  for  classes  of  sequences 
which  exhibit  low  cross -correlation.  Classes  of  orthogonal  binary  se¬ 
quences  have  been  derived  from  the  theory  of  Hadamard  matrices.1  The 
theory  is  presently  being  extended  to  roots  of  unity  sequences. ^  For 

nonorthogonal  sc  uences,  bounds  on  the  number  of  sequences  with  cross - 

ll* 

correlation  bounded  by  an  arbitrary  limit  has  been  derived.  Bounds  on 
the  number  of  sequences,  with  the  cyclic  cross -correlation  function, 
bounded  by  certain  limits  have  been  derived  by  Gilbert.1^ 

b.  Nonexistence  Proofs 

ff'uryn16  has  presented  proofs  on  the  nonexistence  of  perfect 
binary  sequences  beyond  length  4.  In  addition,  he  has  shown1^  that  two- 
level  N-ary  sequences  do  not  exist  for  many  lengths. 

c.  Generation  of  Binary  Sequences 

The  generation  of  many  classes  of  binary  sequences  by  linear 

1  17  18 

feedback  shift  registers  has  been  considered.  *  ’  Same  classes  of  se- 

19 

quences  can  be  generated  by  the  filtering  of  maximal  length  sequences. 
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d.  Applications  of  Sequences 

An  Important  application  of  sequence  theory  has  been  the 

use  of  the  orthogonal  binary  sequences  as  a  signalling  alphabet  for  the 

20 

continuous  Gaussian  Channel.  The  theory  has  been  extended  to  N-ary 

sequences.^  Sciffler^  has  considered  the  use  of  sequences  for  the 

synchronization  of  error-correcting  codes.  Binary  sequences  have  been 

25 

considered  as  carriers  for  binary  multiplex  systems.  ^ 

1.2  Brief  Review  of  New  Developments  Presented  in  the  Thesis 

In  Chapter  2  the  background  for  the  succeeding  chapters  is  pre¬ 
sented.  In  addition,  several  properties  of  nonbinary 
sequences  are  discussed.  A  theorem  on  the  shift  and  add  property  of 
nonbinary  maximal  length  sequences  is  derived,  from  which  the  autocorre¬ 
lation  function  of  ternary  (+1,  -1,  0)  maximal  length  sequences  is 
derived.  A  method  of  signalling  with  roots  of  unity  sequences  as 
carriers  is  outlined.  Information  is  transmitted  by  sending  a  sequence 
or  one  of  its  N-l  complements. 

In  Chapter  3  a  class  of  binary  sequences  with  useful  autocorre¬ 
lation  and  cross -correlation  properties  is  presented.  These  sequences 
are  derived  by  the  interleaving  of  tiro-level  sequences,  or  two-level  and 
complement  two-level  sequences.  It  is  shown  that  if  n  two-level  sequences 
are  interleaved,  the  autocorrelation  function  of  the  resultant  sequence 
has  n(n-l)  minor  peaks  of  height  -  -■  — — .  Fewer  peaks  of  larger 

amplitude  are  obtained  If  several  peaks  occur  at  the  same  location. 

Using  this  interleaving  procedure,  autocorrelation  functions  with  (n-l) 
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peaks  of  different  amplitude,  and  "almost"  tvo-level  autocorrelation 
functions  are  synthesized.  Classes  of  sequences  exhibiting  low  values 
of  cross -correlation  at  all  cyclic  shifts  are  synthesized  by  interleav¬ 
ing.  Autocorrelation  and  cross -correlation  functions  with  positive  and 
negative  peaks  are  synthesized  by  interleaving  the  complement  two-level 
sequence  along  with  the  uncomplemented  sequence.  Interleaved  maximal 
length  sequences  can  be  generated  by  the  nonlinear  filtering  of  a  single 
maximal  length  sequence.  A  number  of  examples  illustrating  this  tech¬ 
nique  are  presented.  The  application  of  the  interleaved  sequences  for 
synchronization  is  considered. 

In  Chapter  U  methods  are  presented  for  the  synthesis  of  N-ary 
sequences  (called  cyclically  orthogonal  sequences),  which  are  orthogonal 
for  all  cyclic  shifts.  All  of  the  binary  sequences,  derived  by  the  tech¬ 
niques  presented  in  this  chapter  are  of  different  least  period.  It  is 
conjectured  that  cyclically  orthogonal  binary  sequences  of  the  same  least 
period  do  not  exist. 

In  Chapter  5  the  application  of  cyclically  orthogonal  sequences 
as  carriers  for  an  asynchronous  linear  multiplex  system  is  considered. 

In  this  system  there  are  k  users  on  a  channel,  with  each  transmitter- 
receiver  pair  using  a  different  sequence  as  a  carrier.  Although  a  par¬ 
ticular  transmitter  is  synchronized  with  its  corresponding  receiver,  all 
of  the  other  transmitter-receiver  pairs  can  be  asynchronous  with  this 
particular  pair.  The  respective  carriers  are  summed  linearly  on  the 
channel.  The  density  function  of  the  interference  is  determined  and  the 
average  probability  of  error  is  evaluated  when  the  received  signals  are 
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corrupted  by  additive  white  Gaussian  noise.  Several  other  classes  of 
functions  are  considered  as  carriers  for  this  system. 

A  diiferent  type  of  multiple::  system  is  considered  in  Chapter  6. 
In  this  system  there  are  k  users  on  the  channel,  all  transmitting  synchro¬ 
nously  with  N-ary  sequences  as  carriers.  Information  is  transmitted  by 
sending  a  sequence  or  one  of  its  N-l  complements.  However,  the  carriers 
are  summed  and  then  passed  through  a  hard-limiter  prior  to  transmission. 
The  set  of  cyclically  orthogonal  sequences  are  shown  to  minimize  the 
average  probability  of  error. 
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2.  PRELIMINARY  DERIVATIONS 

In  this  chapter  some  basic  definitions,  mappings,  and  theorems, 
which  will  be  required  in  the  succeeding  chapters,  are  established. 
Several  types  of  sequences  are  described  and  the  autocorrelation  and 
cross -correlation  functions  of  the  mappings  of  these  sequences  are  de¬ 
fined.  The  complements  of  sequences  are  defined  and  the  application  to 
signalling  is  noted.  A  theorem  concerning  the  shift  and  add  property  of 
the  maximal  length  sequence  is  derived,  and  is  then  used  to  derive  the 
autocorrelation  function  of  a  mapped  ternary  (+1,  -1,  0)  maximal  length 
sequence . 


2.1  Basic  Definitions  and  Mappings 

An  N-ary  cyclic  sequence,  Sa,  of  length  L,  is  defined  in  this 
thesis  as  a  vector,  Sa  =  (a0,ax,  ...,  a]^) ,  where  a^,  i  =  0,  . ..,  L-l, 
can  assume  the  values  0,1,  ...,  N-l.  Sa  can  also  be  represented  as  a 
polynomial  Sa(x) , 

Sa(x)  =  a0  +  axx  +  agX2  +  ...  +  sl.jX1'"1  (2.1) 

where  x  is  an  indeterminate  and  xL  ~  1,  mod  L.  For  practical  applica¬ 
tions  the  symbols  of  Sa  Eire  mapped  onto  quantities  for  which  the  usual 
addition  and  multiplication  operations  are  defined.  Three  useful  map¬ 
pings  Eire  described  below. 

a.  Binary  Mapping 

For  N  =2  the  mapping  of  the  aj/s  is  0  -*  +1,  1  -  -1. 


b.  Complex  Roots  of  Unity  Mapping* 

The  N  possible  values  of  ai  are  mapped  onto  the  Nth  complex 
roots  of  unity,  0  -  e^0,  1  -  eJl^Jt/NJ  ...,  (N-l)  -  eJ(N-l)2*/n.  In  the 
determination  of  correlation  functions  this  mapping  will  be  shown  to  be 
equivalent  to  the  mapping  0  -  cos  tut,  1  -  cos(cut  +  2n/N) ,  ...,  (N-l)  -» 
cosCcut  +  (N-l)  2jt/ll] . 


c.  Mapping  Onto  Positive  and  Negative  Integers 

As  a  generalization  of  the  binary  mapping  the  N  values  of 
ai  can  be  mapped  onto  positive  and  negative  integers.  For  N  odd  0  -  0, 
1-1,  ...,  (N-l)/2  -  (N-l)/2,  (N+l)/2  -  -(N-l)/2,  ...,  (N-l)  -  -1.  The 
only  result  obtained  for  thi3  mapping  is  for  N  =  3* 


The  normalized  cross-correlation,  Pab(^)>  between  two  mapped 
sequences,  (Sa)^  and  (S^)^  delayed  by  X  digits,  where 


and 


(Sa)(u) 
(8b)  W 


=  (a(m)  aT(m)) 

\ao  >  •  •  •  9  aL-i' 


=  (b 


_(a) 

ai 


(m) 


Jmh 

bW 

.(m) 


; 


cure  the  mappings  of  aj.  and  bj_,  is  defined  as:++ 

^For  N  =  2  the  Nbb  complex  roots  of  unity  mapping  and  the  binary  mapping 
are  identical. 

++Tiie  general  time  cross -correlation  PabCf),  of  two  functions  Sa(t)  and 
Sb(t)  of  period  T  is  defined  as 

T 

Pab<T>  -  "f-  J  Sa(t)s£(t-r)  dt 
0 


For  sequences  the  integral  reduces  to  a  summation. 
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PabCU  =  -f  Lf  4m)  <4-1) *  (£-S) 

i=0 

where  the  subscripts  of  b^  are  taken  mod  L,  ard  (b^\)#  indicates  the 

f  \ 

complex  conjugate  of  b',  If  a  =  b,  then  (2.2)  reduces  to  the  autocor¬ 
relation  function  of  Sa. 

For  the  binary  and  roots  of  unity  mappings,  the  cross-correlation 
can  be  determined  from  the  po'.momia?.  reprssentr.nvons  of  Sa  and  Sb.  For 
the  binary  case  form  a  polynomial  Sc(x)  “  3a(x)  x^S^(x)  mod  (xL-l) 

(where  &  indicates  modulo  2  sum) .  Then  pa-0(  X)  is  seen  to  be : 

pab(X)  =  “X~  |_No-  of>  °’s  ^  Sc(::)  -  No.  of  l's  in  Sc(x)"|  (2.3) 


If  Sa  and  Sb  are  N-ary  sequences  mapped  onto  the  complex  roots 
of  unity,  then  form  the  polynomial  Sc(x)  ^  Sa(x)  -  x^Sb(x),  mod  (xL-l), 
(where  the  symbol  subtraction  is  mod  N) .  Then  pa^(  X)  is  seen  to  be 

N-l 

PabU)  =4"  I  [No*  0f  d'S  in  Sc(x)]  ej(2nd/N)  (2.4) 

d=0 


It  can  be  shown  that  the  mapping  onto  the  phase  modulated  cosine 
functions  is  related  to  the  mapping  onto  the  complex  roots  of  unity. 

The  cross -correlation  function,  pa^(X),  for  the  roots  of  unity 
mapping  is 


L-l 

-r  l  exp 

1=0 


2jc 

N 


(ai"bi-x)] 
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Under  the  mapping  onto  the  phase  modulated  cosine  functions  the  cross - 
correlation  function  becomes 

L-l  (i+l)ti 

Pab^)  -  l  J  cob  (art  +  aj)  •  cos  (cut  ♦  « 

i=0  itj- 

L-l 

’“ar  I  cos  [~: h*  (ai'bi-x)] 

i=0 

L-l 

+  ^  {sin  [(i+l)mti  +  ai+1]  •  cos  [( i+l)cuti  +  bi-x+i 

i=0 


-  sin  {iojti  +  a*]  •  cos  {loti  +  bj.x]  } 

where  ti  is  the  duration  of  one  time  slot.  The  second  summation  can  be 
dropped  if  mti  »  0,  and  thus  the  correlation  function  under  the  cosine 
function  mapping  is  proportional  to  the  real  part  of  the  correlation 
function  for  the  roots  of  unity  mapping. * 


♦The  analysis  for  the  phase  modulated  cosine  functions  assumes  that  the 
phases  of  the  carriers  of  the  two  functions  which  are  cross -correlated 
are  identical.  However,  if  the  phases  are  random,  then  the  typical 
correlation  term  becomes 

_A_  cos  [v  +  -Sf-  (H-bi-x)] 

where  cp  is  uniformly  distributed  between  0  and  2jt.  For  the  applications 
which  will  be  considered,  it  will  be  assumed  that  cp  is  zero. 
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It  can  be  shown  that  all  correlation  functions  of  mapped  se¬ 
quences  (real  and  complex)  are  linear  between  integral  cyclic  shifts, 
indicating  that  the  determination  of  the  correlation  function  at  in¬ 
tegral  cyclic  shifts  completely  specifies  the  correlation  function.* 

A  technique  for  binary  signalling  with  binary  sequences, 

(called  coherent  phase  shift  keying,  PSK)  is  to  transmit  the  mapping 
of  a  sequence  or  its  complement  (negative).  If  roots  of  unity  (or 
phase  modulated  cosine)  mapped  sequences  are  used  as  carriers,  then 
a  generalization  of  the  PSK  case  can  be  specified  by  defining  ccmple- 

f  Qm  I 

ments  of  N-ary  sequences.  The  r  complement  sequence,  of  Sa 

is  defined  as: 

si^  =  (a 0-r,  ax-r,  ...,  aL-1-r)  mod  N  (2.5) 

In  polynomial  representation  the  complement  polynomial,  S^^(x),  is 
seen  to  be: 

Sa^U)  =  Sa(x)  -  r(x)  (2.6) 

where  r(x)  =  r  +  rx  +  ...  +  rx^"1. 

The  cross -correlation  function  Paar(°)  between  a  mapped  root 
of  unity  sequence  and  its  rth  complement  is  Paar(0)  *  e~J(2«r/N)# 

From  a  recent  report  by  Reed  and  Scholtz,21  integral  expressions  for  the 
average  probability  of  error  for  Nth  roots  of  unity  complement  signal¬ 
ling  can  be  derived. 


*The  proof  for  binary  sequences  is  a  special  case  of  Theorem  5*1 
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The  complement  sequence  concept  will  be  used  in  Chapter  4  for 
the  synthesis  of  classes  of  N-ary  sequences  which  are  orthogonal  for 
all  cyclic  shifts.  The  complement  signalling  concept  will  be  discussed 
further  In  Chapter  6  when  the  N-ary  "hard"  limiting  multiplex  system  is 
analyzed. 

In  the  following  section  some  new  results  are  presented  on  the 
characteristics  of  maximal  length  p-nary  (N  =  p  where  p  is  a  prime  number) 
sequences . 


2.2  Maximal  Length  p-nary  Sequences 

The  best  known  of  the  p-nary  sequences  Is  the  maximal  length 

Pli 

p-nary  sequence.  Peterson  has  shown  that  the  polynomial  representation 
h^x)  of  a  maximal  length  sequence  can  be  derived  from  a  primitive  irre¬ 
ducible  polynomial,  gp(x),  over  GF(p).  If  gp(x)  is  of  degree  r,  then  the 

coefficients  of  h^x),  which  is  of  degree  pr-.l-r 

r 


bm(x)  =  — — r"\"  1 

rn(  x) 


(2.7) 


I* 

specify  the  first  p  -r  digits  of  the  maximal  length  sequences  of  length 
L  =  pr-l,  the  last  r-1  digits  being  zero.  The  sequence  corresponding  to 
^(x)  can  be  generated  by  an  r  stage  p  level  shift  register  with  feedback 
connections  prescribed  by  gp(x).  It  has  been  shown  that  in  a  maximal 
length  sequence,  the  field  elements,  1  through  p-1,  appear  exactly  pr“~ 

r-i 

times  and  the  field  element  0  appears  exactly  p  -1  times.  It  has  also 

2h  x* 

been  established  that  the  p  -1  r-tuples  formed  by  taking  r  successive 
digits  of  a  maximal  length  sequence,  aQ,a1,  ...,  a^j.;  ax,a2,  ...,  arJ 
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l 


t 


. ..;  apr_i-,ai>  •••>  ar..2f  are  tne  (pr‘l)  distinct  nonzero  r-tvples  over 
GP(p). 

k  ,  . 

It  is  well  known  that  maximal  length  [!l  -  p  =  2)  se¬ 

quences  exhibit  the  shift  and  modulo  2  add  (or  equivalently  the  shift 
and  modulo  2  subtract)  property.  That  is,  if  a  maximal  length  binary 
sequence  is  delayed  an  integral  number  of  digits  and  the  delayed  se¬ 
quence  is  added  (modulo  2)  to  the  original  sequence,  the  resultant  se¬ 
quence  will  itself  be  a  delayed  version  of  the  original  sequence. 

The  shift  and  add  property  may  be  written  as: 


,  CxT(X)H,(x)  \  =  1,  ....  L-l 

hm(x)©xhm(x)  =  j  mod(x^-l) 

(0(x)  X  =  0 


(2.8) 


It  can  be  shown  that  I(X)  is  unique  for  given  X,  and  for 
Xf  /  Xj  mod  L ,  I(XjL)  f-  I(Xj)  mod  L.  These  properties  of  I(X)  will  be 
used  in  Chapter  3  when  the  nonlinear  filtering  of  maximal  length  se¬ 
quences  is  discussed. 

Wolf  et  al  ^5  Using  a  set  of  tables  computed  by  Elspa^6  de 
scribed  a  method  for  determining  I(X)  as  a  function  of  the  primitive 
polynomial  gp(x),  over  GF(2) . 

The  following  theorem  describes  the  shift  and  modulo  p  add 
property  of  p-nary  maximal  length  sequences  (p  >  2). 


THEOREM  2.1  If  a  p-nary  maximal  length  sequence  (p  >  2)  of  length 
L  =  pr-l  is  delayed  an  integral  number  of  digits,  and  the 
delayed  sequence  is  added  (mod  p)  to  the  original  sequence, 
the  resulting  sequence  will  itself  be  a  delayed  version  of 
the  original  sequence,  except  for  a  shift  of  {■&-!)/ 2  =  L/2 
digits,  in  which  case  the  resulting  sequence  is  the  all.  zero 
sequence . 
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24  V 

Proof?  Peterson  has  shown  that  the  p  -1  nonzero  terms  In  the  ideal 
generated  by  hj^x)  in  the  algebra  of  polynomials  modulo  x^-1  correspond 
to  the  j^-l  cyclic  shifts  of  the  maximal  length  sequence.  Since  the 
ideal  is  an  additive  subgroup,  the  addition  (mod  p)  of  h^x)  and  a  shifted 
version  of  hm(x)  yields  a  unique  shifted  version  of  l^(x)  except  for  the 
summation 

hjnCx)  +  (p-lJhpjCx) 

which  will  yield  the  polynomial  0(x)  corresponding  to  the  all  zero 
sequence . 

It  will  now  be  shown  that  the  polynomial  (p-l)hm(x)  is  congru¬ 
ent  to  xL/2hm(x)  modulo  (xL-l);  the  latter  polynomial  corresponding  to 
a  cyclic  shift  of  L/2  digits.  We  know: 

gp(x)hm(x)  a  0  mod(xL-l)  (2.9) 

Thus  in  order  to  shew  that 

(P“l)  hjji(x)  =  xL/2^n(x)  mod(xL-l)  (2.10) 


it  is  sufficient  to  demonstrate  that 


(xL/2-p+l)  -  (xL/2+l)  a  0  mod  gp(x) 


(2.11) 


However,  for  (2.11)  to  be  true,  it  is  nocessary  to  show  that 


gp(sc)  l(xL/2+l) 


(2.12) 


o 

We  know  that  gp(x)  divides  (x  -1),  for  no  value  of  £  xess  than 
p37-!  =  L  since  gp(x)  is  primitive.  Thus 
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gp(x)|(xL/2tl)(xL/2-l) 


Since  Gp(x)  cannot  divide 


1, 


gp(.x)  |(xL/2+l) 


(2.13) 


establishing  the  theorem. 


The  following  corollary  to  Theorem  2.1  pertains  to  the  shift 
and  subtract  of  p-nary  maximal  length  sequences. 


COROLLARY  2.1  If  a  p-nary  maximal  length  sequence  is  delayed  an 
integral  number  of  digits  and  the  delayed  sequence  is  subtracted 
(mod  p)  from  the  original  sequence,  the  resultant  sequence  will 
itself  be  a  delayed  version  of  the  original  maximal  length  se¬ 
quence  for  all  nonzero  shifts. 


Proof:  Since  the  subtraction  from  hj^x)  of  any  term  in  the  ideal  gen¬ 
erated  by  l\n(x),  yields  a  unique  nonzero  term,  except  for  the  operation 
hjufx)  -  hjnCx) ,  the  shift  and  subtract  property  is  evident. 


From  Theorem  2.1,  Corollary  2.1,  and  the  properties  of  maximal 
length  sequences,  the  autocorrelation  functions  of  these  sequences, 
transformed  by  two  mappings,  are  now  presented. 


THEOREM  2.2  The  autocorrelation  function,  p^(X),  of  a  ternary  (p  =  3) 
maximal  length  sequence  transformed  by  the  mapping  0  -*  0, 

1  ~  1,  2  -  -1  is: 


,  (2-3r'^)/(3r-l) 
PA(X)  =j -(2.3r'1)/(3r-l) 
0 


X  =  0 
X  =  L/2 
X  /  0,  L/2 


(2.14) 
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The  proof  of  Theorem  2.2  is  presented  in  Appendix  A.  A  plot  of  Pa(X)  as 
given  by  Theorem  2.2  is  shown  in  Figure  2.1. 

The  autocorrelation  function  of  a  maximal  length  p-nary  se¬ 
quence  (all  p)  transformed  by  the  roots  of  unity  mapping  is 


PaM 


1  X  *  0 

-l/L  X  t  0 


(2.15) 


The  result  is  derived  using  Corollary  2.1  and  Equation  (2.4). 

A  plot  of  the  autocorrelation  is  shown  in  Figure  2.2.  This  type  of 
autocorrelation  function  is  known  as  a  two-levr  -©correlation  func¬ 
tion  fluid  the  maximal  length  sequence  is  known  as  a  two-level  sequence. 

Two-level  binary  sequences,  besides  the  maximal  length  se¬ 
quence,  have  been  derived  for  the  foilwing  values  of  L. 

a.  L  -  4  (Perfect  Sequence) 

b.  L  »  4n-l  is  prime  n  =  1,2,  ...  (Quadratic  Besidue  Sequence)^ 

g 

c.  L  =  p(p+2)  where  p  and  p+2  are  prime  (Twin  Prime  Sequence) 

d.  L  =  4n2+27  is  prime  ( Hall  Sequence)^ 

The  out-of -phase  autocorrelation  of  the  first  sequence  is  0, 
emd  it  is  the  only  known  binary  sequence  to  exhibit  this  property.  The 
out-of -phase  autocorrelation  function  of  the  remaining  sequences  is 
-l/L.  The  two-level  binary  sequences  with  out-of -phese  autocorrelation 
equal  to  -l/L  will  be  used  in  Chapter  5  to  derive,  by  interleaving, 
larger  classes  of  biimry  sequences  with  various  autocorrelation  and 
cross -correlation  properties. 
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AUTOCORRELATION  FUNCTION  OF  TERNARY  MAXIMAL  LENGTH 
SEQUENCE  TRANSFORMED  BY  MAPPING  O-^O,  l-M.2-^-1 


P 

A 


FIG  . 

AUTOCORRELATION  FUNCTION  OF  P  -  NARY  MAXIMAL  LENGTH 
SEQUENCE  TRANSFORMED  BY  ROOTS  OF  UNITY  MAPPING 


ie 


Two-level  perfect  N-ary  sequences,  mapped  onto  the  roots  of 
unity,  are  known®' ^  for  l  =  N2.  However,  except  for  the  maximal  length 
sequences,  the  general  theory  of  two-level  sequences  has  not  been  satis¬ 
factorily  generalized  to  N-ary  sequences.  For  N  =  5  a  two-level  root  of 

unity  mapped  sequence  (01221)  has  been  found  by  the  author  and  it  has 
27 

been  indicated  that  no  two-level  ternary  roots  of  unity  mapped  se¬ 
quences  exist  for  lengths  other  than  8*  9  up  to  25. 


3.  SEQUENCES  DERIVED  BY  THE  INTERLEAVING  OP  TWO-LEVEL  BINARY*  SEQUENCES 


3*1  Introduction 

In  Chapters  1  and  2  the  prior  research  in  the  derivation  of  two- 

12 

level  binary  sequences  vas  'U^inarized.  Titsvorth  has  derived  limited 

classes  of  "almost"  two-level  sequences  (the  autocorrelation  function  is 

close  to  zero  for  nonzero  shifts),  by  forming  either 

the  tensor  or  term-by-term  pioducts  of  the  known  two-level  sequences. 

20 

From  the  theory  of  cyclic  codes,  in  particular,  Bose -Chaudhuri  Codes, 
bounds  on  the  autocorrelation  functions  of  some  additional  sequences  can 
be  found.  Emphasis  has  been  placed  on  the  derivation  of  two-level  se¬ 
quences  because  they  provide  simple  and  accurate  methods  for  tracking, 
ranging  and  synchronization.1 

With  the  major  emphasis  placed  on  the  derivation  of  two -level 
and  "almost"  two-level  sequences,  the  general  problem  of  synthesizing 
sequences  exhibiting  arbitrary  realizable  autocorrelation  functions  has 
remained  unsolved. 

In  this  chapter  a  large  class  of  previously  unknown  sequences 
are  synthesized  by  the  arbitrary  n-fold  interleaving**  of  the  same  two- 
level  sequence.  It  will  be  shown  that  the  autocorrelation  function  of 
the  interleaved  sequence  is  similar  to  that  of  the  original  sequence 


*  The  discussion  to  follow  is  limited  to  the  interleaving  of  binary  se 
quences  since  many  more  binary  two -level  sequences  are  known  than 
general  N-ary  two-level  sequences.  However,  the  theory  presented  is 
applicable  to  the  interleaving  of  N-ary  sequences. 

**The  interleaving  method  has  been  used  previously^  to  derive  error- 
correcting  codes  with  multiple  burst  error  correction  capability. 


except  for  the  addition  of  n(n-l)  intermediate  minor  peaks  of  height 
(L-n+l)/nL.  A  number  of  these  peaks  can  occur  at  the  same  location  to 
produce  fever  peaks  of  larger  amplitude.  A  formula  for  calculating  the 
location  of  the  minor  peaks  is  derived.  Using  the  interleaving  tech¬ 
nique  ,  a  new  set  of  "almost1*  two-level  sequences  are  synthesized.  Also 
autocorrelation  functions  with  (n-1)  peaks  of  different  amplitudes  are 
presented.  The  cross -correlation  function  between  pairs  of  sequences, 
with  these  autocorrelation  properties,  is  considered.  The  theory  is 
extended  to  sequences  formed  by  the  interleaving  of  two-level  and  com¬ 
plement  two -level  sequences. 

An  n  fold  interleaved  sequence  can  be  generated  by  the  gating 
of  sequences  from  n  separate  shift  registers.  However,  it  is  shown  that 
a  class  of  sequences  formed  by  the  interleaving  of  maximal  length  se¬ 
quences  can  be  generated  by  a  relatively  simple  nonlinear  filtering 
technique.  Several  examples  of  sequences  derived  by  nonlinear  filtering 
are  presented. 

The  application  of  interleaved  sequences  for  synchronization  is 
discussed. 

3.2  Derivation  of  Peak  Locations 

It  was  shown  in  Chapter  2  that  the  cross -correlation  function 
of  two  transformed  binary  sequences  or  the  autocorrelation  function  of 
a  transformed  binary  sequence  could  be  determined  from  the  polynomial 
representations  of  the  sequences.  Throughout  the  remainder  of  this 
chapter  h(x)  and  x^x)  will  be  the  polynomial  representations  of  a 
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two-level  binary  sequence,  with  out-of -phase  autocorrelation  of  -l/L, 
and  a  two-level  binary  sequence  delayed  by  a  digits,  respectively. 

Hence,  the  polynomial  h’(x) ,  h'(x)  s  h(x)0x°^i(x)  =  h(x)  [l0xa],  mod 
x^-1,  has  L  zeros  if  a  =  0,  or  (L-l)/2  zeros  and  (L+l)/2  ones  it  a  f  0. 

Two  definitions  concerned  with  general  interleaved  sequences 
are  now  presented. 

Definition:  If  the  polynomial  Sa(x)  corresponds  to  the  sequence 
(a0,  a^j.)  then  the  polynomial  xna*vSa( x11) ,  where  n  is  a  positive 

integer,  a  =  0,  1,  . ..,  L-l,  v  =  0,1,  . ..,  n-1,  corresponds  to  the 
sequence 

(0,0,  ...,  0,ajj_a,  0,0,  ...,  0,aI<_clffl, 
v  .  n-1 

0,0,  ...,  0,0,  ...,  o) 

n-1  n-l-v 

Polynomials  must  now  be  taken  mod  (x  -1). 

Definition:  The  sequence  Sc  resulting  from  the  (2-fold)  interleaving  of 
sequence  Sa,  and  sequence  delayed  by  a  digits  is: 

Sc  =  (ao>^-ct>ai>ki-G>  •••$  aL-i^^L-i-a) • 

The  polynomial  Sc(x)  can  then  be  written  as 

Sc(x)  =  Sa(x2)  +  x2QfflSfc(x2)  mod  (x2^-l). 

The  synthesis  procedure  presented  is  to  interleave  a  two-level 
sequence  with  (n-1)  shifted  versions  of  the  same  two-level  sequence. 
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The  resultant  sequence,  (in  polynomial  representation),  denoted  as  Sr(x) 
is: 

n-1 

Sr(x)  =  h(xn)  £  xnav+y  f  mod  (x^-l)  (3.1) 

v=0 

where  the  ov's,  called  the  interleaving  constants,  can  assume  the  values* 
0,1,  ■  •  . ,  L-l,  and  =  0. 

The  procedure  for  deriving  the  autocorrelation  function  pr(n6+w), 
of  the  transformed  binary  sequence  corresponding  to  Sr(x),  for  an  arbi¬ 
trary  integral  delay  (n6+v),  6  =  0,1,  . ..,  L-l,  w  =  0,1,  , . .,  n-1  is: 

a.  Form  the  polynomial  Sp(x)  =  Sr(x)  0xI1^+vSr(x) ,  mod  (x^-l) 

b.  Apply  Equation  (2.3)  to  Sp(x)  [count  the  number  cf  O's  and 
l*s  in  Sp(x)  ]  replacing  l/L  by  l/nL. 

It  is  seen  that: 

xn6+vSrW  .  h(*n)  J  ^  (xnL.1}  (j.2) 

v=0 

Thus  the  polynomial  Sp(x)  is: 

Sp(x)  =  h(xn)  £  =ctn<Vfv)  [l© xn( 6+<^-'f+'r'“v*M«v) ]  mod  (xnL-X)  (J.J) 
v=0 

where 

(^  1  v  <  w 

hrv  =  , 

(0  v  >  w 

and  subscripts  are  taken  mod  n. 

*For  the  most  general  type  of  interleaved  sequence  Qo  /  0.  Hcwever,  as¬ 
suming  Qp  =  0  does  not  reduce  the  generality  of  the  autocorrelation 
functions. 
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The  exponents  of  the  terms  xn(  k+Qh-v+v-CVHVv)  determine  the 
value  of  the  autocorrelation  function  for  a  particular  delay  (n6+v).  If 
no  exponent  is  equal  to  zero  (for  a  particular  6  and  v)  then  the  sequence 
Sp  consists  of  n  interleaved  sequences  each  containing  (L-l'/2  zeros  and 
(L+l)/2  ones  providing  an  autocorrelation  value  of  pr(n6+v)  «=  -l/L.  How¬ 
ever,  if  a  exponents  r.re  equal  to  zero  (for  a  particular  set  of  6  and  w) 
then  Sp  consists  of  a  all  zero  sequences  plus  n-cr  sequences,  each  con¬ 
tain4  %  (L-l)/2  zeros  and  (l/KL)/2  ones,  providing  an  autocorrelation  value 
of  pr(n6+w)  =  (crL-n+o)/nL.  For  this  condition  the  autocorrelation  func¬ 
tion  is  said  to  have  a  "o"  order  minor  peak,  a  -  0,1,  ...,  n. 

From  (3*3)  it  can  be  seen  that  an  "n"  order  peak  occurs  at  a 
delay  of  0,  ( corresponding  to  unity  autocorrelation) ,  and  no  peaks  occur 
for  delays  of  n,2n,  . ..,  (L-l)n. 

From  (3*3)  it  is  determined  that  peaks  occur  for 

^w,v  s  “  °h-w+v  "  IVv  >  mcx*  ^  (3*^) 

where 

w  =  1,2,  . ..,  n-1 
v  =  0,1,  . ..,  n-1 

and 

ao  *  0 

In  order  to  identify  each  peak,  subscripts  have  been  placed  on  the  para¬ 
meter  6.  For  a  given  value  of  w  there  are  n  peaks  providing  a  total  of 
n(n-l)  single  order  peaks.  However,  a  number  of  single  peaks  can  occur 
at  the  same  point  to  produce  multiple  order  peaks. 

The  n(n-l)  peak  locations  cannot  be  arbitrarily  chosen.  Since 
there  are  (n-l)  independent  parameters  -  the  interleaving  constants 
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0^,0^,  . ..,  Q^1_1  -  only  (n-1)  peak  locations  are  con>^ -,.ilable . 

It  will  be  useful  also  to  lmow  the  cross -correlation  function 
between  two  interleaved  sequences.  For  two  interleaved  sequences, 

S*.l)(x)  =  h(xn)  £  x“4l)+v 
Mid  V“° 

S<.2)(x)  *  h(xn)  f  xno^£)+v  ,  mod  (xnL-l)  (5.5) 

v=0 

t 

cross -correlation  peaks,  6y^v,  occur  for: 

~  w+v  -  *  Jnod  L  (5*6) 

where  v,v  *  0,1,  ...,  n-1  and  a^l\  c/2^  are  the  interleaving  constants 
of  the  respective  sequences.  There  are  n2  single  order  cross-correlation 
peaks  (some  of  which  can  occur  at  the  same  point  to  yield  multiple  order  peaks). 

In  the  following  section  a  number  of  constraint  equations  on  the 
autocorrelation  peak  locations  are  derived  and  a  synthesis  technique  is 
indicated. 


3*3  Constraints  on  Autocorrelation  Peak  Locations 


The  following  constraint  equations  are  derived  from  (3*M* 

a.  If  all  of  the  peak  locations  for  a  given  value  of  w  are 
summed,  mod  L,  it  is  noted  that: 


n-1 


v=0 


w  ^  L  -  w  ,  mod  L 


(3.7) 
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b.  There  is  a  symmetry  relationship  which  the  peak  locations 
satisfy  which  ensures  that  pr(n6+w)  =  pr(-n6-w).  The  re¬ 
sulting  constraint  equation  is: 

&w,v-  -*‘n-w,v-w  “  1  >  nod  L  (3*6) 

c.  The  peak  locations  by  ,v>  v  ~  • •  •  t  n-1,  can  be  easily 

calculated  from  the  peak  locations  for  w  =  1  from  the 
equation: 

v 

-  £  6i,b  i  mo^  L  (3.9) 

bev-w+1 

Equations  (3*7)  and  (3-9)  can  be  used  to  synthesize  sequences 
with  various  autocorrelation  functions.  The  procedure  used  in  the  ex¬ 
amples  to  follow  is  to  prescribe  values  for  •••>  6i,n-2> 

utilize  (3-7)  to  calculate  6i#n-i,  and  utilize  (3-9)  to  calculate  the 
remaining  peak  locations.  Equation  (3»^)  is  then  used  to  determine  the 
interleaving  constants  from  &ifo>  •  &i,n-2* 

It  is  interesting  to  note  that  although  there  are  Ln_1  possible 
combinations  of  interleaving  constants  (i.e.,  Ln_1  different  interleaved 
sequences),  there  are  fewer  distinct  autocorrelation  structures.  If  the 
peak  locations  for  w  =  1  are  written  as  a  linear  array,  (61^0,61^1,  ..., 
£>i,n-i) ,  then  from  ( 3 •  9)  >  any  cyclic  shift  of  the  array  or  any  cyclic 
shift  of  the  reciprocal  array  (backwards  array),  will  yield  identical 
autocorrelation  functions  although  the  interleaved  sequences  will  in 
general  be  different.  The  exact  number  of  distinct  autocorrelation 
structures  is  the  number  of  unequivalent  n-tuples  in  the  ring  of  integers 
mod  L,  with  the  constraint  that  the  sum  of  the  elements  of  the  n-tuple  is 
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-1  mod  L.  Two  n-tuples  are  defined  as  unequivalent  if  and  only  if  one 
n-tuple  cannot  be  derived  from  the  other  by  a  cyclic  shift  or  a  cyclic 
shift  of  the  reciprocal.  The  number  of  such  unequivalent  n-tuples  has 
not  been  determined. 

In  the  following  two  sections  examples  will  be  presented  of 
two  types  of  autocorrelation  functions  which  can  be  derived  by  inter¬ 
leaving  two-level  sequences.  In  the  first  example  the  interleaving  con¬ 
stants  are  chosen  so  that  the  highest  autocorrelation  peak  is  "2"  order, 
providing  an  almoi  t  two-level  autocorrelation  function.  In  the  second 
example,  the  interleaving  constants  are  chosen  so  that  all  autocorrela¬ 
tion  peaks  are  of  different  heights. 

3.4  "Almost"  Two-Level  Autocorrelation  Function 

The  synthesis  procedure  will  be  to  choose  . ..,  6i,n-2 

derive  the  remaining  peaks  from  these  values.  Let  L  be  a  prime*  and  let 
n  =  L+- 1.  Thus  the  length  of  the  interleaved  sequence  will  be  L(IrKL). 
Choose  the  peak  locations  for  v  =  1  as: 


6,  ^  =  L-l,  0,  2  s  d,  3  2  2d,  ...,  ^  3  (L-l)d  ,  mod  L 

(3.10) 


i,o 


where  d  -  1,2,  ...,  L-l,  mod  L. 

Since  L  is  prime,  the  numbers  o,d,  ...,  (L-l)d  are  all  of  the 
distinct  numbers  mod  L.  It  is  easily  verified  that  the  peak  locations 
given  by  (3*10)  satisfy  the  constraint  (3*7)  since 


*Most  of  the  two-level  binary  sequences  are  of  prime  length.  The  quac- 
ratic  residue  and  Hall  sequences  are  always  of  prime  length  and  the 
maximal  length  sequences  of  length  2r-l  are  of  prime  length  if  r  is  a 
Mersenne  prime. 
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L  L-l 

£  6i,v  =  £  v  +  (L-l)  ~  -1  ,  mod  L 
v=0  v=0 

There  are  then  (L-l)  single  order  peaks  and  one  double  order  peak  for 
v  =  1. 

The  remaining  peak  locations  are  calculated  using  (3.9)*  In 
applying  (3*9)  consider  the  (L+l)  peak  locations  for  a  given  value  of 
w  divided  into  two  groups.  The  first  group  contains  those  values  of 
&wfy,  which  are  not  explicit  functions  of  6i^0*  P6^  locations  for 

this  group  will  assume  the  value 

v  v-1 

^v  2  £  6i,b  "  £  vd  "  v(wd)  -  (w+l)  ,  (3-11) 

b=v-w+l  v^-w  mod  L  # 

From  the  above  equation  it  is  noted  that  all  of  the  peak  loca¬ 
tions  in  the  first  group  are  unequal  since 

v(wd)  -  (w+l)  f  v'(wd)  -  (w+l)  mod  L 

for  v  f  v',  mod  L. 

The  second  group  of  peak  locations  contains  those  values  of 
6 w#v  which  are  functions  of  For  these  values  of  v  a  typical  peak 

location  is 

6w,v  -  v(wd)  — (w+1)  +  (L-l)  ,  mod  L  . 

The  peak  locations  in  this  second  group  io  ere  all  distinct.  However, 
a  peak  location  in  the  first  group  can  be  equal  to  a  peak  location  in 
the  second  group.  Hence,  the  autocorrelation  function  will  contain  zero. 
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first  and  second  order  peaks.  The  magnitudes  of  first  and  second  order 
peaks  are  0  and  +l/L  respectively.  The  autocorrelation  function  of  the 
sequence  derived  by  the  (L+l)  fold  interleaving  of  the  length  L  =  7  max¬ 
imal  length  sequence,  with  d  =  1,  is  shown  in  Figure  3-1- 

L  sequences  instead  of  (L+l)  can  be  interleaved  to  provide  auto¬ 
correlation  functions  bounded  by  double  order  peaks  by  choosing  the  peak 
locations  for  w  =  1  as  6lj0  =  L-l,  61;1  =  d,  =  2d,  ...,  6i,l-i 

=  (L-l)d. 

Fewer  than  L  sequences  can  be  interleaved  to  provide  sequences 
with  autocorrelation  functions  bounded  by  double  order  peaks.  If  a  se¬ 
quence  of  length  nL,  n  <  L,  is  to  be  synthesized,  a  synthesis  procedure 
would  be  to  find  a  set  of  n  numbers,  v,v+d,  ...,  v+(n~l)d,  which  sum  to  -1 
mod  L.  There  is  always  a  unique  number  v  for  any  d  and  n  such  that  this 
congruence  can  be  satisfied.  These  numbers  (in  order)  are  then  the  peak 
locations  0,6i  i,  . ..,  6i,n-i>  and-  it  can  shown  that  the  resultant 

autocorrelation  function  is  bounded  by  double  order  peaks.  It  has  been 
found  that  for  certain  values  of  n  the  autocorrelation  f ■unction  derived 
in  this  manner  is  bounded  by  single  order  peaks,  although  no  general 
result  to  this  effect  was  found. 

It  can  be  shown  that  the  cross -correlation  function,  for  all 
cyclic  shifts,  between  two  "almost"  two-level  interleaved  sequences,  for 
n  <  L  and  n  prime,  is  bounded  by  fourth  order  peaks.  The  proof  for 
n  =  L  has  been  presented  elsewhere the  proof  for  n  <  L  is  given  below. 

Consider  an  "almost"  two-level  interleaved  sequence  with  the 
peak  locations  for  w  =  1  chosen  as 


ALMOST  TWO -LEVEL  AUTOCORRELATION 
SYNTHESIZED  BY  8  FOLD  INTERLEAVING 
OF  LENGTH  7  MAXIMAL  LENGTH  SEOUENCE 

WITH  d  =  I 


62- 
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,(i)  (1)  ,(i)  (l) 


1,0 


=  v 


■'1,1 


*(l), 


.(i) 

l.n-i 


=  v 


(i) 


+  (n-l)d 


(i) 


Similarly  a  different  "almost"  two-level  sequence  will  have  peak  locations 


6^  for  w 
v,v 


«  1  chosen  as 


.(2)  (2)  .(2)  (2) 

1,0  *  v  >  6i,i  E  v 


♦  *(2>, 


M) 


• ) 


'i,n-i 


-1  =  v 


(2) 


+  (n-l)d 


(2) 


where  d^1),  d^2)  =  1,2,  ...,  L-l  with  d^1)  t  d(2),  mod  L,  v^1)  and  v(s)  are 
determined  from  (3. 7)  setting  w  =  1. 

The  cross -correlation  peak  locations,  6^v,  are  given  by  (3*6) 


K,,w  -  4lJ  -  -  ‘‘wv  >  mod  1 

where  the  J1)  's  and  o<a>'s  are  the  interleaving  constants  of  the  two 
sequences . 

From  (3.^)  it  is  determined  that 

41)  *  w*1’  +  ,  mod  L  , 

and 

!  (n-v+v)Jz>  +  (n-v*v)(n^v+l},d(a)„  _  ^  L 


Define  a  variable  813  ®w,v  =  6w,v  +  ^wv*  Thus  ®w,v  is  the  cross" 
correlation  peak  location  for  v  >  w,  and  the  peak  location  increased  by 
one  for  w  <  v.  Thus  g^v  is 


+  v(  v+lM^1) 
2 


(n-v+v)'/2)  -  ^ 


mod  L 

In  order  to  determine  the  maximum  multiplicity  of  peaks,  it  is 
necessary  to  determine  the  values  of  v'  /  v  mod  n  for  which  g^v  r-  vi. 
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mod  L.  The  procedure  is  to  set  g^v  3  g^  vt,  mo<*  ^  and  solve  for  the 
values  of  v*  mod  n  which  satisfy  the  congruence.  It  is  apparent  that  a 
quadratic  equation  in  v'  results,  and  if  r  is  prime  then  only  two  solu¬ 
tions  exist,  one  of  which  is  evidently  v*  s  v  mod  n.  Thus  any  value  of 
6^v  can  occur  with  maximum  multiplicity  of  two  for  a  given  value  of  w. 
However,  if  the  effect  of  Mwv  is  considered,  it  is  seen  that  the  first 
w  values  of  g^v  are  reduced  by  1  to  yield  the  cross -correlation  peak 
locations  6^y.  It  is  then  possible  that  four  values  of  y  will  be 
identical  (for  a  given  w) ,  providing  a  cross -correlation  function  bound¬ 
ed  by  fourth  order  peaks. 

The  class  of  sequences  just  described  has  "good"  autocorrelation 
and  cross -correlation  properties.  Thus  these  sequences  can  be  used  as 
a  set  of  code  words  for  a  coranuni  cation  system  with  the  transmitter 
asynchronous  with  the  receiver.  The  "almost"  two-level  autocorrelation 
property  can  provide  a  means  for  synchronization,  and  the  "good"  cross - 
correlation  property  for  all  cyclic  shifts  provides  a  low  probability 
of  deciding  on  the  wrong  code  word. 

As  a  final  remark  on  these  "almost"  two-level  interleaved  se¬ 
quence  it  should  be  noted  that  the  class  of  sequences  of  length  L2  forms 
a  set  of  good  error- correcting  codes,  if  the  code  consists  of  all  cyclic 
shifts  of  all  of  the  (L-l)  interleaved  sequences.  A  total  of  L  (L-l) 

l  2 

code  words  are  formed  with  minimum  distance,-*  dm-tn  =  t(L  -JL+4) .  Thus 
for  large  L  the  code  consists  of  (approximately)  L  code  words  of  length 

*The  minimum  distance  is  easily  evaluated  by  relating  the  maximum  value 
of  cyclic  cross -correlation  (equivalent  to  fourth  order  peak)  to  the 
Hamming  distance. 
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L  with  an  error-correction  capability  of  V~/h.  A  Bose-Chaudhuri  code, 

(the  best  of  the  known  linear  codes)  containing  L3  code  words  of  length 
L  provides  a  guaranteed  error-correction  capability  of  only  L  /2  log2  L. 
Thus  the  code  formed  from  the  interleaved  sequences  provides  protection 
greater  ihan  a  Bose-Chaudhuri  code  of  equal  length  with  an  equal  number  of 
words. 

In  the  following  section  a  class  of  sequences  with  autocorrelation 
functions  containing  intermediate  peaks  all  of  different  amplitude  will  be 
synthesized. 

3-5  Autocorrelation  Function  with  All  Peaks  of  Different  Amplitude 

The  synthesis  procedure  is  to  choose  the  peak  locations  b\r,v  for 
w  e.  l.  Choose  0lfO  -  6^3,  =  ...  =  61>n-2  =  where  v  can  assume  any 
value,  mod  L,  except  values  for  which  the  congruence  vn  5  L-l,  mod  L,  is 
satisfied,*  where  n  is  any  positive  integer.  The  choice  of  allowable 
values  of  v  provides  an  "(n-l)"  order  peak  In  the  autocorrelation  function 
at  (nv+l).  There  is  then  a  single  order  peak  at  (n61;n-i+l)  where 
6i,n-i  H  (L-l)  "  (n-l)v,  mod  L. 

The  remaining  peak  locations  are  calculated  from  ( 3  •  9)  •  It  is 
determined  that  6W^W._1  =  6V^W,  =  ...  =  5w,n-2  s  yV}  mod  L'  6w.n-i 
=  =  ...  =  ^wiW-2  ~  (L-l)  “  (n-w)v,  mod  L.  Hence  there  is  an  "(n-w)" 

order  peak  of  ( nwj+w)  and  a  "w"  order  peak  at  [n(L-l)  -  n(n-l)v  +  w], 
mod  L.  The  autocorrelation  function  thus  contains  (n-l)  peaks,  (exclud¬ 
ing  the  symmetric  peaks),  all  of  different  order.  If  n  <  L  and  nv  and  L 

^Interleaved  sequences  for  which  v  satisfies  the  congrusnce  vn  E  L-l, 
mod  L,  reduce  to  n  repeated  versions  of  a  cyclic  shift  of  the  original 
two-level  sequence. 
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are  relatively  prime,  the  different  order  peaks  correspond  to  distinct 
values  of  6.  This  condition  ensures  that  the  (n-1)  peaks  axe  not 
"bunched''  together. 

The  autocorrelation  function  of  the  interleaved  sequence  derived 
from  the  length  15  maximal  length  sequence  by  choosing  n  =  7>  v  =  1  is 
shown  in  Figure  3*2. 

It  will  now  be  shown  that  the  cross -correlation  function,  at 
all  cyclic  shifts,  between  two  "different"  peak  sequences,  derived  from 
different  values  of  v  for  L  a  prime  and  n  <  L,  is  bounded  by  double 
order  peaks. 

It  can  be  seen  that  two  different  "different"  peak  autocorrela¬ 
tion  functions  result  for  interleaving  constants  > 


(l) 

=  Ji) 

% 

“  0, 

and 

42) 

where 


=  v^,  a  2v^^,  . ..,  s  (n-l)v^  ,  mod  L  , 

=  0,  a[2^  -  v^,  ...,  s  (n-l)v^  ,  mod  L  , 


n  <  L 


vW, 


=  1, 


L-l 


,(i)  f  v(s) 


mod  L 


Cross -correlation  peaks,  6^v,  then  occur  for 

^w, v  =  ^  “  °n-w+v  "  ^wv  mo<*  ^ 

Defining  again  the  variable  gv,v  =  &w,v  +  Mw/  it  is  seen  that 

g -  (n-w+v)v^  ■  (w-n)v^  +  v(v^-v^)  mod  L 

Hence,  for  a  given  value  of  w,  g^v  f  g^vi,  mod  L  for  v  f  v*  mod  L. 


DIFFERENT  PEAI  AUTOCORRELATION  FUNCTION 
SYNTHESIZED  FROM  7  FOLD  INTERLEAVING  OF  LENGTH 
15  MAXIMAL  LENGTH  SEQUENCE  WITH  v  =  I 


} 
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Thus  the  values  of  g^v  are  all  distinct  for  a  given  w,  indicating,  after 
considering  the  Pyy  function,  the  maximum  multiplicity  of  cross -correlation 
peaks  is  2. 

The  application  of  a  simple  "different"  peak  sequence  for  syn¬ 
chronization  is  considered  in  Section  3.8. 

3.6  Complement  Sequence  Interleaving 

A  method  of  synthesizing  sequences  with  useful  autocorrelation 
properties  "by  the  n-fold  interleaving  of  the  same  two-level  sequence  was 
presented  in  the  previous  sections.  Additional  sequences  can  be  synthe¬ 
sized  by  proper  interleaving  of  a  two-level  sequence  along  with  the  com¬ 
plement  of  that  sequence  uch  sequences  will  exhibit  autocorrelation 
functions  with  positive  ana  negative  peaks. 

It  is  recalled  that  the  complement  sequence,  h^^(x),  of  a 
two-level  sequence,  h(x)  is,  h^^(x)  =  h(x)0l(x).  The  polynomial, 
h"(x),  h"(x)  £  hUX+Jx0!^*11^)  s  h^^(x)  ©x°^t(x),  mod  (xL-l),  has  L 
ones  If  a  =  0,  or  (l/+i)/2  zeros  and  (L-l)/2  ones  if  a  f  0.  The  poly¬ 
nomial  h,M(x),  h"'(x)  =  h^^(x)  ©x°^^^(x) ,  mod  (xL-l)  is  identical 
to  the  polynomial  h'(x)  defined  previously. 

The  synthesis  procedure  used  Is  the  n-fold  interleaving,  with 
arbitrary  shifts,  of  a  two-level  sequence  of  length  L  and  its  complement. 
The  resultant  sequence  (in  polynomial  representation),  denoted  as  Sr»(x) 
is: 

n-1 

Sr.(x)  s  ^  ,  mod  (x^-l)  ,  (3.11) 

v=0 

where  h-^x)  is  either  h(x)  or  h^l^(x). 
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In  order  to  find  the  autocorrelation  function,  it  is  necessary 
to  form  the  polynomial,  Spt(x), 

Spi(x)  5  Sr'(x)  0  x11^+vSr«(x)  ,  mod  (x^-l)  . 

Thus 

n"1  hv(xn)©xn(6+0h-w+'rfllv<,,wv)hnWxn)]  , 

mod  (x^-l)  (3.12) 

Peaks  will  occur  in  the  autocorrelation  function  for 

6  +  °b-wfv  "  +  V\rr  =  0  ,  mod  L  . 

However,  the  peak  can  be  either  positive  or  negative.  A  positive  peak 
occurs  if  hv(x)  and  hn_y+v(x)  are  both  uncomplemented  or  both  comple¬ 
mented  polynomials.  A  negative  peak  occurs  if  one  of  the  polynomials 
is  uncomplemented  and  the  other  complemented.  The  peak  locations  satis¬ 
fy  (3*7) t  (3.8)  and  (3*9).  No  peaks  occur  for  w  =  0,  (excluding  of 
course  v  =  0),  and  the  autocorrelation  function  for  these  delays  is 
-l/L. 

It  can  be  shown  if  a  weight  of  +1  is  assigned  to  a  positive 
peak  and  a  weight  of  -1  to  a  negative  peak,  then  the  total  weight  of 
the  peaks  of  a  sequence  interleaved  from  m  two-level  sequences  and  m 
complement  sequences  (m+n2  =  n)  is  (m-n 2) 2  -  n. 

In  Figure  3.2  an  autocorrelation  function  with  (n-l)  peaks  of 
different  amplitude  was  presented.  By  interleaving  complement  sequences, 
autocorrelation  functions  with  positive  and  negative  peaks,  all  of  dif¬ 
ferent  amplitude,  can  be  synthesized.  The  procedure  is  to  choose  hy.(x) 

=  h(x)  for  v  =  0,2,4,  ...,  and  hy(x)  =  h^^(x)  for  v  =  1,3,  .  .  .  The 


Sp.(x)  *  l  x“^v  ‘ 


v--0 
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peak  locations  fcr  w  =  1  are  again  chosen  as  6:  Q  =  ...  =  ^i,n-2  ~  v* 
Negative  peaks  occur  for  w  =  1,3,  . ..,  and  positive  peaks  for  v  =  2,4, 
....  A  positive  and  negative  "different"  peak  autocorrelation  func¬ 
tion  for  n  =  J,  L  =  13,  v  =  1  is  shewn  in  Figure  3*3* 

In  the  following  section  it  is  shown  that  sequences  generated 
by  the  interleaving  of  maximal  length  sequences  can  be  generated  by  the 
simple  nonlinear  filtering  of  a  maximal  length  sequence. 

3«7  Generation  of  Interleaved  Sequences  by  Nonlinear  Filtering 
19 

Raphael  has  described  a  technique  for  synthesizing  binary  se¬ 
quences  (which  are  recognized  as  interleaved  maximal  length  sequences), 
with  autocorrelation  functions  containing  minor  peaks,  by  the  nonlinear 
filtering  of  a  maximal  length  sequence.  The  filter*  (shown  in  Figure 
3*4)  forms  the  mod  2  sum  of  a  maximal  length  sequence,  of  length  L, 

(with  each  symbol  repeated  n  times,  n  =  2,3,  . ..),  and  this  sequence  de¬ 
layed  by  ny+y  digits,  y  =  1,2,  ...,  L-l,  y  =  1,2,  ...,  n-1.  Raphael 
found  the  autocorrelation  function  of  the  output  sequence  for  a  few  val¬ 
ues  of  L,  n,  7  and  y,  by  computer  simulation.  The  generalization  of 
these  results  are  presented  a a  an  example  at  the  end  of  this  section. 

First  a  multiple  nonlinear  filter  which  can  be  used  to  generate 
a  wider  class  of  interleaved  maximal  length  sequences  is  presented.  The 
filter,  as  shown  in  Figure  3»5>  consists  of  n  delay  elements,  delaying 
the  input  repeated  maximal  length  sequence  n70,  n7i+l,  . . . ,  n/n-j+n-l 

*The  filter  operation  is  nonlinear  because  the  mod  2  sum  of  binary  se¬ 
quences  is  isomorphic  to  the  product  of  mapped  binary  sequences. 


I  •  I  I  +  +  + 
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digits  respectively,  (yy  =  0,1,  . ..,  L-l,  y  =  0,1,  . ..,  n-1) ,  and  a  mod  2 
summer.  There  is  also  the  provision  that  particular  delay  elements  can 
be  switched  out  of  the  circuit  by  the  operation  of  switches  SW^SV^, 
SWn.1.  The  opening  of  switch  SVy  is  equivalent  to  adding,  (mod  2),  an 
all  zero  sequence  in  place  of  the  input  sequence  delayed  nyy+y  digits. 

The  analysis  of  the  multiple  nonlinear  filter  is  as  follows. 

The  input  maximal  length  sequence,  with  each  symbol  repeated 
n  times,  is  represented  as  the  polynomial  S^x), 

n-1 

SnU)  -  hnU11)  V  x°  ,  (J.1J) 

<7=0 


where  h^x)  is  the  polynomial  representation  of  a  maximal  length  sequence 
of  length  L. 

The  output  of  the  yth  delay  element,  £y(x) ,  in  polynomial  repre¬ 
sentation,  is 

n-1 

Sy(x)  E  xnyy+ya y(xn)  ^  x°  ,  mod  (x^-l) 

a=0 


where 


■  ay(xn)  l  ,  nod  (x"1-!) 

v=0 


= 


v  <  y 
v  >  y 


(5.1*0 


if  SWy  is  closed 
if  SVTy  is  open 


and 
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Then  the  polynomial  representation  of  the  output  of  the  mod  2 
summer,  S0^(x),  is 

Sqt(x)  =  S0(x)  ©SiCx)  ®  . ..  ©Sn.!(x)  , 

n-1  n-1 

s  1  x"  la  <cn(,,y+,lyv)ciy(xn)  >  mod  t*"1-1)  '  (5-15) 

v=0  y=C/^ 

where  £0  indicates  mod  2  summation. 

Since  the  mod  2  summation  of  n  arbitrarily  shifted  maximal 
length  sequences  yields  either  a  shifted  maximal  length  sequence  or  an 
all  zero  sequence,  S0t(x)  can  be  represented  as 

n-1 

S0t(x)  =  £  xnav+Vav(xn)  ,  (mod  xnL-l)  ,  (3.16) 

v*=0 

where  the  o^'s  and  a^x)  's  are  functions*  of  7y  and  ay(x)  ,  y  =  0,  . .., 
n-1.  Hence  the  output  sequence  is  an  interleaving  of  shifted  maximal 
length  sequences  and  all  zerc  sequences.**  There  are  L  possible  states 
of  the  delay  element  ( n/y+y) ,  corresponding  to  the  L  values  of  7y  with 
SUy  closed,  and  one  extra  state  corresponding  to  SWy  open,  providing  a 
total  of  (L+l)  states.  Hence  there  are  (L+l)n  states  of  the  nonlinear 
filter.  However,  there  are  also  (L+l)n  possible  n-fold  interleaved 


*av(x)  is  either  h^x)  or  0(x) 

**The  autocorrelation  functions  of  sequences  interleaved  from  binary  two- 
level  sequences  and  all  zero  sequences  were  not  considered  in  the  pre¬ 
vious  sections.  However,  it  can  be  shown  that  these  autocorrelation 
functions  have  similar  minor  peak  structures  as  those  sequences  consid¬ 
ered  previously  except  for  the  occurrence  of  an  "i"  order  peak  at  delays 
of  n,2n,  ...,  (L-l)n,  where  &  is  the  number  of  all  zero  sequences  inter¬ 
leaved. 
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sequences,  suggesting  the  possibility  that  every  interleaved  sequence 
can  be  generated  by  the  multiple  nonlinear  filter.  It  is  proved,  in 
Appendix  B,  that  this  is  indeed  the  cafe  for  n  =  2,  by  demonstrating 
that  each  filter  state  yields  a  different  interleaved  sequence  as  an 
output.  It  is  conjectured  that  the  result  is  true  for  all  n,  although 
the  method  of  proof  appears  involved  for  n  >  2. 

Two  examples  of  sequences  derived  by  nonlinear  filtering  are 
now  given.  In  the  first  example  the  general  solution  to  the  filter 
(Figure  3-4),  consisting  of  a  single  arbitrary  delay,  ny+y,  is  given. 

In  the  second  example  the  synthesis  of  a  certain  type  of  "different" 
peak  sequence  is  presented. 

Example  3*1  With  a  single  delay  element,  n7+y,  the  output  sequence, 
soT(x)>  111  polynomial  form  is 

S0t(x)  s  Sja(x)  ©x^^Sj^x)  ,  mod  (x^-l)  , 

or  y_1  n_1 

SoT(x)  sh(xn)  {xnI(?+1)  £  xb  +  xnI(7)  l  x*} 

b=o  c=y 

Thus  the  interleaving  constants,  Oy,  are 

O^q  =  0£j_  =  •••  =  Gy — ^  =  l(7^l) 

°y  -  Vi  =  •"  =an-i  "  I(7) 

Applying  the  theory  of  Section  3*2,  it  is  found  that  autocorre¬ 
lation  peaks  occur  at  shifts  of  n6+w,  with  the  following  peak  orders. 


(3.17) 


mod  (x1^-!)  . 


(3- 18) 


a. 


For  w  <  min[y,  n-y]  the.  peak  orders  are: 


"v"  order  peak  at  6  =  1(7+1)  -  1  -  1(7) 
"w"  order  peak  at  6  *  1(7)  -  1(7+1) 
"(n-2v)"  order  peak  at  &  =  0 


b.  For  y  <  v  <  n-y,  with  y  <  n/2,*  the  peak  orders  are: 


V  order  peak  at  6  =  1(7+1)  -  1  -  1(7) 
"y"  order  peak  at  6  *»  1(7)  -  1(7+1) 
"(n-y-w)"  order  peak  at  6  =  0 
"(y-v)"  order  peak  at  6  =  L  -  1 


c.  For  w  >  max  [n-y,  y]  the  peak  orders  are: 


"(n-w)"  order  peak  at  5  =  1(7+1)  -  1  -  1(7) 
"(n-w)H  order  peak  at  6  =  ?(y)  -  1(7+1) 
"(2w-n)"  order  peak  at  6  **  L  -  1 


It  is  noted  that  the  peak  locations  1(7)  -  1(7+1)  and  1(7+1) 
-1-1(7)  correspond  to  symmetrical  peaks.  A  plot  of  the  autocorrela¬ 
tion  function  of  SQj  is  shown  in  Figure  3.6  for  arbitrary  values  of  n, 

( n  even) ,  and  L,  and  three  values  of  y,  yx  =  1  <  y2  <  y3  =  n/2 .  The 
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case  y  =  2n  is  the  situation  investigated  by  Raphael. 


Example  3.2  Consider  the  multiple  nonlinear  filtering  network.  Figure 
3»5>  with  all  switches  closed  and  70  =  7X  =  •••  =  7n-x  =  0.  For  n  odd, 
S0t  becomes 

n-1  n-2 

S0t(x)  =  h(xn)  {  £  xb  +  x11  £  xb}  (3.19) 

b°0  b=l 

b  even  b  odd 


♦For  y  >  n/2,  then  y  is  replaced  by  n-y  in  all  relationships  in  case  (b) . 


AUTOCORRELATION  FUNCTION  OF  NONLINEAR 
FILTERED  REPEATED  MAXIMAL  LENGTH  SEOUENCE 
WITH  SINGLE  DELAY  ELEMENT  ny+y 
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Thus  the  Interleaving  constants,  Oy,  are 

°lo  *  "  •"  =  °n-l  ■  0 

Ofi  =  O3  =  . .  .  =  Qjj.g  =  1  . 

The  peak  locations,  6y^v,  for  w  =  1  are  then 

*i,o  =  *1,2  =  •••  =  6i,n-i  =  L_1 

^1,1  3  *1,3  =  “  ®i,n-2  3  1 

There  is  a  "[(n+l)/2]"  order  peak  at  a  delay  of  n(L-l)+l  and  a  "[(n-l)/2]" 
order  peak  at  (n+1).  Then  applying  (3-9 )j  the  peak  locations  &v,v  for 
arbitrary  v  <  n/2  are: 

v  odd:  6W  Q  =  6  x  =  ?>Vf2  -  . . .  =  6W^W_1  =  L-l 

*w,v  3  *w,v+2  =  ~  \r,n-2  3  1 

^w,w+i  *  ,w+3  =  *  •  •  “  &w,n-i  * 

There  is  then  a  "[(n+v)/2]"  order  peak  at  n(L-l)  h/  and  a  "[(n-w)/2]" 
order  peak  at  (n+w) 


v  even:  „ 

=  _ 

1 

«o 

II 

« 

11 

=  L-2 

v,o 

W,2 

W,W-2 

*v,i 

=  6w,3 

=  . . .  =  ^v,W-3 

=  0 

*V,V-1 

=  ^w,w 

=  *w,v+i  =  • • • 

s  6v,n 

There  Is  then  a  "(v/2)"  order  peak  at  n(L-2)+v,  and  a  "[n-v/2]"  order 
peak  at  v.  This  is  a  "different"  peak  autocorrelation,  although  the 
peaks  are  "bunched"  together.  The  plot,  for  n  =  7;  L  =  7  is  shown  in 
Figure  3*7* 


AUTOCORRELATION  FUNCTION  OF  SEOUENCE 
SYNTHESIZED  BY  NONLINEAR  FILTERING 


o 
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3.8  Synchronization  vith  Sequences 


In  order  to  detect,  by  correlation  techniques,  the  information 

carried  by  most  signals,  the  correlator  receiver  and  transmitter  must  be 

synchronized.  Synchronization  information  can  be  transmitted  on  the 
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channel  carrying  the  signed,  or  on  a  separate  channel.  For  separate 
channel  synchronization  a  two-level  sequence  (or  "almost"  two-level 
sequence) ,  whose  period  is  equal  to  the  duration  of  the  information 
signed,  can  be  transmitted  on  the  separate  channel.  The  two-level  prop¬ 
erty  of  the  sequence  provides  unique  synchronization  with  a  low  probabil¬ 
ity  of  error. 

A  two-level  sequence  can  also  be  used  for  synchronization  when 
the  synchronization  information  is  carried  on  the  channel  with  the  signal. 
In  this  case  the  two-level  sequence,  used  as  a  subcarrier,  is  phase 
modulated  by  the  signal. 

Several  techniques  for  utilizing  the  two-level  sequence  (lv 
either  the  separate  channel  or  same  channel  method)  are  possible,  assu^i 
ing  that  symbol  synchronization  is  known  for  the  sequence. 

a.  The  transmitted  sequence,  with  unknown  delay,  is  cross- 
correlated  with  all  cyclic  shifts  of  itself,  and  all  of  the  correlator 
outputs  are  stored.  The  shift  which  provides  the  largest  output  deter¬ 
mines  the  synchronization  decision.  The  probability  of  correct  synchro¬ 
nization  is  determined  by  noting  that  the  set  of  L  sequences  resulting 

from  L  cyclic  shifts  of  two-level  sequence,  form  a  set  of  simplex  codes. 
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The  error  probability  for  these  codes  has  been  derived. 


b.  The  transmitted  sequence  is  cross -correlated  with  successive 
cyclic  shifts  of  itself.  The  shift  which  provides  a  correlator  output  ex¬ 
ceeding  a  predetermined  threshold  specifies  the  synchronization  decision. 

Fewer  correlation  operations,  on  the  average,  are  lequired  for 
case  (b).  The  number  of  operations  can  be  reduced  further  by  using  inter¬ 
leaved  sequences  as  synchronization  carriers,  and  "locking-onto"  a  minor 
peak  or  the  peak  at  zero  delay.  The  average  probability  of  correct  syn¬ 
chronization  is  now  derived,  when  a  two-level  sequence  is  used  with  the 
technique  specified  in  (b).  and  when  two  classes  of  interleaved  sequences 
are  used. 

Consider  the  application  of  the  two-level  sequence,  of  length  L, 
with  the  autocorrelation  function  shown  in  Figure  3*8,  for  synchronization. 
Assume  the  transmitted  sequence  is  initially  delayed  by  b  digits,  b  =  0,1, 

. . . ,  L-l,  with  respect  to  a  reference  sequence .  The  first  shift  which 
yields  a  correlator  output  above  the  threshold  Mi  determines  the  synchro¬ 
nization.  If  the  average  sequence  power  is  S,  the  time  duration  of  the 
sequence  T,  and  additive  white  Gaussian  noise  with  zero  mean  and  two-sided 
power  spectral  density  N0/2  is  assume^  then  the  average  probability  of 
correct  decision,  PCl,  is: 


-  erf  ((Mr-1)  JSC)] 

(3-20) 


erf(y)  =  -fa 


dt 


(3.21) 


F  /  6  .  3.  a. 
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It  is  necessary  to  find  the  threshold  value,  Mi,  which  maximizes  ^Ci* 
For  high  values  of  ST/N0,  the  quantity  [l  +  erf  (Mi  v&T/N0)  ],  can  be 
approximated  by 


1  +  erf 


-  2  - 


Then  retaining  only  the  first  two  terms  in  the  expansion  of 


(3.22) 


2  - 


Mi 


/  itST 


exp 


V  N, 


the  expression  for  PCl  becomes 

“"ari1  - erf  C(Mi_1)  ftB  [L  -  Mp  G*  -f-)] 

J  ^  VnT  (3.23) 

For  given  values  of  L  and  ST/N0,  (3.23)  can  be  maximized  by  a  trial  and 
error  procedure. 

Nov  consider  the  application  of  two-fold  interleaved  sequences 
for  synchronization.  The  sequence  with  polynomial  representation, 

R(x)  *  hfx2)  +  xt'h(4l)(*a)  exhibits  the  autocorrelation  function  shown 
in  Figure  3*9.  If  the  delay  of  the  transmitted  sequence,  b,  is  ,  b  >  L', 
then  synchronization  can  be  achieved  by  "locking-onto"  the  negative  peak 
at  a  delay  of  L*.  If  b  <  L',  then  synchronization  is  achieved  by  "locking- 
onto"  the  positive  peak  at  zero  delay.  For  thresholds  iMe,  the  average 
probability  of  correct  synchronization,  PC2  Is  derived 


b=0 


(5.24) 
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Then  PC£  can  be  approximated  as,  for  large  signal- to-noise  ratios, 


If  L'  ^  L/2,  then  the  expression  for  *Cl  and  are  identical,  although 
the  use  of  the  interleaved  sequence  requires  only  half  the  number  of  de¬ 
lays  (on  the  average)  for  synchronization. 

As  a  final  example  consider  the  interleaved  sequence  R(x)  *=  h(x  ) 
+  x2a+1h^^(x2) .  If  a  is  properly  chosen,  the  autocorrelation  function 
shown  in  Figure  3*10  can  be  realized.  For  this  sequence  synchronization 
can  be  achieved  by  "locking -onto"  either  of  the  negative  minor  peaks. 
However,  when  a  particular  minor  peak  is  "locked- onto",  it  is  necessary 
to  perform  one  additional  operation  to  achieve  synchronization. 

The  synchronization  procedure,  referring  to  Figure  3.10,  is  as 
follows.  If  the  initial  delay  corresponds  to  region  III,  the  reference 
sequence  is  shifted  until  a  minor  negative  peak  is  detected.  Then  the 
reference  sequence  is  shifted  L"+l  digits,  and  the  correlator  output 
should  indicate  no  peak.  If  the  initial  delay  corresponds  to  region  II, 
the  reference  sequence  is  shifted  until  a  minor  peak  is  detected.  Then 
the  reference  sequence  is  shifted  L"+l  digits,  and  the  correlator  output 
should  indicate  a  large  positive  peak.  For  region  I  the  reference  se¬ 
quence  is  shifted  until  a  large  positive  peak  is  detected.  With  thresh¬ 
olds  Nb  and  -M4  the  average  probability  of  correct  synchronization  PC3 
for  large  signal -to-noise  ratios  is  approximated  by: 


AUTOCORRELATION  FUNCTION  WITH 
TWO  NEGATIVE  MINOR  PEAKS 
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The  number  of  operations  required  to  synchronize  with  this  se¬ 
quence  is  approximately  one -third  the  number  required  with  a  two-level 
sequence,  or  two- thirds  the  number  required  with  the  single  negative 
peak  sequence.  However,  a  significant  decrease  in  average  probability 
of  correct  synchronization  is  observed  with  the  reduction  in  synchroniza¬ 
tion  time.  For  ST/N0  =  100,  L  “  63,  L’  =  31,  L"  =  20,  the  values  of  PCl 
and  PC2,  maximized  over  Mi  and  M2  are  each  1-6*  10'11.  The  value  of 
?c3,  maximized  over  M3  and  M4  is  only  0.999* 

The  general  "different”  peak  sequence,  analyzed  in  Section  3*5 
can  be  used  to  further  reduce  synchronization  time  by  "locking-onto"  a 
minor  peak,  although  there  will  be  a  corresponding  decrease  in  the  proba¬ 
bility  of  correct  synchronization. 
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4.  N-ORTHOGONAL  AND  CYCLICALLY  N-OBTNOaONAL  SEQUENCES 
4.1  Introduction 

In  this  chapter  sequences  with  symbols  mapped  onto  the  com- 
plex  loots  of  unity,  which  are  N- orthogonal  and  cyclically  N-orthogonal 
will  be  derived. 

Definition;  A  class  of  N-ary  sequences  is  N-orthogom!  If  given  any  two 
sequences  In  the  set,  the  cross-correlation  function  at  zero  delay, 
between  the  mappings  of 

a.  the  sequences 

b.  one  of  the  sequences  and  any  complement  of  another  sequence 

c.  any  conqplement  of  one  sequence  and  any  complement  of  another 
sequence 

Is  zero. 

For  the  roots  of  unity  mapping  sequences  Sa  and  S^, 

*  (®o>°4 1  •  •  )  )  ®b  “  » •  •  *>^11-1 ) 

are  N-orthogonal  If 


1  r  e^§”^ai"r)  “  (^-t)] 

L  L 
1*0 


0  , 


(4.1) 


th  th 

where  r,  t  ■  0,  1,...,  N-l  to  yield  the  r  and  t  complements  of  Sft 
and  ^  respectively. 

Definition:  A  class  of  N-ary  sequences  is  cyclically  N-orthogonal  if  the 
class  is  N-orthogonal  for  all  cyclic  shifts  of  all  sequences. 
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Then  Sa  and  8,.  are  cyclically  N-orthogonal,  for  the  roots  of 
unity  mapping,  If 

.  A  I?ai(2Vn)«vr>  -  o  ,  (4.2) 

for  all  X,  r  and  t. 

In  this  chapter  techniques  will  be  presented  for  the  derivation 
of  sets  of  cyclically  N-orthogonal  sequences.  A  subset  of  the  set  of 
cyclically  N-orthogonal  sequences  form  a  basis  for  a  set  of  N-orthogonal 
sequences.  The  relationship  between  N-orthogonal  sequences  and  the 
generalization  of  Hadamard  matrices  is  indicated. 

All  of  the  techniques  for  the  synthesis  of  binary  cyclically 
orthogonal  sequences  produced  a  class  of  sequences  all  of  different  least 

period.  It  is  conjectured  that  cyclically  orthogonal  binary  sequences  of 

’  #  ) 

the  same  least  period  do  not  exist.  This  conjecture,  along  with  the 
Fourier  analysis  of  binary  sequences,  is  discussed. 

4.2  Derivation  of  Cyclically  N-orthogonal  Sequences 

In  this  section  some  theorems  will  be  derived  on  the  construction 
of  cyclically  N-orthogonal  sequences,  with  symbols  mapped  onto  the  roots 
of  unity.  Most  of  the  proofs  will  be  established  by  noting  the  symbol 
structures  of  a  sequence  which  is  the  mod  N  difference  between  the  two 
sequences  which  are  to  be  shown  cyclically  orthogonal.  Modulo  N  arith¬ 
metic  is  used  for  all  arithmetic  operations  with  the  symbols  of  sequences. 

A  lemma  is  now  presented  on  the  symbol  structure  of  N-ary  sequen¬ 


ces 
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Lenina  4.1  Given  any  N-ary  sequence,  S„  of  length  L,  then  the  sequence, 

Sa'  of  length  NL  which  consists  of  a  cooq>lement  of  Sa  followed 
by  all  other  distinct  complements  of  Sa,  and  the  sequence  8. 
Itself  (in  any  order)  has  each  symbol  appearing  exactly  L  times. 

Proof:  The  sequence  S’  can  be  denoted  as 

where  Sa^ij^  is  the  ij  complement  of  Sa  and  •  for  j/j*.  If  the  ith 
symbol  of  S  is  a  ,  then  the  i,  b*- i,  2L+i,  ...,  (N-l)lrH,  symbols  of  S 

&  X  S 

are  s^-Iq,  a^-ij,  . a^-l^,  which  are  all  distinct  modulo  N.  Hence 
each  symbol  appears  exactly  L  times  in  S'  and  the  lenma  is  true. 

A  method  of  constructing  two  cyclically  N- orthogonal  sequences 
is  given  by  the  following  theorem. 

THEOREM  4.1  Given  any  two  N-ary  sequences  Sa  and  S^,  of  length  L,  then 
the  sequences  Sa'  and  S^’,  of  length  NL, 

V  '  [Sa(V,  Sa<»i,  >,  .... 

Sb'  ‘  [V  Sb’  ®b]  ’ 
are  cyclically  N-orthogonal. 

Proof:  Form  the  sequence  Sc'  ■  (Sa')^r^  -  where 

xnI?^(8b,)(<lt)  represents  the  tth  conqplement  of  S^'  shifted  by  nl &l  digits. 

Noting  that  is  Just  the  sequence  repeated  N  times, 

S  1  reduces  to 
c 


sc'  -  [(8a-x\)(lV‘"t),  (S.-xi8b)('1il«-t), 

(vjt\)(<liH-i+r"t)] 

Thus  from  Lemma  4.1  each  symbol  in  8C*  appears  exactly  L  times  indicating 
that  Sa'  and  S^*  are  cyclically  N-orthogonal. 

0  .  * 

It  is  noted  that  the  sequences  S&’  and  8^’  defined  in  Theorem  4.1 
are  of  different  least  perl'  1.  Theorem  4.2  gives  a  procedure  for  construc¬ 
ting  non-binary  sequences  of  the  same  least  period  which  are  cyclically 

,t  t 

N-orthogonal. 


THEOREM  4.2  Given  the  sequence  8  ", 

fit 


-  [Sa#  Sa(<5l)' 


noting  that  the  complements  of  Sa  are  taken  in  a  specific  order  - 
not  the  arbitrary  order  of  Saf.  Then  the  sequences  8Sa"  and  3,Sa", 
p, 8**0, 1, ...,N-1,8^P’,  are  cyclically  N-orthogonal. 


Proof:  Form  the  sequence  Sc", 

Sc"  -  S(Sa")(<1r)  .  x^B’CSa")^*)  . 

Then  the  i,  L+i,  . ..,  (N-l)LM.,  (i  <  A)  symbols  of  S  "  are 

C 

P(ai.r)  -  e,[»1.r(s-n-1)-t]>  -8’[  ai_LnlI-r.-l)-t-l], 

Btaj-r+l)  -  B’ [a±_ Jt-(N-n-l)--t+l j  , 


respectively.  These  terms,  using  ring  of  integer  mod  N  operations,  reduce 
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to  X+0(p'«3)#  X+l(p*-0),  IM|  X+(N-l)(0*-0),  where  X^aj-r)  -  0'fa^j- 
(N-n-l)-tl.  In  order  to  complete  the  proof  it  is  necessary  to  show  that 


e 


kN-1 


2nX  I  2n(p'-p) 


0  . 


31 

However,  Van  der  Waerdan  shows  that  the  sun  of  all  powers  of  any  root 
of  unity  (except  l),  is  zero,  thus  establishing  the  theorem  for  1  <  X. 
The  proof  for  1  >  l  Is  similar  and  hence  Is  omitted. 


It  Is  of  Interest  to  note  which  of  the  sequences  08  ",  all  0,  are 

81 

of  the  same  least  period.  TWo  sequences  0S  "  and  0'S  "  are  of  the  same 

A  a: 

least  period  If  the  roots  of  unity  eu  N  and  e*rN  are  of  the  same  order. 

If  N«p  Is  prime  all  of  the  non-real  roots  are  of  the  same  order  and  thus 

all  of  the  sequences  08  ",  0-1,  ...,  p-1  are  of  the  same  least  period,  HL. 

If  N/  prime  then  only  the  values  of  0  corresponding  to  primitive  roots, 

yield  sequences  of  least  period  NL. 

From  Theorems  4.1  and  4.2  a  general  class  of  .(N«»l)k+1  cyclically 

k 

N« orthogonal  sequences*  of  length  H  L  can  be  derived  from  an  arbitrary 
sequence  Sa  of  length  L.  The  set  of  sequences  are:  the  all  zero  sequence 
plus  0jSa^,  pj-1,  ...,  N-l,  1-0,  ...,  k-1  where  Sa^  is 


Sa  -  ...,  •  ••*  Sa^i^,  ...,  ),  «•*,  Sa(^if.jL)J  , 

N*  N*  I^- 

repeated  N*"1"3,  times. 


♦The  least  period  of  all  of  the  sequences  is  not  N^L. 
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An  Important  special  set  of  (N-l)k+l  cyclically  N-orthogonal 
sequences  of  length  I!**  Is  derived  by  letting  S  «(0).  Then  the  sequence 

S_  reduces  to 

ai 

‘  i 

'  s*i  ■  [0>  0>  N~1>  •••»  frl>  •••'  •••*  -1] » 

h1  i?  i?- 

repeated  N^"1”1  times. 

It  Is  easily  shown  that  the  set  of  reciprocal  sequences,  are 
also  cyclically  N-orthogonal.  Examples  of  cyclically  N-c  .hogonal  sequen¬ 
ces  are  presented  below. 

» 

Example:  N«3,  k«2,  Sa*(0) 

-  (021021021) 

2Sq^  -  (012012012) 

S.  -  (000222111) 

2S^  *  (000111222) 

Example:  N-U,  k-1,  8a*(0) 

S.  »  (0321) 

2So  «  (0202) 

38^  -  (0123) 

It  can  be  proven  that  the  set  of  sequences,  Sai,  1*0,  1,  ..., 
k-1  are  a  basis  for  a  sot  of  N-orthogonal  sequences,  For  the  special  case 
of  8a*(0),  the  set  of  N^  N-orthogonal  sequences  of  length  N^  result.21 
For  N*2  this  set  reduces  to  the  set  of  2^  bi-orthogonal  sequences  known 
as  the  Walsh  functions.^1  The  set  of  binary  sequences  Sa^  will  hereafter 
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be  denoted  as  the  Basic  Walsh  functions. 

The  set  of  sequences  Sa^,  [with  Sa»(0)],  are  not  the  only  set  of 
cyclically  N-orthogonal  sequences  in  the  set  of  N-orthogonal  sequences. 
The  set  of  sequences  0^^  ',  i-1,  ...,  k-1. 


i 


0,  1,  ...,  N-l, 


»,.1"  -  1,  2,  N-l  (4.3) 


are  also  cyclically  N-orthogonal.*  It  can  be  shown  that  the  set  of  sequen¬ 
ces  S&  ‘  form  a  basis  for  the  set  of  N-orthogonal  sequences. 

It  is  shown  in  the  following  theorem  that  there  are  no  additional 
v 

sequences  in  the  set  of  N  N-orthogonal  sequences  which  are  cyclically 
orthogonal  with  each  of  the  sequences 


THEOREM  4.3  There  is  no  set  of  cyclically  N-orthogonal  sequences  in  the 
set  of  1 r  N-orthogonal  sequences  which  contains  more  than 
(N-l)k+l  sequences. 


Proof:  The  proof  will  demonstrate  that  there  are  no  sequences,  (except 
the  all  zero  sequence),  in  the  set  of  N^  N-orthogonal  sequences,  which  are 
cyclically  N-orthogonal  to  PjS^,  [with  Sa«(0)].  The  more  general  proof 
which  demonstrates  that  there  are  no  sequences  which  are  cyclically  orth¬ 
ogonal  to  each  of  the  sequences  PjS^1  is  similar  to  the  proof  which 
follows,  except  for  a  more  complicated  notation.  Hence,  this  more  general 


#lPor  i,  ...,  i-1  and  P^'-l,  the  set  S  1 


reduces  to  S_  . 

ai 


tiSZR*.  P9  il 
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proof  is  not  given. 

Consider  the  sequence  8^  generated  by  linear  combinations  of  the 
basis  elements  8^.  Then  Sd  is 


S, 


1 

J-l 


1j<V»  - 


lhe  sequence  S^,  which  has  the  same  least  period  as  >  iB  of>  the  fonn 


sd-  [So,  <Ser)(<lpijt),  (S (S0)(4(H-l)Bix)] 

N1^ 

repeated  times,  where  8^  is  a  sequence  which  has  the  property 

xu 

that  the  sum  of  the  sequence  symbols  (mapped  onto  the  N  complex  roots 

of  unity)  is  zero,  and  the.  first  symbol  of  S0  is  zero  also. 

Consider  the  cross  correlation  function  between  S,  and  (L  SRj 

d  ai i 

shifted  by  one  place.  This  is  done  by  considering  the  structure  of  the 
sequence  SQ  *  Sd  -  ~\>  which  is  found  to  be: 


Se  *  ,  Sty SaJ  +  J~8^^(N-l),  0, ,,,,  0,  P^(N-l),  0, ,,,, 

r?x 


O;  3^^(lT-l),  0,  ,  «,,  0] 

with  the  8,  (N-l)  term  occurring  in  the  1,  N^l, N^-N^l  places. 
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Thus  the  structure  of  Sft  is  identical  to  the  structure  of  repeated 

N^"ii  tines  except  that  the  first  term  of  S0  is  changed  from  0  to 

0 j  (N-l)/0,  Thus  the  sum  of  the  symbols  of  mapped  out  the  roots  of 

unity,  is  not  zero  and  S.  is  not  cyclically  orthogonal  to  S.  ,  indi- 

d  *X 

eating  that  there  are  no  sequences  cyclically  orthogonal  to  every  sequen¬ 
ces  in  the  set  0<Sa  . 

In  Section  4.3  the  relationship  between  N-orthogonal  sequences 
and  generalized  Hadamard  matrices  is  discussed. 


U.3  N-orthogonal  Sequences  and  Generalized  Hadamard  Matrices 

It  has  been  shown'1’  that  the  theory  of  bi-orthogonal  sequences 

(fc*2)  is  closely  related  to  the  theory  of  Hadamard  matrices.  A  Hadamard 

matrix  is  a  square  mxm  matrix  whose  elements  are  ones  and  minus  ones,  and 

whose  row  vectors  are  mutually  orthogonal  and  whose  column  vectors  are 

also  mutually  orthogonal.  The  dimension  of  a  Hadamard  matrix  must  be* 

bp*1,  2,  Ur,  r-1,  2,  ...  For  every  known  Hadamard  matrix  of  dimension  m 

there  is  a  set  of  m  binary  bi-orthogonal  sequences  of  length  m.  Many 

32 

techniques  for  the  construction  of  Hadamard  matrices  exist. 

The  concept  of  Hadamard  matrices  can  be  logically  generalized 

to  square  matrices  whose  elements  are  the  N^  roots  of  unity.  The  theory 

of  these  generalized  Hadamard  matrices  is  at  this  point  practically  non- 

13 

existent,  although  some  theorems  have  been  suggested  by  related  work. 
Throughout  the  remainder  of  this  section  the  elements  of  the  generalized 
Hadamard  matrlce  will  be  indicated  to  be  the  ring  of  integers  mod  N, 


♦Although  the  dimension  of  a  Hadamard  matrix  must  be  1,  2,  Ur,  Hadamard 
matrices  are  not  known  for  every  value  of  r. 
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although  It  la  understood  that  the  elements  should  be  the  mapping  of  these 
Integers  onto  the  complex  roots  of  unity. 

Generalized  Hadamard  matrices  of  dimension  exist  since  It  has 
been  shown  that  a  set  of  N-orthogonal  sequences  of  length  exist. 

As  an  example  consider  two  possible  matrices  [Hi  ],  [Ha  ],  for  N«p*3>  k*l. 

P 000“] 

[Hi]  -  012  [Hg]  “ 

l02lj 


If  generalized  Hadamard  matrices  [Hi]  of  dimension  mi,  and  [Ha] 
of  dimension  Oq,  are  known  then  a  generalized  Hadamard  matrix  [Ha]  of 
dimension  ma-mi'nigis  derived  by  forming  the  Kronecker  product  written 
[H3]«[Hi]x[Ha],  of  [Hi]  and  [Ha],  The  Kronecker  product  matrix  Is  formed  by 
substituting*  [Hg ]  for  each  zero  In  [Hi],  the  first  complement  of  [1^] 
for  each  one  in  [Hi],  ...,  the  (N-l)  complement  of  [Ha]  for  each  N-l  in 
[Hx].  As  an  example  consider  the  Kronecker  product  of  the  two  three  by 
three  matrices  indicated  above.  Then 


[H3]  = 


OOOlOOO'OOO 
021; 021 1 021 
012i012'012 


000; 222; 111 
021' 210 ! 102 
0121 201 !  120 


OOO! Ill '222 

0211102 1210 

012;i20'201 
•  1  •  -* 


> 


*The  proof  of  this  result  follows  immediately  from  the  proof  for  the 
ordinary  Hadamard  matrices^2  and  hence  is  omitted. 


1 
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where  the  Kronecker  product  formation  is  indicated  by  the  dotted  partition 
lines.  It  is  noted  that  the  Kronecker  product  of  matrices  of  dimension 
N^1  and  N^*  to  yield  ia  matrix  of  dimenstion  N^l+^*  does  not  yield  a  new 
class  of  N-orthogonal  sequences  since  a  class  of  N-orthogonal  sequences 
is  known  for  N*.  all  k. 

However,  new  classes  of  N-orthogonal  sequences  can  be  generated 
from  generalized  Hadamard  matrices  derived  from  certain  two-level  N-aiy 
sequences.  From  a  two- level  V  ry  sequence  of  length  L,  (for  the  roots 
of  unity  mapping)  with  out  of  phase  autocorrelation  -VL,  a  generalized 
Hadamard  matrix  of  dimension  L+l  is  derived  whose  rows  are  the  all  zero 
row,  plus  all  cyclic  shift  of  the  sequence  with  an  extra  zero  at  the 
beginning.  For  example  consider  the  matrix  [H*],  derived  from  the  two 
level  sequence  (01221), 


CHc] 


000000 

001221 

010122 

021012 

022101 

012210 


The  set  of  sequences  derived  by  forming  the  Kronecker  product  of  the 
above  matrix  with  the  matrices  of  dimension  N*  will  not  be  of  length  N*' . 

All  of  the  Kronecker  products  considered  previously  were  between 
two  matrices  each  with  elements  in  the  same  ring  (same  value  of  N).  How¬ 
ever,  some  additional  classes  of  N-orthogonal  sequences  can  be  derived  by 
forming  the  Kronecker  product  of  generalized  Hadamard  matrices  with  ele¬ 
ments  contained  in  the  ring  of  integers  mod  Nx,  Na  respectively.  The  re¬ 
sultant  matrix  will  contain  elements  in  the  ring  of  integers  mod 
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‘a 


[L.C.M,  (Nx,Na)].  For  example  the  generalized.  Hadamard  matrix  [Ife],  shown 
below,  with  elements  In  the  ring  of  Integers  mod  6,  is  derived  by  forming 
the  Kronecker  product  of  [H*],  (shown  previously)  and  [Hb]  with  elements 
in  GF(2), 


Thus  [^]  is 


[ft] 


ooooocxxxxxxf 

030303030303 

000044222244 

0303412525^1 

004400442222 

034103412525 

002244004422 

032541054125 

002222440044 

032525410341 

004422224400 

034125254103 


4.4  Cyclically  Orthogonal  Binary  Sequences  of  the  Same  Least  Period 
In  Section  2  of  this  chapter  som«i  techniques  were  presented  for 
the  synthesis  of  sets  of  cyclically  orthogonal  sequences.  The  sets  of 
binary  cyclically  orthogonal  sequences,  synthesized  by  these  techniques, 
contained  sequences  all  of  different  least  period.  For  example  the 
periods  of  the  three  Basic  Walsh  functions  of  length  8, 

»  01010101 
Sax  «  00110011 
=  00001111  , 

are  2,  4,  and  8  respectively.  Ihe  all  zero  sequence,  which  is  cyclically 
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orthogonal  to  each  of  the  Basic  Walsh  functions,  is  of  period  1. 

It  is  conjectured  thut  cyclically  orthogonal  binary  sequences 
of  the  same  period  do  not  exist.  This  conjecture  is  discussed  by  noting 
the  Fourier  series  of  mapped  binary  sequences.  The  following  theorem 
gives  necessary  and  sufficient  conditions  that  two  periodic  functions 
must  satisfy  if  they  are  to  be  cyclically  orthogonal. 

THEOREM  4.^  Two  periodic  functions  are  cyclically  orthogonal  if  and  only 
if  th  respective  Fourier  series  of  the  two  functions  have  no 
terms  in  common. 

Proof:  Let  of  period  (not  necessarily  least  period)  T,  be  repre¬ 

sented  in  the  complex  Fourier  series 


v**®  ,,  x  .v2rrt 
V/JT" 


C  e4  . 

y  > 


(M) 


Y=  -oo 


with  Fourier  coefficients  Cy.  Similarly  fv  '(t)  is  represented  in  the 
series 


f(S)(t)  =  l 


v2TTt 


0.5) 


yo  _  oo 


The  cross-correlation  function,  px 3 (th  ) between  these  two  periodic 
functions,  at  a  delay  of  ts  is  then 


Pia(Ta)=fr 

Jo 


00  r  _*  w2ttt£ 


W=  -oo 


(4.6) 


However,  Equation  (4.6)  is  equal  to  zero  for  all  t2  if  and  only  if  either 

(l)  r  (2)-]* 

Cy  or  Cw  is  zero  for  eveiy  w.  Since 

(a ) 

C  =0,  the  theorem  is  established. 

W 


WT  ■ » 


if  and  only  if 


From  Theorem  4.4  it  is  immediately  established  that  if  two 
binary  sequences  are  cyclically  orthogonal  then  at  least  one  of  the  sequen¬ 
ces  must  have  an  equal  number  of  l's  and  0's  since  otherwise  both  sequences 
would  possess  D.C,  (w=0)  components. 

It  is  clear  that  the  theorem,  as  presented,  does  not  prove  that 
cyclically  orthogonal  functions  of  the  same  least  period  do  not  exist. 

The  least  period  T  of  a  function  with  a  Fourier  series  gives  by  4.4,  with 

L 

non-zero  Fourier  coefficients  Cy  ,  Cy  ,  ...,  Cv.  ,  is 

1  2  K 

Tl  =  T/[G.C.D  (vi,  vs,  Vk)]  (4.7) 

(1).  \  (2),  v 

Thus  for  example  the  functions  f  (t)  =  sint,  f  (t)  =  sin  2t  +  sin  3t 
of  the  same  least  period  are  cyclically  orthogonal. 
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Consider  the  binary  sequence  of  length  L=2a  which  is  j  ones 
followed  by  a  zeros.  If  we  set  T=2ati  (where  tx  is  the  duration  of  one 
time  slot),  then  it  is  easily  verified  that  the  non- zero  Fourier  coeffic¬ 
ients  of  the  mapping  of  S  are  C , , ,  C,_.  ...  Then  any  sequence  cyclically 

O  ±A  Is 

orthogonal  to  Sa  can  only  have  even  Fourier  coefficients  indicating  that 

its  period  is  an  even  fraction  (£,  ...)  of  the  period  of  S^.  Thus 

there  is  no  sequence  cyclically  orthogonal  to  S  with  least  period  2otx 

a 

In  order  to  consider  more  general  sequences  it  is  necessary  to 
derive  the  Fourier  coefficient  of  an  arbitrary  mapped  binary  sequence. 
Consider  the  mapped  binary  sequence  3, 

S  =  [®oj  aij  •••>  aL-il  t  ai  =  *1  • 

Then  the  Fourier  coefficient  Cy  is 


Lti  /  \  iV2TTt 
cv  =  T  f^  (tje"^  Ltx 

LtlJo 


dt 


L-l  (b+i)t1 


1  b=0  0 


,  <V2TTt'|1. 

eX1l 


**  2ttJv  [  ^  V 


r  ■.  bv2n 

l.1  -rr*- 


b=0 


T  bv2n 

I  ab-ie  ] 

b=l 


(4.8) 
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Equation  (4.8)  can  be  simplified  to 


1 

2ttJv 


L-l 


-Jb’ 


2tt 


(4.9) 


In  attempting  to  prove  the  existence  or  non-existence  of  cyclically 
orthogonal  sequences  it  is  necessary  to  establish  only  ’  hether  a  particu¬ 
lar  Fourier  coefficient  exists.  Hie  following  remarks  on  the  existence  of 
Fourier  coefficients  are  noted. 

a.  A  particular  coefficient,  C  is  zero  if  and  only  if 

■  «  •  *  • 

L-l  bv2Tr 

^  abe  ^  =  0 

b=0 


b.  If  Cy=0,  then  Cv+mT.— 0 «  m=0,  1,  ... 

Thus  it  is  only  necessary  to  consider  the  values  of  V,  mod  L. 


c. 


Assume  Cv=0.  Then  from  (4.9)  it  i3  noted  that 


,,  2tt 

c  -  jbv— •  _ 

)  ab  e  L  - 0 


all  b  mod  L 
such  that 
ab=+1 


all  b  mod  L 
such  that 


ab="1 


(4.10) 


Substituting  the  indicated  values  of  a.  in  (4.10), 


i 


e 


-Jb' 


2tt 


1 


-Jb 


2n 


a  particular 
subset  Q  of 
the  integers 
mod  L 


the  subset,  R, 
of  integers 
mod  L  disjoint 
with  Q 


(4.11) 


I 


n 


However,  from  the  property  of  the  sum  of  all  powers  of  any  root  of  unity 
(see  proof  of  Theorem  4.2)  it  is  noted  that 


y 

all  beQ 


-Jb  v£H 
e  L  + 


l 

all  beR 


-jb  v£H 
e  l 


(4.12) 


But  (4.11)  and  (4.12)  contradict  each  other  unless  each  of  the  summations 
of  (4.10)  are  zero.  Thus  is  zero  if  and  only  if 

-jbv£l 

®be  I*  •  f 

all  b 
such  that 

v+1 

indicating  that  only  the  components  of  a  mapped  sequence  which  are  +1 
need  be  considered  in  determining  the  existence  of  a  particular  Fourier 
coefficient. 

Clearly  if  two  sequences,  Si  and  Sg  of  length  L  are  to  be  cyc¬ 
lically  orthogonal,  then  for  at  least  one  of  the  sequences  the  funda¬ 
mentals  Ci  (or  C.i),  must  be  zero.  However,  if  C.i  is  zero  then 


'l  eJbTi  =  0  (4.13) 

all  beQ 

If  Equation  (4.13)  is  to  be  satisfied  a  sum  of  a  subset,  Q,  of  the  L 
complex  roots  of  unity  must  be  zero.  For  the  following  situations 
Equation  ( 4 . 13 )  is  satisfied. 
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a.  Q  contains  all  of  the  Integers  mod  L.  Then  the  components 
of  mapped  sequence  associated  with  the  subset  Q  are  all  +1. 

b.  Let  I>pi*pa  where  px  and  pa  are  different  prime  numbers. 

Then  ( U .  13 )  is  satisfied  if  Q  contains  b,  b+px,  b+2px,..., 

b+(pa-l)px. 

The  least  period  of  the  sequence  derived  from  (a)  is  1.  The  least  period 

of  the  sequence  derived  from  (b)  is  px.  In  case  (c),  below,  a  technique 

for  synthesizing  sequences  with  least  period  L  for  which  C_x=Cx=0  is 


presented 

M 

c.  Let  ,  where  the  p^'s  i=l,  . ..,  M,  are  primes  not 

all  the  same.  Then  if  px  and  pm  are  different  prime  factors 
of  L,  let  Q  contain  the  integers  of  the  following  2  subsets 
(b,  b+L/Pjj  ,  b+2L/p£  ,  b-l/px);  (V,  b'+l/pm  ,  b'+2L/pm 
. ..,  b'-l/pja);  with  b/b',  and  b+ip^b'+mpm  ,  mod  L  for  all 
integers  l  and  m. 

As  an  example  of  the  procedure  outlined  in  (c),  let  L=2*2*3=12,  indicating 
px=3,  1^=2.  Then  let  b=0  and  b'=l;  the  integers  contained  in  Q  are  0,  kf 
8;  1,  7.  The  corresponding  sequence  (unmapped)  is  001101100111.  It  is 
easily  shown  that  the  above  sequence  is  of  least  period  L,  and  C_.x=Ci=0. 

.  It  appears  that  the  synthesis  techniques  outlined  in  (a),  (b) 
and  (c)  are  the  only  technique s^which  will  yield  (mapped)  binary  sequences 
which  do  not  possess  fundamentals.  The  technique  outlined  in  (c)  cannot 


be  used  if  L=2px  since  the  two  subsets  are  not  disjoint.  The  subset 

and  (b1,  b’+-i  ,  b'+r^  t  •••>  b’-i  )  respectively  contain 
Pi  2pi  Pi 

one  even  integer  and  one  odd  integer,  and  all  of  the  even  integers  or  all 

of  the  odd  integers  mod  L.  Hence,  it  appears  that  cyclically  orthogonal  bi 

nary  sequences  of  the  same  least  period  of  length  L=2p,  do  not  exist. 

M 

If  J>2  IT  p-j.  ,  where  M  >  1  (eg.  L-12)  then  the  synthesis  tech- 


*Additi'onal  subsets  of  integers',  bC^),  b(r)+L/pr  .  ...,  b(r)-L/pr  ,  can 
be  added  to  the  two  subsets  already  in  Q  (case  c)  provided  all  of  the 
subsets  remain  disjoint. 


73 


nique  outlined  in  (c)  will  yield  a  sequence,  Sx,  of  least  period  L,  for 

which  Ci=0.  However,  for  Si  it  can  "be  shown  that  the  Fourier  coefficients 

C^\  (v=mpm  m=l,  2,  ...,  and  v=£p^  ,  4=1,  2,  . ..)  will  not  be  zero.  Hence 

fa ) 

if  a  sequence  Sg  is  to  be  cyclically  orthogonal  to  S*  then  C'  '  must  be 

P  SL 

zero.  It  was  found  that  if  the  components  of  Sg  were  chosen  such  that 
(a) 

C  =0  and  the  sequence  was  of  least  period  L,  then  the  Fourier  coeffi- 

^  (® ) 

cients  of  Sg,  Cy  ,  v=rpr  (where  pr  is  a  prime  which  could  be  equal  to 

\  f1)  fa) 

P-J  were  not  zero.  Thus  the  Fourier  coefficients  C '  '  .  C'  '  are  not 

®  PjPr 

zero  indicating  that  Si  and  Sg  are  not  cyclically  orthogonal. 

As  a  final  remark  it  should  be  noted  that  classes  of  sequences, 
all  of  the  same  least  period,  with  "good"  cyclic  cross-correlation  prop¬ 
erties  can  be  derived  by  choosing  the  sequence  components  such  that  the 
mapped  sequences  have  few  Fourier  coefficients  in  common.  Utilizing  this 
technique  a  set  of  3  sequences  of  length  24  were  derived  which  exhibited 
a  cyclic  cross-correlation,  (between  any  pair  of  sequences)  bounded  by 
±1/4. 


These  sequences  are 


( 11110100101001113 1000000 ) 
( 101100101101001010101100 ) 
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5  BINARY  ASYNCHRONOUS  SIGN.iLLING 
5. 1  Introduction 

Significant  attention  is  presently  being  directed  towards  the 
derivation  of  techniques  for  the  simultaneous  transmission  of  many  bi¬ 
nary  messages  on  the  same  channel  Most  of  these  multiple-access  systems 

under  consideration  operate  as  follows  There  are  k  users  on  the  channel 

* 

with  each  transmit  ter- receiver  pair  assigned  a  different  carrier  In¬ 
formation  is  transmitted  by  sending  the  carrier  or  its  negative.  Correla¬ 
tion  detection  is  usually  specified,  with  the  polarity  of  the  correlator 
output  determining  the  decision  on  whether  the  carrier  or  its  negative 
was  transmitted  The  only  external  noise  is  assumed  to  be  additive  white 
Gaussian  noise  with  mean  zero  and  two- sided  power  spectral  density  Nq/2. 
Several  types  of  systems  with  different  constraints  can  be  specified. 

a.  Linear  Summation  of  Carriers,  All  Carriers  Synchronous 

The  carriers  from  the  different  transmitters  are  summed  on 
the  channel,  and  the  transmitters  are  all  synchronized  with  each  other, 
and  with  the  corresponding  receivers.  The  optimum  set  of  binary  carriers 
for  this  system  are  the  bi-ort'nogonal  carriers  specified  by  rows  of 
Hadamard  matrices  The  performance  is  specified  by  the  usual  binary 
PSK  probability  of  error  relationship  since  there  is  no  interference 
between  carriers. 


*In  this  chanter  a  carrier  will  be  the  unmodulated  signal  emanating 
from  the  transmitter.  The  carrier  can  be  the  mapping  of  a  binary 
sequence. 
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b«  Hard  Limiting  of  Carriers,  All  Carriers  Synchronous 

For  satellite  repeater  applications  the  linear  summation 

of  carriers  is  inefficient  because  of  the  nonlinear  characteristic  of 

travelling  wave  tubes  -  the  usual  satellite  transmitter.  The  linear 

sum  of  the  carriers,  for  this  application,  is  usually  hard  limited 

prior  to  transmission  over  the  channel.  This  problem  with  the  carriers 

33  * 

synchronized  has  been  studied  in  great  detail. 

c.  Linear  Summation  of  Carriers,  All  Carriers  Asynchronous 
Each  transmitter  is  synchronized  with  the  corresponding 

receiver  but  the  set  of  k  transmitter- receiver  pairs  are  asynchronous 
with  each  other.  The  carriers  are  summed  linearly  on  the  channels. 
Interference  among  carriers  then  arise  from  two  sources: 

1.  if  the  carriers  are  not  mutually  bi-orthogonal 
for  all  cyclic  shifts, 

2.  the  complementing  or  '  flipping"  of  an  asynchronous 
carrier. 

d.  Hard  Limiting  of  Carriers,  All  Carriers  Asynchronous 
This  is  a  combination  of  cases  (b)  and  (c). 

In  this  chapter  the  communications  system  described  under 
case  (c)  will  be  analyzed.  First,  some  preliminary  derivations  are 
presented  and  then  the  system  performance  is  analyzed  for  four  dif¬ 
ferent  choices  of  carriers  -  mapped  cyclically  bi-orthogonal  Basic 
V/alsh  functions,  sinusoids  of  different  periods,  mapped  randomly 

*  In  Chapter  6  some-  aspects  of  the  problem  of  the  hard  limiting 
of  mapped  N-ary  sequences  are  considered. 
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chosen  binary  sequence  and  mapped  cyclic  error- correcting  codes. 

Consider  the  asynchronous  linear  multiplexing  of  k  carriers, 
SQ(t),  Si(t),  . ..,  S]c_,(t),  of  period  T.  The  receiver  synchronously 
correlates  with  S^t)  and  the  phase  of  the  other  carriers  with  respect 
to  S^(t)  is  random.  (As  previously  stated  a  carrier  or  its  negative 
is  sent  by  each  transmitter.  )  A  block  diagram  of  the  system  is  shown 
in  Figure  5.1.  Then  the  received  signal,  r^(t)  at  the  input  to  the 
i^1  receiver  is 

k-1 

rt(t)  =  St(t)  +  [apS^t-T^.)  +  dpS^t-T,)  +  n(t)  (5.1) 

/.= O 

Jf./i 


where:  n(t)  is  white  Gaussian  noise  with  zero  mean  and  two-sided 
power  spectral  density  N  / 2,  is  a  uniformly  distributed  random 

variable  over  the  interval  0  <  Tr  <  T,  and  zero  elsewhere;  t.  and 
t,,  are  statistically  independent  for  1  4  ■  '  ( is  by  definition 
zero), 

and  /  i  1  with  equal  probability  over  the 

,  \  interval  0  <  t  <  t. 

a 1  .)  -  -  i 

«  ) 

'  /  0  elsewhere  , 


!±  1  with  equal  probability  over  the 

interval  t.  <  t  <  T 

- 

I  0  elsewhere, 


and  d 


(a') 

r 


are  statistically  independent  for 
.*  /  //  or  o/o', 


d^  =  ±  1  with  equal  probability. 


ASYNCHRONOUS  COMMUNICATION  SYSTEM  MODEL 
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If  dp/  dp*  then  it  will  be  said  that  carrier  S^,(t)  has  been  com¬ 
plemented  or  "flipped".  The  detection  procedure  is  to  cross -correlate 
r^(t)  with  S^(t)  with  the  polarity  of  the  correlator  output  determining 
the  decision  on  d^  The  output  of  the  i^  correlator  is  then 

T 

Zi  =  f  J  Si^^  dt 

o 

T  T 

-  i  J"  St(t)  n(t)  dt  [  dt  s“(t)  dt 

o  o 

k-1 

+  *  l  {  J  41)si(t)s/t-TPdt 

1=0  o 

H  i 

T 

+  |  d^3)  S.(t)  S^(t-T^)  dt  }  .  (5.2) 

T£ 

Equation  (5.2)  can  be  written  in  a  more  simplified  notation  as  follows: 

k-1 

Zi  =  pin  +  S  I  dJ l  pi£  ^5' 

f,=o 


where 


in 


T 

i  r 

T 


V 

O 


Si(t)  n(t)  dt  +  ^ 


T 

r  di  sa(t)  dt 

o 


is  Gaussian  distributed  with  means  ±  S  and  variance  NqS/2T,  where  S, 
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the  average  signal  power  (assumed  equal)  for  all  carriers,  is 


T 

5  f  sj(t)  at 


(S.4) 


T 

■ f  r  dt> 

o 

with  probability  l/2, 


»U)(Tt) 


T£ 

i  I  s±(t)  at 


i  i  st(t)  sf(t-Tx)  at  , 


with  probability  l/2 


(5.5) 


The  cross-correlatioi:  functions  pj^C^)  and  o^(t^)  are  ’tlie"un" 
flipped"  and  "flipped"  cross- correlation  functions  between  S^(t) 
and  S^(t)  delayed  by  seconds,  respectively,  and  d^  =  ±  1  with 
equal  probability. 

The  probability  density  function  of  P^(T£)  is  dependent  upon 
the  characteristics  of  the  carriers. 


If  carriers  S  (t)  and  S^(t)  are  cyclically  orthogonal,  then 

n(u)  _  o 

pi  l  ■  0 


1  K> 

*!  f  8  (t)  s/t-r,)  at 


(5.  5a) 


Bo 


Assuming  a  zero  threshold,  the  average  probability  of  binary 
error  in  the  1th  channel,  P(e^),  is 

o  oo 

p(e1)  =|  j'  p(zi|di=+l)  dz  +  [  p(zi|d1=-l)  (5.6) 

-CD  O 

where  p(ZjJd^=+l)  and  p(z^|di=-l)  are  the  conditional  pdfs  of 
fch 

the  output  of  i  correlator  given  that  di  =  +1  and  -1  respectively. 
These  conditional  density  functions  are  given  by  the  convolution  of 
the  pdfs  of  the  k  terms  of  (5.3). 

Most  of  the  carriers  to  be  considered  in  the  following  sec¬ 
tions  are  mapped  binary  sequences.  It  was  stated  in  Chapter  2  that 
the  cross -correlation  function  of  mapped  binary  sequences  between  in¬ 
tegral  delays  is  linear.  It  is  now  shown  that  the  "flipped"  cross¬ 
correlation  function  between  integral  delays  is  also  linear. 


THEOREM  5. 1  Consider  the  mapping  of  two  binary  sequences,  A  and 
B,  of  length  L,  where,  A  =  (a0,  aa ,  ...,  ar-i )  and 
B  =  (b0,  bj,. . .bL-!  )  and  aj_,  bj.  =  0,1.  (The  mapped 
elements  are  a\m'  and  b£m).  )  The  "flipped"  cross-correlation 
function,  t*+c)j  (X=0,1,  ...,  L-l,  0<  e  <  1)  is  linear 
between  integer  values  of  X  +  e» 


Proof:  The  "flipped  cross -correlation  function  at  a  delay  of  (X+r) 

Is 

X  X-l 

Jf)(,  )  .£  V  (»)  k(»)  ,  (l-c)  r  (m)  (m) 

°ab  U el  -  L  Z_aJ  Vuj  L  L  i  L->.+J+l 


j=o 


J=o 


■i  H 


L-l  L-l 

(m)  Jm)  (1-c)  \  „(m)  .  (ra) 


L-X+j  L 
J='+l  J=X 


)  aV‘/  b'  . 

L  J  L-X+ j+i 


p 


+ 


(x+l) 


Thus  the  "flipped1  cross-correlation  function  is  linear  between  in¬ 
tegral  values  of  X  +  c  . 


The  noise  or  interference  (in  the  determination  of  d^) 
due  to  the  presence  of  the  other  carriers  on  the  channel  is  related 
to  the  values  of  vX^)  and  (l^).  In  the  selection  of  binary 
sequences  for  the  binary  asynchronous  signalling  application  the 
goal  is  to  find  sequences  for  which  the  "unflipped"  and  "flipped" 
cross-correlation  functions  are  low  for  all  cyclic  shifts.  Since  it 
is  difficult  to  deterministically  derive  sequences  for  which  both 
types  of  crcss-correlation  functions  are  low,  the  basis  for  choosing 
the  sequences  analyzed  in  the  following  sections  was  that  the  value 
of  p^(\.)  be  lew  for  all  i,  l,  and  «. 

1  At  At 


5.  2  Signalling  With  Basic  Walsh  Functions 

In  this  section  the  performance  of  the  mapped  cyclically 
bi-orthogcnal  Basic  Walsh  functions,  as  carriers  for  the  asynchro¬ 
nous  signalling  system,  is  determined.  The  "flipped"  cross -correla¬ 
tion  function  is  determined  from  which  an  integral  expression  for  the 
probability  of  binary  error,  P(ei),  is  derived.  3y  application  of  the 
Central  Limit  Theorem  P^)  is  determined  for  k  -  oo,  and  all  i.  The 
probability  of  error  for  finite  k  is  found  by  asymp-i  otic  series  ap¬ 
proximation  and  computer  evaluation  of  the  probability  of  error  in¬ 


tegral. 


It  was  shown  in  Chapter  4  that  the  set  of  k  sequences, 


0,1, 


k  -  1  of  length  2^, 


where  SQ  is 
ai 


repeated 


k-  i  - 1 

2*  1  times. 


are  cyclically  bi- orthogonal.  Hence  for  the  mapping  of  these  se- 

f\X) 

quences  P^v^)  =  0,  i  /  l,  for  all  X^. 

It  is  now  necessary  to  derive  an  expression  for  the  "flipped" 

( f) 

cros 8 -correlation,  procedure  is 

a.  form  the  sequence  Sg  =  S&i0x^  S^p,  mod  (x^-1), 
where  Sap  is  a  sequence  which  is  the  last  X  ^ 

digits  of  Sajt  followed  by  the  first  L  -  \  ^  digits  of 
Sa^  complemented. 

b.  count  the  number  of  zeros  and  ones  in  the  first  in  the 

first  X  e  digits  of  S  . 

_ _ _  g 

Then  from  Equation  (5.5a),  the  "flipped"  cross-correlation  is 

o^)(X  )  —  ^°*  zeroB  in  ^irs't  i  digits  of  - 

^  ^  No.  of  ones  in  first  X,  digits  of  S  ]  (5.7) 

( f ) 

Ihe  derivation  of  (^4)  ^or  the  Basic  Walsh  Functions 
is  outlined  in  Appendix  C.  A  general  plot  of  ^ ^  ( X ^ )  for*  ^  >  i 

is  shown  in  Figure  5.  2.  It  is  thus  seen  that  the  maximum  value  of 

i+i-k 

the  "flipped"  cross-correlation  for  l  >  i  is  2  .  Thus,  the 

interference  a  particular  usercan  expect,  due  to  the  presence  of  the 

*  For  i  >  l  the  plot  is  similar  except  the  roles  of  i  and  l  are 
interchanged. 
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"flipping"  of  other  carriers  is  dependent  upon  the  sequence  he 
is  using.  The  longer  the  period  of  the  sequence  (i.e.  the 
higher  the  value  of  i),  the  larger  is  the  interference. 

Before  determining  the  average  probability  of  error  in  the 
presence  of  white  Gaussian  noise,  it  will  be  useful  to  determine  the 
noiseless  channel  performance  of  these  carriers. 

Without  the  external  additive  Gaussian  noise,  the  output  of 
the  i^*1  correlator  z^  is 

k-1 

»i  ■  dis  +  s  l  d£  eu<x£>  •  <5-8> 

1=  0 

Vi 

It  is  assumed  that  d.  =  -*1  and  the  minimum  value  of  z. ,  (z. )  .  , 

l  l  l  m3  n 

under  the  "worst" flipping  condicions  is  found.  An  error  will  occur 

( f ) 

if  (zj  <0.  From  the  maximum  value  of  o:.'(X„)  indicated  in 
'  i  min  i/.  f 

Figure  5.2,  (zi)min  is 

i-1 

■  S  -  2l'1'S  l  *  '  S  •  (5-9) 

f-o 

«  “  1  ^  ^  IcH*  i 

The  interference  terra,  2  "  S  £  2  +  (k-l-i)^  S,  is  maximized 

for  i  =  k  -  1  or  k  -  2.  Substituting  i  =  k  -  1  in  ( 5. 9 ) >  (zk-l^min 
is 

k-2 

I  ^  -  $  -  S(l-21-k)  >0 
.e=o 


(5.10) 


Hence,  an  error  cannot  occur,  even  under  the  "worst  flipping"  condi¬ 
tions  for  a  noiseless  channel. 

Now  the  average  probability  of  error  in  the  presence  of  vhite 
Gaussian  noise  will  be  determined.  In  Section  1  of  this  chapter,  it 
was  shown  that  the  average  probability  of  error,  P(e^)  when  correlating 
synchronously  with  the  carrier  associated  with  Sa^  is  given  by  Equation 
(5.6).  However,  since  all  of  the  noise  phenomena  are  clearly  symmetric, 
and  the  apriori  probabilities  that  d^  =  +1  and  d^  =  -1  (all  A)  are  the 
same,  the  expression  for  P(e^)  is  simplified  to 

o 

P(q)  -  [  p(zi|a1-+i)  a*t  (s.n) 

-00 


where  p(z^|d  =+l)  is  given  by  the  k  fold  convolution  of  the  proba¬ 
bility  density  functions,  p( Pin )  and  pfS^p^],  f  =0,  1,  ...,  k-1, 
But, 


p( ' in)  - 


mrjr 


exp  - 


(5. 12) 


p[Sp  ]  is  the  probability  density  function  of  S*  p  assuming  t 

If  f. 

is  uniformly  distributed  between  0  and  T,  and  Prob  ["flip"  in  Sa^]  = 
Prob  [no  "flip"  in  Sa^  ]  =  .  We  can  then  represent  p[S#p^]  as  the 

sum  of  the  two  density  functions 


p[s-pt.]  =  |p[S  •p^)(t(,) )  + 

Since  =  0 ,  all  r.t  then 

fj  L 


(5.  13) 
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<*i»  •  6'0i/,J 


(5. 14) 


The  density  function  of  the  "flipped"  cross- correlation  function  for 

i+i  ~  k 

uniformly  distributed,  is  uniformly  distributed  between  +2 
and  -2i+1”k.  Thus 


1  9-i-a+k 
S  * 


p[s-p[J4t  )]  - ) 

“  I  1  _-/!-»+ k 


I2 


-21+1-k  S  <  p  <21+,'kS 
elsewhere 

-2/>1_k  S  <  pu  <  2f+1'k  S 
elsewhere 


i  <  l 


i  >  l 


(5. 15) 


Substituting  (5.15)  and  (5.14)  into  (5.13)  the  expression  for 
p[S*pi/)]  becomes 


p(s.pu)  = 


lc  /  »  1  -i-3+k  „i+i-k  9i+i-k 

26(plP+S2  -2  S5Plje<2  S 


elsewhere 


i  <  i 
(5. 16) 


or 


p(s-pu)  = 


l*/n  W  1  ?-^-3+k  o/+i-k  q  <’.+i-k 

26(pU)+S2  '2  S±P1J!-2  S 


elsewhere 


i  >  f 


(5. 17) 


The  result  of  convolving  the  k  density  functions, 
p^A  =0,  1,  ...,  k-1,  A  f  i,  to  form  p(z|di=+l)  and  then  per¬ 
forming  the  integration  of  p(z|d^=+l)  to  derive  P(e^)  is  not  easily 
derived  in  closed  form  for  k  >  2.  The  derivation  of  P(eQ)  =  P(ej  ) 


°7 


for  k  =  2  by  the  convolution  method  is  presented  in  Appendix  D.  It 
is  shown  that 


(5.  lfl) 

A  convenient  integral  expression  for  P ( e ^ )  with  k  >  2  is  ob¬ 
tained  using  the  characteristic  functions  of  the  cross-correlation 
functions.  The  characteristic  function  M^(  jv)  of  a  variable  x  with 
pdf  p(x)  is  defined  as 

oo 

P 

M  (jv)  =  exp  (jvx)  p(x)  dx  (5.19) 

-00 


The  characteristic  function  M 

variable,  p.  is 
’  in 


Pin 


(j 


v) 


of  the  Gaussian  distributed 


(jv) 


=  exp  |^jvS  - 


,3  N3  s3l 

21®  J  ’ 


(5.20) 


Similarly  the  characteristic  function  Mpi^,(  jv)  of  the  random  variable 

S*  p.  4  is  for  i  <  L, 
it 


-i-3+k 

MPu^>  ■  2  +  ST— 


L 


exp  (jvS2i+1”k) 


exp  (-JvS2i+1"k)  j  .  (5.21) 


For  i  >  (*,  M.  (jv)  is  obtained  from  (5.21)  by  replacing  i  with  *. 
pi/ 

The  characteristic  function  M_  (jv)  of  the  correlator  out- 

zi 

put,  z ^  given  that  d^  =  +1,  is  given  by  the  product 


VJv)  '  X  Mpir,(Jv)l  K/.  (Jv)]k'1'1 


(5.  22) 


t<  i 


i> i 


Then  the  conditional  pdf  p(zitdi=+l)  is  derived  from  M^Jv)  by  the 
integral 

oo 

p(zi|di=+l)  =2^  f  \  (Jv)  exp  (-JvZi)  dv  .  (5.23) 

-oo 


Thus,  the  expression  for  the  probability  of  error,  when  correlating 
th 

with  the  i  carrier  becomes 

r  V2  ig  S31 

)  exp  [jvS - ^ — J  •  exp  (-Jvz^ 

[exp  (jvS2/+1"k)  -  exp  ( - JvS2£+1  "k)]j 

(jvS2i+1'k)  -  exp  (-JvS2i+1_k)]|  .  (5.24) 

This  integral  can  be  transformed  by  noting  some  simplifications  and 


P(e1)  = 


-OO 


f  4zi  h 


-oo 


i-i  r,  P-i*+2+k 
n  {!♦ 


L-o 


Jv 


{ 


i  2'i-2+k  r 

I  exp 


*  i2  +  Jv  L 


substitutions.  First  it  should  be  noted  that 


(5.25) 


O 


9 


o 

f  dzt  exp  (-jvz^  --  -  ~  +  TT6(v) 
-oo 


Then  the  exponentials  with  complex  arguments  can  be  replaced  by  sine 
and  cosine  functions,  the  odd  function  terms  of  the  integrand  deleted. 
Then,  P(e^)  reduces  to 


P(ei) 


CD 


|’[£i2j»  ]  [eXp  (.  £)]  V  ♦  1 


.<■.=0 


2f+1  -k0v 


] 


f  sin  21+1  ~kCv  ",k"i"1 


L  21+1  *kGv 


+  1 


dv 


(5.  26) 


where  G  =  /  —  . 

'  No 


It  is  clear  from  the  above  integral  that  P^)  is  a  function 
of  (a)  the  carrier  in  question,  i,  (b)  the  number  of  carriers  on  the 
channel,  k,  and  (o)  the  signal-to-noise  ratio  G2  =  ST/Nq.  It  is  also 
noted  that  =  P^..,). 

No  technique  was  found  for  evaluating  (5.26)  in  closed  form. 
Some  computer  derived  curves  which  give  the  probability  of  error  as  a 
function  of  the  parameters  G,  i,  and  k,  will  be  presented  later.  How¬ 
ever,  first  some  indications  of  the  performance  can  be  noted  by  evalu¬ 
ating  the  asymptotic  behavior  of  P(e^). 

a.  For  large  signal-to-noise  ratio  (G  -  oo) 

Gv 

First  evaluate,  lim  sin  —  considered  as  a  generalized  func- 
G-  co 


tion.  Thus 


90 


lim 
Q-  nr> 


aln  Gv 

V 


d_ 

dv 


lim 

G-«oo 


sin  Gv  \ 

J  -57- d(0v) 


d_ 

dv 


lim 

n-»oc 


nv 

r  sin  x 

.1 

o 


*  1 1? Ssn  v 


£ 

2 


6(v) 


Substituting  this  result  into  (5.26)  it  is  found  that 

00 

lim  P(e±)  =  |  -  [  •  6(v)  *  2k_1  dv  =  0 

G-  00  2 

o 

b.  For  small  signal- to-noise  ratio  (G-o) 

It  is  seen  that 


li  1 

G-o 


sin  Gv 
v 


=  0 


Thus 


lim  P(ei)  =  j 
Q~o 


c.  For  an  unbounded  number  of  carriers  on  the  channel  (k  —  00), 
when  we  are  correlating  with  carrier  associated  with  Sa^,  i  finite 


it  is  seen  that 


lim 
k—  00 


sin  21 


Gv 


+£-k 


=  1. 


Gv 


Thus 


lim  P(e  ) 
k-*oo 


1 

2 


1 

n 


00 


sin  Gv 


exp  (-  v2 /k)  dv 


V 

O 


(5 


.27) 
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The  integral  in  Equation  (5.27)  is  noted  to  be  —  erf  (G).^ 

Thus 

lim  P(ei)  =  erf  (g)  (5.28) 

k-oo 

Summarizing  the  results  derived  in  (a),  (b)  and  (c),  it  is  seen  that 
the  limiting  values  of  P(e^)  for  high  and  low  values  of  signal- tc- 
noise  ratio  (any  k,  i)  are  0  and  respectively.  As  the  number  of 
carriers  on  the  channel  becomes  unbounded,  then  P(ei),  (i  finite),  ap¬ 
proaches  the  probability  of  error  value  that  would  be  realized  if  there 
were  no  interfering  signals  on  the  channel.* 

For  large  values  of  k  and  relatively  small  values  of  i,  an 
asymptotic  series  approximation  to  the  integral  expression  (5.26)  can 
be  derived.  In  Appendix  E  an  asymptotic  series  approximation  to  P(eQ) 
which  yields  accurate  results  for  k  >  4  is  derived. 

_t  was  shown  that  the  probability  of  error  when  correlating 
with  the  ith  carrier  (i  finite)  approaches  the  PSK  error  probability 
as  k  -  oo.  Using  the  Central  Limit  Theorem,  the  limiting  probability 
of  error  curves  for  i=k-E,  k-oo  and  E  =  1,  2,  ,  can  be  derived. 

As  k  -  oo,  the  correlator  output  z ^  =  z^  E  becomes  Gaussian 

distributed  with  mean  m  and  variance  o3  , 

zi  zi 


=  m_ 


in 


+ 


(5.  29) 


#  This  is  referred  to  as  the  coherent  binary  PSK  error  probability. 
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where  m  ,  the  mean  of  p,  ,  is  +  S 
pin  ln 


and  nu  ,  the  mean  of  5*p,.  is  0 
P1jC 


Thus 


in  =  S, 


3  =  a3 


k-1 


in 


+  1  °s% 


i=0 


U 


where 


a  =  N  S/2T 

Pj  o  ' 
in 


and 


S-  p 


ijC 


S  _2i-ak-*a 
3 

S2  p3-S”2  k+i 
3 


i  <  L 


<  i 


Performing  the  summation  indicated  by  (5.30),  a2  =  a2  is  found 

Zi  Zk-E 

to  be  for  E  =  1,  2 


‘k-E 


oo 

V  .  C  ?2(k-c)-2k+x 
2T  3  i 
c=2 


oc 

_  2S3  , 1  ,  r  g-sd  V 
"  “  (IE>  l  2  +  zT 


d=o 


3  S2  N0S 

\  '  15  +  W  • 

k-E 


and  for  E  =  3,  4,  ... 


(5.30) 


(5.  31) 
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k-E 


N03 

2T 


T 


^(k-Ej-zk+i  + 


oo 


).  2 

c=E+l 


a(k-c)-2k+i 


3  M  £. 

\-E  =  2T  +  5 


[(E-l)  2'3E<1  M^HZ-815-1)]  • 


(5.  32) 


The  conditional  pdf,  pCz^Jd^+l),  for  unbounded  k  becomes 


lira  p(z1|di=+l)  = 
k-oo 


2 n  a 


exp  -  ] 


(5.  33) 


Hence  the  error  probability,  PCe^  =  P(ek  £),  for  k-oo  is 


Xlm  P(ek  )  = 
k-oo 


u 

f  p(\-Eidk-E=tl)  dZk-E  *  h  -  \  erf  [-JZ - ] 


-00 


'k-E 


(5.  34) 


i  ST 

Plot 6  of  lim  P(e^  ^),  as  a  function  of  G  =  ,1  >  for  E  =  1 

lc-  oo  ’  o 

(or  2),  3,  4,  5,  and  oo,  are  shown  in  Figure  5.3  for  ranges  of 

probability  of  error  down  to  10"5. 

The  curves  shown  in  Figure  5.  3  indicate  the  maximum  proba¬ 
bility  of  error,  P(e^  ^),  maximized  over  k  for  given  values  of  E  and 
G.  The  minimum  of  P(e  ),  (also  for  a  given  E  and  G)  occurs  for 

K-£j 

k  =  E.  Plots  of  P(eQ),  as  a  function  of  G,  for  E  =  k  =  2,  3,  4,  5,  oo, 
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PROBABILITY  OF  ERROR  VS  SIGNAL  TO  NOISE  RATIO 
(PER  BIT)  /NOISE  POWER  DENSITY  WHEN 
CORRELATING  WITH  S^  ,  S^^.5^,  . . 

WHEN  ALL  SEQUENCES  ARE  OF  LENGTH  Z*  AS  K-^CX> 
(ALL  CURVES  THEORETICALLY  DERIVED) 


95 


are  shown  in  Figure  5.4.  The  values  of  P(c0>  for  k  =  2  are  calcu¬ 
lated  from  Equation  (5.18).  The  values  for  k  =  4,  5  are  calculated 
from  the  asymptotic  series  approximation  derived  in  /\ppendix  E.  The 
values  for  k  =  3  were  obtained  by  numerical  evaluation  of  the  integral 
expression  (5.26)  on  the  PDP5  and  IBM  7094  computers. 

The  curves  shown  In  Figures  (5.3)  and  (5.4)  are  the  limiting 
curves  of  P(e^  g).  In  Figure  (5.5),  P(eic-i)  =  P(ek-2^  *s  plotted  f°r 
k  =  2,  3,  4,  5,  6,  oo.  The  curves  for  k  =  2,  oo  were  taken  from  Figures 
(5.4)  and  (5.3).  The  probability  of  error  values  for  the  remaining 
curves  were  obtained  by  numerical  evaluation  of  (5.26).  In  Figure 
(5.6),  P(e  )  is  plotted  for  k  =  3,  4,  5,  oo.  In  Figure  (5. 7)> 

k-3 

P(e  )  is  plotted  for  k  =  4,  oo. 

k-4 

The  error  probability  curves  not  included  in  these  plots, 

P(e,  _),  E  =  5,  6,  ...,  were  omitted  because  the  curves  (for  a  given 
E)  for  k  =  E  and  k  =  oo  (and  hence  the  curves  for  all  k  between  E  and 
oo)  were  indistinguishable.  From  the  given  plots,  it  is  possible  to 


determine  the  error  probability  curves  for  any  carrier 


)  with  any 


number  of  interfering  carriers  on  the  channel. 


It  was  noted  in  the  previous  analysis  that  the  highest  error 


probability  occurs  when  we  are  correlating  synchronously  with  carriers 


associated  with  S  or  Sn,  as  k  -*  oo.  If  the  carrier  associated 

ak-i  .  ak-2 

with  S  is  removed  from  the  channel  (providing  k-1  carriers  on  the 
ak-i 

channel),  then  the  maximum  error  probability  occuring  for  S  or 


S  will  be  significantly  lower.  Similarly,  if  S  and  S 
ak-3  ak-i  ak-n 

are  removed  (providing  k-2  carriers  on  the  channel),  then  the  highest 


LRROR  PROBABILITY 
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PROBABILITY  OF  ERROR  VS  SIGNAL  ENERGY 
(  PER  BIT  )  /  NOISE  POWER  DENSITY  WHEN 
CORRELATING  WITH*  SIGNAL  Sc  >  ( 10  1 0 10  10. . . .  ) 
CURVE  FOR  k  ■  2  DERIVED  FROM  CLOSED  FORM  SOLUTION 
CURVE  FOR  k  «  3  DERIVED  BY  COMPUTER 
CURVES  FOR  k  «  4.5  DERIVED  BY  ASYMPTOTIC  METHOD 
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PROBABILITY  OF  ERROR  VS  SIGNAL  ENERGY 
(PER  BIT) /NOISE  POWER  DENSITY  WHEN  CORRELATING 
WITH  Sh  4 

CURVE  FOR  k  *  4  ASYMPTOTICALLY  DERIVED 
FROM  SERIES 

CURVE  FOR  k  «  oo  THEORETICALLY  DERIVED 
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error  probability  occurs  for  or  S„  •  The  maximum  error 

K-3  “k-4 

probability  curves  when  there  are  k,  k-1,  k-2  carriers  of  length 
2  ,  with  k  -  co  o  1  the  channel  are  shown  in  Figure  5-8. 

For  fut>ire  comparison  of  the  performance  of  the  k  Basic 
Walsh  function  carriers  with  randomly  derived  mapped  binary  sequences 
for  asynchronous  signalling,  the  probability  of  error,  P  ,  averaged 
over  all  k  carriers  will  be  determined.  It  is  shown  in  Appendix  F  that 

r  =  !  -  f  erf  (G)  ,  (5.34) 

for  all  values  of  G2  <  9- 


5.  3  Asynchronous  Binary  Signalling  With  Sinusoids 

In  the  previous  section  the  performance  of  the  k  Basic  Walsh 
functions  cf  length  2  ,  as  the  carriers  for  the  binary  asynchronous  com¬ 
munications  system,  was  analyzed.  In  this  section  another  set  of  cy¬ 
clically  orthogonal  signals  -  the  set  of  sine  functions  whose  fre¬ 
quencies  are  all  different  and  all  multiples  of  some  fundamental  fre¬ 
quency  -  will  be  analyzod  for  this  application,  and  the  asymptotic 
performance  will  be  compared  with  the  asymptotic  performance  of  the 
Basic  Walsh  functions. 

Consider  the  set  of  sinusoids  S^t)  of  period  T, 


where 


S  (t)  =  /is  sin  i  co  t  ,  0  <  t  <  T  (5.35) 

(jo  =  2n/T 
i  =  1,  2,  ...  , 

S  is  the  average  signal  power  assumed  equal 
for  all  i. 


* 
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PROBABILITY  OF  ERROR  VS  SIGNAL  TO  NOISE  RATIO 
(PER  BIT)  /  NOISE  POWER  ,  MAX  I  MUM  PROBABILITY  OF 
ERROR  CURVES  FOR  k,  k-l,  k-2  SEQUENCES 
ON  CHANNEL  ALL  OF  LENGTH  2*  ,  k-*-®® 

(ALL  CURVES  THEORETICALLY  DERIVED) 
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A  set  of  these  sinusoids  are  considered  as  carriers.  It  is 
clear  that  S^(t)  and  S„(t),  i  /  /  are  cyclically  bi- orthogonal.  The 
"flipped"  cross-correlation  function,  p^'(t„)  between  S^t)  and 
S  (t)  at  a  delay  of  t  ,  0  <  t,  <  T,  is  derived  as  follows.  Since 
Si(t)  and  S^(t)  are  cyclically  orthogonal  the  "flipped"  cross-  correla¬ 
tion  is,  from  (5.5a), 


»i?<v 


2 

T 


C 

O 


S^t)  SJ&(t- 


Tx} 


dt 


TX 

=  ~  (sin  i  cjt][sin  X  u)(t-x^,  )]  dt  .  (5.36) 

o 


Performing  the  indicated  integration,  it  is  seen 


sin[(i-X^u>t  +  X  o>  t  ] 
2u)(  i-X) 


8in[(i+X)urt;  -  i  o  xj  ..  ‘X 


Mi+X) 


I 


or 

(  f )  2S 

pi/(0  =  -  [i  sin  X  cu  t  -  X  sin  i  o>  t  ]  (5.37) 

U  1  k(  i3-X2)  1  X 

Pbr  noiseless  channel  consideration,  we  are  interested  in  the 
maximum  interference  resulting  from  the  presence  of  other  carriers  on 
the  channel,  when  correlating  synchronously  with  S^t).  It  is  thus 
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necessary  to  find  the  maximum  value  of  the  flipped  cross-correlation, 
maximized  over  t  .  Differentiating  (5.37),  the  following  expression 
is  obtained: 

d  P$J?(T  )  2S  w  1  t 

— aT~"~  *  «(1°-C°)  tcos  Tx  -  c°6  1  m  TJ  (5-38) 


Thus  local  extremum  points  are  found  for 


cos  j!  oi  t,  -  cos  i  0)  t  =  0  , 

L  & 


or  equivalently. 


2  sin  -|{i+£)^  *sin  -|(,*-i)cu  t  =  0 


(5.39) 


Equation  (5.39)  is  satisfied  for 


Tf  =  i+j5 


2nut 

"  h  =  FT  ’ 


(5.  4o) 


m  =  Oj  •  •  • 


It  is  still  necessary  to  determine  vhich  value  of  x  ,  or  equivalently 

Xf 

~  t^,  given  by  (5.40),  absolutely  maximizes  p^T^).  This  can  be  de¬ 
rived  by  substituting  the  values  indicated  by  (5.  40)  into  (5.37). 


Case  a. 


U>  T  = 


(  f  )  (  f  ) 

The  local  extremum  values  of  p^/(t(()  written  as  pj '(t„) 

lXr  0j  IX  /v 


are  then 


( f ) /  \  2S  «  JC2mit>\  .  , 

pj ,  (t  J  =  -  ism  (  -r—  )  -  i  sin 

i£  £  ext  *(i3-Xs)  L  V  1+i  ' 


(t5)1 


2S(  i+X)  ,  f  Um.  2S 

■*TiS-1=7  Bin  k  ~  J  = 


2S  (  jf2mn 

n(  i-,C)  in  V  i+j*  >  * 


(5.41) 
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Prom  (5.  4l),  it  is  seen  that  the  maximum  absolute  value  of  "flipped" 

cross- correlation,  !p^,;(t^)|  ,  is  achieved  for  the  value  of  m, 

1  max 

for  a  given  i  and  £,  for  which  the  angle  — is  "closest"  to  n/2  or 

—  Aj 

ln/2.  For  (i+/)  >  10  the  maximum  "flipped"  cross-correlation  for 
2mn  . 

^  TX  =  i+i  is 


!4l)(^)i 


2S 


max 


Tt.  i-J 


(5.42) 


i  +  /  >  10 


Case  b. 


2mn 

^  Tx =  n 


U)t 


For  this  case  the  local  extremum  values  of  pj  .  ;(t.)  are 

±f  f 


(f)/  \l  2S  r,  .  £2mn  N  ^  i2mn  M 

*  *rr^Fj  L1  sin  v  -n  >  - J  sln  C — j\ 

ext 


2S(i-jfi) 

/■:C  .4  v  sin 

«(i2-£2 ) 


/"  £2mit 
\  jt-i 


2S 


J  ~  «(i+j6) 


sin 


Cfr) 


(5.  45) 


As  in  Case  (a)  the  maximum  absolute  "flipped"  cross-correlation  is 
achieved  for  ~^£  "closest"  to  n/2  or  3n/2.  For  i  +  £  >  10,  the  maxi- 


mum  absolute  value  for  o>  t  =  -7-7 

L  £-i 


is 


2S 

n( i+£) 


(5.44) 


max 
i  +  i  >  10 

Clearly  for  i  +  f,  >  10,  the  absolute  maximum  value  indicated  by  (5.42) 
is  greater  than  the  value  indicated  by  (5.44).  However,  it  can  be 


verified  that  for  all  i,  jC,  the  maximum  value  achieved  by  Case  (a)  is 
always  greater  than  that  achieved  by  (b). 

It  is  now  of  interest  to  compare  the  noiseless  channel  per¬ 
formance  of  k  Basic  Walsh  functions,  k  oo  with  a  similarly  chosen 
set  of  sinusoids.  The  frequencies  of  the  Basic  Walsh  functions, 

r*  l 

S  ,  S  ,  ,  S  ,  are  proportional  to  1,  2,  ...  ,  2r~  ,  re- 

ak-i  ak-2  VE 

spectively.  Thus,  a  similarly  chosen  set  of  sinusoids  would  be  Sa(t), 
S2(t),  ...  ,  S^p(t)  (F  =  0,  1,  ...  ,  k  -  1,  k  -  oo)  where  S  (t)  is 
given  by  (5.35). 

Without  any  external  noise,  the  output  of  the  i  correlator 
z^,  given  by  (5.8),  can  be  represented  as 

zi  =  z2F  =  ^F’S  +  I2F  >  (5.45) 


where  I^p  is  the  interference  terra  due  to  the  presence  of  the  other 
sinusoidal  carriers. 

The  maximum  value  of  the  interference  ll  p|  for  the 

53  max 


sinusoid  carriers  Is 


‘n? 


max 


k-1 

y 

F,=0 

F’/F 


(5.46) 


max 


Sr) 


For  F  >  4,  the  value  of  p  vi(t  )  specified  by  (5.42)  can  be 
”  3  '3  -  „„„ 


max 
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(5.  47) 


A  lower  bound  on  1 1  r»|  ,  denoted  as  which  is 

1  2  'max  s 


an 


accurate  representation  for  F  >  6  is  found  as 


i<P 

2 


23  r  ?  1  .  __1_  V  .-C-]  2S  F+l 

“  n  L  L  "f  0}  +l  L  J  "  «  0T 


(5.48) 


d=l 


2“  2 


c=o 


An  upper  bound  for  |I  p|  ,  denoted  as  1^,  is  found  as 

max  2 


2  ! 


00 


T(u)  _  2S  f  V  _1_  +  _1  \  2-c-l  _  2S  _TJ± 

AaF  "  n  L  L  2F-1  2F  L  J  "  *  2^-1  ‘ 


2r-x  2F 

d=l  c=o 


(5.49) 


* 


For  the  Basic  Walsh  function  carriers  the  maximum  interference 

S  i2 

has  been  determined  (as  a  modification  of  5-9)  as  ■=•  y,  when  we  are  cor- 

2E-1 

relating  with  the  carrier  associated  with  sequence  S  ,  E  =  1,  2,  . . .  , 

\-E 

k,  k  -*  00.  A  comparison  of  maximum  interference  of  the  two  types  of 
carriers,  for  several  small  values  of  E,  and  for  large  E,  where 
E  =  F  +  1,  is  presented  in  Table  5. 1. 
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TABUS  5.1 

Maximum  Interference  For  Basic  Walsh 
Function  and  Sinusoidal  Carriers 


1 

2 

3 

>  >  0 


|l  pi  -Walsh 
1  8* 'max 


S 

S 

0.75  •  S 
SE 


I  -sinusoids 
a**  'max 


0.92  *  S 

1.04  •  S 
0.73  ‘  S 

2SE 


2E‘: 


it  2 


E-l 


It  is  concluded  then  that  the  sinusoid  carriers  provide  slightly 
less  maximum  interference  for  large  E.  However,  the  interference  for 
the  sinusoid  carriers  can  introduce  an  error  in  the  absence  of  addi¬ 
tive  noise,  when  the  'eceived  asynchronous  signal  is  correlated  syn¬ 
chronously  with  Sa(t). 

In  order  to  determine  the  average  probability  of  error,  when 
correlating  with  sinusoid  S^(t)  in  the  presence  of  white  Gaussian  noise, 
it  is  necessary  to  determine  p(zildi=+l)  and  then  apply  (5.11).  However, 
the  determination  of  p(z^|di=+l),  for  a  finite  number  of  carriers,  re¬ 
quires  the  derivation  of  the  pdf  of  the  cross-  correlation  function  be¬ 
tween  two  sinusoids.  This  derivation  appears  difficult. 

However,  as  the  number  of  carriers  increase  without  limit,  then 

p(z^|d^=+l)  approaches  a  Gaussian  distribution  with  mean  S,  and  variance 

a2  as  derived  below. 
zi 
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The  variance,  c2/f\  of  the  "flipped"  cross-correlation  func- 


>Lf) 


tion  between  sinusoids  S^(t)  and  S^,(t),  is 


Vf) 

u 


(5.  50) 


( f ) 

Substituting  for  '(t^),  it  is  seen  that 


°p<*> 

U 


ksz 


t?  T^(  i2-Jl2) 


i3  sin2£  cu  t  +  .C2  sin3  &  03  x 
l  £ 


-  21  £.  si/l  i  03  • 

_  2S^  #  ( i2  +J? ) 
n3  (i2-*2)2 


sin  jC  a)  t. 

£ 


(5.  51) 


If  it  is  considered  that  the  probability  of  a  "flip"  in  S^,(t)  is  1/2, 
then  the  variance  of  the  correlation  function,  a2^,  between  S^t)  and 
S^(t)  (considering  both  "flipped"  and  "unflipped"  cross-correlation)  is 
then 


0*  .ij 

U  2  It*  (i 3-l3f 


Then  the  variance 


c2  is 
zi 


all  x/i 


(5.  52) 


(5.  53) 


*Note  that  r*  is  uniformly  distributed  over  the  interval  0-T  and  the  moan 
of  the  flipp'd  cross-correlation  function  is  zero. 
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If  we  set  i,  i  =  2  ,  F  =  0,  1,  ...  ,  k  -  1,  k  -  oo,  then  for  F  <  <  k  -  1, 


NS  q3 
o2  =  o2  *  +  S. 

Zi  a?  2T  n2 


r  ^  22^  +  2s  P-2  ^ 

L  Zi  24P  +  24*'~4C  -  24^”ac+1 

c=l 


oo 


^  j  o4 


gsF+ad  +  g3  F 


d=l 


24F+cd  +  24F  _  24F4«d+i 


]  .  (5.54) 


A  cougar ison  of  the  corresponding  variances  of  the  Basic  Walsh  Func¬ 
tion  Carriers  (see  5.32)  and  the  sinusoid  carriers  (see  5.54)  for 
E  =  F  +  1  =  1,2,  3  and  4,  is  presented  in  Table  5.  2. 

Table  5.2 

Values  of  Corresponding  Variances  for 
Basic  Walsh  Functions  and  Sinusoids 


E 

Variance- Walsh 

Variance- 

-Sinusoid 

1 

0.055  •  S2 

+  N0S/2T 

0.066  •  S2 

+  NqS/2T 

2 

0.055  *  S2 

+  NqS/2T 

0.073  •  S2 

+  NqS/2T 

3 

0.024  •  S2 

+  NqS/2T 

0.026  •  S2 

+  NqS/2T 

4 

0.0087-  S2 

+  NqS/2T 

0.0081-  Sa 

+  NqS/2T 

The  average  probability  of  error,  given  by  (5.34)  is  plotted 
for  sinusoids  S^t),  Sa(t),  S4(t)  and  Se(t)  in  Figure  (5.9).  Comparing 
these  curves  with  the  corresponding  curves  of  Figure  (5.3)  it  is  seen 
that  for  k  -*  oo  the  Basic  Walsh  functions  provide  a  lower  probability 


ERROR  PROBABILITY 
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F  /  G.  5.9. 

PROBABILITY  OF  ERROR  (PER  BIT)  /NOISE 
POWER  DENSITY 
CORRELATING  WITH  SINUSOIDS 
S  (t),  S  (t),  S  (t),  S  (t) 


Ill 


of  error  than  the  sinusoids  for  the  carrier  of  longest  period,  al¬ 
though  the  error  probabilities  are  approximately  equal  for  the  re¬ 
maining  cases. 


5. 4  Asynchronous  Binary  Signalling  V/ith  Random  3inary  Sequences 
In  this  section  the  average  performance  of  the  binary  asyn¬ 
chronous  signalling  system  will  be  determined,  when  the  carriers  are 
mapped  randomly  chosen  binary  sequences.  The  results  will  be  compared 
with  the  results  derived  in  Section  5.  2  for  the  Basic  Walsh  functions. 

The  cross-correlation  function  p^(X  )  between  the  mappings  of 
randomly  chosen  binary  sequences  and  of  length  L  (considering  the 
possibility  that  is  "flipped")  will  now  be  determined. 

Let  S<m)  and  S(“>  be  random  rhosen  mapped  binary  sequences. 


!1  *  [a(o1)'  -(i) 


a'  ' 
®1  * 


...  ,  a 


L-i  J 


=  [a(->,  a<‘>,  ...,aW], 


U)“| 


where  a^,  a^  =  ±1,  o  =  0,  1,  ...  ,  L  -  1,  with  equal  probability 


.(i) 


M 


and  a^  ’  and  a^,  are  statistically  independent,  for  o  =  o’,  i  /  / 
and  for  o  /  o',  all  i  and  l.  Then  the  cross-correlation  function  be¬ 
tween  and  sj^  at  a  delay  of  ^  is 


Xd-1 

p»w-V  )  41)a(t) 


i  r  v 


L-X.+b 


b=o 


+ 


L 


l-i 

) 

b=\e 


a(') 

L-  X^+b 


(5.  55) 
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(!)  (3) 

where  S  is  the  signal  power, d  ,  d  =  _l1  with  equal  probability  and 

(1)  (s)  » 
d^  and  d^  are  statistically  independent. 


is 


The  mean  value  of  p^(  \  )  is  zero.  Hence  the  variance,  o^ 

^  1  At 


11 


s3 

F 


s3 

L3 


V1 

VI 

y 

b=o 

> 

C--0 

L-l 

L-l 

.  \  ^ 

L 

/ 

b=X< 

ii 

o 

k(i)  a(i)  a(jl)  eiZ)  1 
c  L-X^+b  if-  X^+cJ 


fJO  a(i)  a(?)  ^ 

LFb  c  L-X.+b  ^-X,+cJ 


dj^  d^  S2 

£  £ 


\rl  L-l 


>  l  Ls J  ac  - 

b=o  c=X 


(i)  (l)  U)  U) 


L-\„+b  L-X. 

^  Jl 


[+c]  (5.56) 


Since  the  symbols  of  the  carriers  were  chosen  randomly,  each  of  the  ex¬ 
pectations  taken  over  the  ensemble  of  carrirers,  is  non-zero  only  if 

b  =  c.  Thus  oa  reduces  to 
Pi£ 


a  S2 

O3  =  73- 

pli  L 


V1  x£-i 

I  I  *•  + 

b=o  c=o 


S3 

be  IF 


L-l  L-l 

y  y  v 

b=X  c=X 
£  £ 


d(o  d(=)  s2  V1  L;x 

~  1?  2  }  6bc  ' 


(5.57) 


b=o  c=X 


'£ 


*It  is  recalled  that  if  d^1 )  and  d(2)  are  of  opposite  sign  then  has 
"flipped"  at  a  delay  of  X^  digits. 
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where  6.  is  the  Kronecker  delta  function.  Thus 
be 


=  rT 


[X.  +  L  -  X;.  +  o)  .  £_ 


(5.58) 


It  can  be  shown  that  X  )  is  in  general  binomially  dis¬ 
tributed,  but  for  large  L  the  binomial  distribution  approaches  Gaussian. 

th 

The  output,  Zj  of  the  i  correlator  for  L  -  oo,  is  then  Gaussian  dis¬ 
tributed  with  mean  ±  S,  and  variance  a2  , 

zi 


N  S  k‘1 


o  (k-l)Sa 
2T  L 


(5.59) 


The  average  probability  of  error,  P^)  is  then 


P^ei^  2  “  2  erf 


2(-fc-U  +  A 

L  G 


(5.60) 


Ihus,  if  the  ratio  of  the  number  of  sequences  on  the  channel 
(k)  to  sequence  length  (L)  approaches  zero  as  L  -*  00,  then  the  aver¬ 
age  probability  of  error  approaches  the  binary  PSK  error  probability. 
Although  the  Basic  Walsh  functions  provide  an  average  probability  of 

error  (averaged  over  all  k  carriers)  approaching  the  binary  PSK  re- 
/ST 

suit  (for  G  =  y  jj-  <  3),  only  log®  L  carriers,  L  -•  co,  are  used. 

Equation  ( 5. 60)  indicates  an  "average"  set  of  randomly  derived  mapped 
binary  sequences  will  provide  this  error  probability  although  permitting 
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many  more  than  logs  L  carriers  on  the  channel. 

37 

Some  questions  have  been  raised  by  Wolf  and  Elspas"  on  the 

validity  of  determining  the  performance  of  random  carriers  in  the 

manner  specified  in  this  section.  Ihe  analysis  appears  to  indicate 

2  N0S  (k-l)S2 

that  z.  is  Gaussian  distributed  with  known  variance  0  =  — =-* — 

1  Z  Cl  JL> 

However,  the  variance  of  the  interference  due  to  the  presence  of 

the  other  carriers  is  not  a  known  factor  but  a  random  variable  with 

mean,  0*  *  and  variance  •  For  a  large 

number  of  carriers  on  the  channel,  the  fluctuations  around  the  mean 

_ _  f 

of  o*  are  small  since  the  ratio  of  £(c^-o®)  J  to  o|  is  [2/(k-l)]^  • 

"th 

ftius,  we  can  approximate  o*  by  its  near  and  then  the  output  of  the  i 
correlator  is  approximately  Gaussian. 


5.  5  Asynchronous  Signalling  with  Linear  Error- Correcting  Codes 
In  this  section  the  noiseless  channel  performance  of  the 
asynchronous  signalling  system  with  mapped  cyclic  error- correcting 
codes  as  carriers  will  be  determined.  The  "unflipped"  cross-correla¬ 
tion  function  will  be  shown  to  be  related  to  the  minimum  distance  of 
the  code.  A  bound  on  the  "flipped"  cross -correlation  function  will  be 
related  to  the  solid  burst  error-correction  capability  of  the  code. 

The  theory  will  then  be  used  to  derive  the  performance  of  several 
Bose-Chaudhuri  codes,  the  largest  known  class  of  cyclic  codes,  for 
the  asynchronous  application.  A  brief  description  of  binary  cyclic 
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codes  is  first  given. 

A  linear  code  C  is  called  a  cyclic  code,  if  for  each  code 

word  Ci  =  (aQ,  ax,  ...  ,  )  in  C,  the  code  word  Ca  =  (aL_i, 

aQ,  ...  ,  aL_2 )  also  in  C. 

Every  cyclic  code  can  be  generated  by  a  raonic  polynomial, 

g(x),  in  the  algebra  of  polynomials  modulo  x^-1,  where  L  is  the  length 

of  each  code  word.  If  g(x)  of  degree  r,  divides  xL-l,  then  g(x)  gen- 

* 

erates  a  cyclic  code  of  dimension  k’  =  L  -  r.  Each  of  the  2  code 

words  is  expressible  in  the  form  g(  x)  [bQ  0biX  0 . . .  0  b^,  x  x  ] 

where  b^,  i  =0,  1,  ...,  k'-l,  is  either  0  or  1. 

24 

It  has  been  shown  that  each  code  word  can  be  generated  by 

a  k1  stage  shift  register  with  feedback  connections  corresponding  to 

the  polynomial  h(x)  =  (x  -l)/g(x),  with  the  proper  initial  loading. 

For  the  asynchronous  signalling  application  we  are  only  interested  in 

code  words  which  are  inequivalent  under  any  cyclic  permutation.  The 

number  of  "inequivalent"  code  words  (hereafter  designated  as  root  words) 

l8 

in  a  cyclic  code,  is  related  to  the  cycle  set  structure  of  h(x). 

A  cycle  of  the  polynomial,  h(x),  is  the  set  of  polynomials 
x^  hi(x),  mod  h(x)  where  d  =  0,  1,  ...,  and  hi(x)  is  any  polynomial  of 
degree  less  than  the  degree  of  h(x).  The  length  of  a  cycle  is  the 
smallest  integer  v,  (v  always  divides  2  -l),  such  that  x  hi(x)  H 

hi(x),  mod  h(x).  [The  root  code  word,  considered  as  a  periodic  se¬ 
quence,  generated  by  the  k'  shift  register  with  feedback  connections 
prescribed  by  h(x)  and  initial  loading  corresponding  to  xd  hx(x),  will 
have  a  minimum  period  or  length  v.  ] 


As  an  example,  consider  the  cycles  of  the  polynomial 
h(x)  “  x4  +  x3  +  x2  +  x  +  1,  over  GF(2),  shown  in  Table  5.3. 


Table  5.3 

Cycles  of  h(x)  =  x4  +  x3  +  x2  +  ::  +  1 


hx(x) 

0 

1 

x  +  1 

x  hx(x) 

0 

X 

a 

X  +  X 

x3hx(x) 

a 

X 

3,3 

X  +  X 

x3hx(x) 

3 

X 

X2  +  X  +  1 

x4hx(x) 

X3  +  X8  +  X  +  1 

3  2 

X  +  X  +  x 

x5hx(x) 

1 

X  +  1 

x3  +  1 
x3  +  xa  +  1 
x2  +  1 

3  , 

X  +  X 
X3  +  X  +  1 

x3  +  1 


From  Table  5.  3  it  is  seen  that  x4  +  x3  +  x3  +  x  +  1  has  4 
cycles  (l  of  length  1  and  3  of  length  5).  Ihe  cycle  set  structure  is 
denoted  as  l(l)  +  3(5).  From  thi3  cycle  set  analysis  it  can  be  deter¬ 
mined  that  the  cyclic  code  of  length  15  with  generating  polynomial 
g(x)  =  (x15-l)/(x4  +  x3  +  x3  +  x  +  1)  =(x.+  l)  (x2  +  x  +  1)  (x4  +  x  +  1) 
(x4  +  x3  +  l)  has  4  root  words,  one  of  which  is  the  all  zero  code  word. 

In  general,  an  irreducible  polynomial  of  degree  k'  will  have 
1  cycle  of  length  1  and  u1  of  length  v',  denoted  as  l(l)  +  u'(v‘)> 
where  u'  and  v'  ore  in  general  difficult  to  determine.  However,  if 
h(x)  is  a  primitive  irreducible  polynomial  of  degree  k’,  then  u'  =  1 

k»  ~[fl 

and  v'  =  2  -1.  Elspas  has  shown  that  a  polynomial  h(x)  vhich  is 

the  product  of  two  different  irreducible  polynomials,  hi(x)  and  h3(x), 
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with  cycle  set  structures  l(  1)  +  ui(vi)  and  l(l)  +  ^(va)  respectively, 
will  have  a  cycle  set  structure  of  l(l)  +  ui(vi)  +  uc(va)  +  m(v),  where 

M  =  Hi  Me  [G.C.  D.  (vx>  vs )  ] 
v  =  L.  C.  M.  ( Vj ,  Vg ) 

This  technique  is  easily  extended  to  derive  the  cycle  set 
structure  of  a  polynomial  which  is  the  product  of  an  arbitrary  number 
of  different  irreducible  polynomials. 

If  the  minimum  distance  of  a  cyclic  code  is  known  then  the  mini¬ 
mum  distance  betveen  cyclic  shifts  of  root  words  iB  known.  However,  for 
binary  signalling  a  root  word  or  it6  complement  is  transmitted,  and  only 
if  the  complement  of  each  code  wor'1  is  in  the  code  will  the  minimum  dis¬ 
tance  between  cyclic  shifts  of  any  root  word  and  the  complement  of  another 
root  word  be  know.  Hie  following  theorem  gives  the  condition  for  a  cyc¬ 
lic  code  to  contain  the  complement  of  every  code  word. 

THEOREM  5.2  Hie  complement  of  every  code  word  of  a  binary  cyclic  code, 
i6  a  code  word  if  the  generating  polynomial,  g(x)  does  not  con¬ 
tain  1  +  x  as  a  factor. 

Proof:  Since  the  code  is  a  linear  code  the  theorem  can  be  established 

by  proving  that  the  all  zero  and  all  one  code  words  are  in  the  code. 

If  1  +  x  is  not  a  factor  of  g(x),  then  1  +  x  is  a  factor  of  h(x)  = 
(x^-l)/g(x).  Hie  polynomial  1  +  x  has  two  cycles  of  length  1.  Hence, 
by  the  product  rule  of  cycle  sets,  h(x)  contains  two  cycles  of  length  1. 
However,  the  lengths  of  distinct  cycles  are  the  lengths  of  distinct  root 
words,  and  the  only  two  root  words  of  length  1  are  the  all  zero  and  all 


one  code  words. 


lift 


It  will  be  shown  in  the  succeeding  analysis  that  the  magnitude 

of  the  "flipped"  cross-correlation  between  mapped  root  words  can  be  re- 

38 

duced  by  using  repeated  mapped  code  words  as  carriers.  Hie  following 
theorem  relates  to  the  structure  of  a  repeated  cyclic  code. 


THEOREM  5. 3  The  code  formed  by  repeating  m  times  each  code  vord  of 
a  cyclic  code  of  length  L  with  generating  polynomial  g(x)  is 
also  a  cyclic  code,  but  of  length  m  L  with  generating  polynomial 

m-1  .  T 

g'(x)  =  ^  x  g(x) 
b=o 

The  number  of  root  words  in  the  repeated  code  is  identical  to 
the  number  in  the  original  code.  If  the  original  code  had 
minimum  distance  d,  then  the  resultant  code  has  minimum  dis¬ 
tance  md. 


It  will  now  be  shown  that  the  "flipped"  cross-correlation 
function  is  related  to  the  solid  burst  error- correct  ion  capability  of 
the  cyclic  code.  Consider  the  "flipped"  cross-correlation  funct' 

Pijfc  between  the  mappings  of  code  words  (in  polynomial  represtnca 

tion)  C^x)  and  C^(x)  at  a  shift  of  X  digits.  Hien  is  deter¬ 

mined  as 


pW(V  =  ■jkno.  of  -s  in  A(x)  -  no.  of  l's  in  A(x)]  (5. 6l 

where 

A(x)  e  Ci(x)  0  x  ^  C£(x)  0x^(1  +  x+  ...  +  xL-1  mod  (xL-l) 

(5.  62) 


Since  C^(x)  610(1  C^(x)  are  the  polynomial  representations  of  code  words 
from  a  cyclic  code,  A(x)  can  be  reduced  to 


A(x)  =  Cq(x)  0x  ^ [1 


+  x^"1  &],  mod  (xL-l) 


+ 


»  •  i 


(5.63) 
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where  C^(x)  corresponds  to  a  code  word  of  minimum  weight  d  and  maximum 

weight  L  -  d,  provided  1  +  x  is  not  a  factor  of  g(x). 

It  is  thus  seen  that  A(x)  is  formed  by  adding  the  polynomial 

* 

representation  of  a  solid  error  burst  of  length  L  -  X,.  A  theorem  on 
the  correction  of  solid  bursts  follows. 


THEOREM  5.4  A  double  error  correcting  cyclic  code  ( d  >  5)  can  correct 
all  solid  bursts  that  do  not  "wrap"  around  the  end  of  the  code, 
provided  the  generating  polynomial  g(x)  does  not  contain  1  +  x 
as  a  factor. 


The  above  theorem  is  generalized  to  multiple  solid  burst  error 
correction  by  the  following  theorem. 


THEOREM  5.5  At  error  correcting  cyclic  code  (t  even)  can  correct  all 
error  patterns  that  consist  of  t  or  fewer  solid  bursts,  each 
starting  at  the  beginning  of  the  code  word,  provided  g(x)  does 
not  contain  1  +  x  as  a  factor. 


Proof:  Since  the  code  can  correct  t  or  fewer  errors,  then  the  error  pat¬ 

tern  E(x), 

E(x)  =  Z(+/X  +  X  f  ^  9  (5.64) 


where  t'  <  t 

Xf  <  L  -  1  , 

is  correctable  by  reducing  the  received  vector  E(x)  +  C  (x)  (in  poly- 

m 

nomial  form),  modulo  g(x).  However,  E(x)  can  be  represented  as 


*A  solid  error  burst^  is  defined  as  an  error  burst  where  every  digit 
in  the  burst  is  in  error. 
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t* 

E(x)  =  V^(l  +  x)(!  +  x  + 


♦  xS 


(5.65) 


Since  1  +  x  and  g(x)  are  relatively  prime,  the  error  pattern  E'(x), 


E'(x) 


E(x)  = 
1  +x 


t' 


+  x  +  . . . 


(5. 66) 


is  correctable,  indicating  that  t  or  fewer  solid  bursts  starting  at 
the  beginning  of  the  code  word  are  correctable. 


Prom  the  preceding  theorem,  upper  and  lower  bounds  on  the  weight 
of  a  vector  formed  by  adding  (mod  2)  a  solid  burst,  starting  at  the  be¬ 
ginning  of  the  code  word,*  to  a  code  word  (except  the  all  zero  or  all 
one  code  words)  is  now  derived. 


THEOREM  5.6  Consider  the  vector  which  is  formed  by  the  mod  2  addition 
of  a  code  word  [from  a  t  error-correcting(t  even)  cyclic  code 
for  which  1  +  x  is  not  a  factor  of  g(x)]  and  a  solid  burst 
starting  at  the  beginning  of  the  code  word.  ^The  minimum  and 
maximum  weights  of  the  resultant  vector  are  -jr  +  1  and  L  -  £  -  1 
respectively. 


Proof:  Jbr  t  even,  it  has  been  6hown  that  t  or  fewer  "beginning"  solid 

bursts  are  correctable.  Hence,  the  vector.  A,  corresponding  to  a  code 
word  (which  is  not  che  all  zero  or  all  one  code  words)plus  t  "beginning" 
solid  bursts  must  be  different  from  the  vector  corresponding  to  the  all 
zero  code  word  plus  t  or  less  "beginning"  solid  bursts.  However,  any 

♦The  bound  when  a  solid  burst  starts'  in  the  middle  of  a  code  word  and 
proceeds  to  the  end  will  be  identical  to  the  bound  which  is  derived 
for  the  "beginning"  solid  burst. 


I 
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vector  of  weight  t/2  can  be  formed  by  the  proper  mod  2  summation  of  t 
"beginning"  solid  bursts.  Hence,  for  A  to  be  different  from  the  all 
zero  code  word  plus  any  t  or  few  "beginning"  solid  bursts,  A  must  have 
weight  at  least  ~  +  1.  If  the  conqaarison  is  made  with  the  all  one  code 
word,  then  the  maximum  weight  of  A  is  derived  as  L  -  |  -  1. 

From  the  preceding  analysis  it  is  now  possible  to  specify  the 
bounds  on  the  "unflipped"  and  "flipped"  cross-correlations  for  a  cyclic 
code  of  length  L,  for  which  1  +  x  is  not  a  factor  of  g(x),  with  a  minimum 
distance  d  =  2t  +  1,  t  even. 

The  "unflipped"  cross-correlation  function  is  bounded  by  ± 

From  Theorem  5. 6  it  can  be  shown  that  the  "flipped"  cross -correlation  is 

bounded  by  ^  ^  “ . 2  . 

L 

For  the  repeated  cyclic  code.  Theorem  5.3,  the  "unflipped"  cross- 

L  •  2d 

correlation  is  still  bounded  by  i:  — ^ - .  However,  the  bounds  on  the 

"flipped"  cross-correlation  are  lowered  when  the  cyclic  code  is  repeated 
m  times.  The  maximum  absolute  value  of  "flipped"  cross-correlation  occurs 
when  the  relative  delay  of  the  code  words,  is  either  0  <  \^  <  L,  or 
(m  -  1)L  <  ^  <  m  L.  It  is  then  determined  that  the  "flipped"  cross-cor¬ 
relation  is  bounded  by 


For  large  values  of  m,  the  "flipped"  cross-correlation  approaches  the 
"unflipped"  value. 

From  the  bounds  on  the  "flipped"  cross-correlation,  the  maxi¬ 
mum  number  of  carriers  which  can  be  asynchronously  multiplexed,  on  a 
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noiseless  channel  with  zero  probability  of  binary  error  can  be  derived. 

It  should  be  noted  if  t  is  odd  then  t  is  replaced  by  t  -  1  in  all  of 
the  cross-correlation  bounds. 

As  eui  illustrative  example  consider  the  (31, 11)  Bose-Chaudhurl 

co<Li  for  which  t  =  5  fluid  d  =  11.  The  generating  polynomial  for  this 
40 

code  is 

g(  x)  =  ( 1+x3  +xB )  ( l+xa  +x3+x4  +xB )  ( l+x+xa  +x4  +xB )  ( 1+x+x2  +x3+xB )  (  5 

The  shift  register  polynomial  is  then 

h(x)  °  *  (l+xHl+x+v3+x4+xB)(l+x3+xB)  (5.68) 

Hie  number  of  root  words  (not  including  complements)  are  the  number  of 
cycles  of  the  polynomial  h’(x), 

h 1  ( x)  =  ( 1+x+x3  +x4  +xB )  ( 1+x3  +x6 ) 

Since  each  of  the  factors  of  h’(x)  are  primitive  polynomials,  it  is 
immediately  apparent  that  the  cycle  set  structure  of  h'(x)  is 

1(1)  +  1(31)  +  1(31)  +  31(31) 

Hence,  there  are  34  root  words  in  the  code. 

Under  the  no  "flipping"  condition,  the  bound  on  the  maximum 
number  of  ccuriers  is  4.  Under  the  "flipping"  condition,  the  maximum 
number  which  ceui  appear  is  2.  If  the  repeated  code  is  used  with  m  =  3 
(code  length  =  93)  then  3  code  words  csui  appear.  Hie  code  length  must 
be  incree^ed  to  (ll)(3l)  =  34l  to  enable  4  code  words  to  appear  with 
zero  probability  of  error.  In  Table  5.  4  the  bounded  asynchronous 
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performance  of  several  other  Bose-Chaudhuri  codes, which  give  the 
"optimum"  performance  for  the  given  length,  are  tabulated. 

Table  5.4 

Asynchronous  Performance  of  Root  Words 
of  Bose-Chaudhuri  Codes 


Code 

t 

No.  of 
Root 
Words 

Max.  No. 
for  "un¬ 
flipped" 
Case 

Max.  No. 
for 

"flipped" 

Case 

Length  for  3 
Words 

Length  for  4 
Words 

(15,7) 

2 

6 

4 

2 

3*15  **  45 

00 

(31,11) 

5 

34 

4 

2 

3*  31  =  93 

11-31  =  341 

(63,10) 

13 

ID 

8 

2 

2-63  =  126 

4-63  =  252 

The  performance  characteristics  listed  in  Table  5.4  indicate 
that  the  noiseless  channel  performance  of  the  mapped  Bose-Chaudhuri  codes 
is  inferior  to  the  performance  of  the  Basic  Walsh  functions,  under  the 
"flipping"  condition,  in  that  fewer  carriers  can  appear  on  the  channel 
with  the  assurance  that  the  probability  of  error  is  zero.  However,  an 
advantage  of  vsing  Bose-Chaudiurl  codes  is  that  root  words  can  be  chosen 
which  have  least  period  equal  to  the  length  of  the  code.  (The  Basic 
Walsh  functions  do  not  exhibit  this  property. )  Since  the  autocorrelation 
function  of  the  mapped  root  words  is  low  for  non-zero  delays  (it  is 
bounded  by  ±  — ) ,  synchronization  of  a  particular  transmitter- 

receiver  pair  can  be  gained. 

It  should  be  noted  that  the  bound  on  the  "flipped"  cross-cor¬ 
relation  is  a  pessimistic  bound.  By  proper  choice  of  the  root  words. 
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the  absolute  maximum  value  of  "flipped"  cross- correlation  will  be  sig¬ 
nificantly  reduced.  It  was  found  that  by  choosing  a  preferred  set  of 
root  wrods  of  the  (15,7)  code,  it  was  possible  to  place  3  code  words 
on  the  channel  when  the  word  length  was  30,  instead  of  the  value  of 
45  indicated  in  Table  5.4. 


6.  HARD  LIMITING  MULTIPLEXING  OF  N-ARY  SEQUENCES 


6. 1  Introduction 

12  23 

Titsworth  *  has  described  a  binary  synchronous  multiplex 
system  which  is  applicable  for 

a.  the  transmission  of  several  signals  from  several 
different  transmitters  over  the  same  channel. 

b.  the  simultaneous  transmission  of  messages  to  several 
receivers  in  a  single  sequence  stream. 

There  are  k  users  on  the  channel,  each  using  a  different 

lc 

binary  mapped  Basic  Walsh  function  of  length  2  as  a  carrier.  Infor¬ 
mation  is  sent  by  transmitting  either  ihe  mapping  of  the  sequence  it¬ 
self,  or  the  mapping  of  the  complement  sequence.  (Hence,  1  bit  of 
information  is  sent  by  each  carrier  in  2  time  slots. )  The  multiplex¬ 
ing  operation  is  performed  by  operating  on  the  sum  of  the  k  mapped  se¬ 
quences  (with  the  provision  that  some  may  be  complemented)  with  a  Boolean 
logic  device.  The  output  of  the  logic  device  is 

a.  plus  one  if  there  is  a  majority  of  plus  ones  in  the  k 
mapped  sequences  during  a  particular  time  slot, 

or 

b.  minus  one  if  there  is  a  majority  of  minus  ones  in  the 
k  mapped  sequences  during  a  particular  time  slot. 

For  an  equal  number  of  plus  ones  and  minus  ones, the  output  is  arbitrary. 

For  k  odd,  the  logic  operation  can  be  performed  by  a  hard  limiter.  (The 

use  of  a  hard  limiter  is  particularly  desirable  for  multiple  access  to 
33 

a  satellite  because  of  the  nonlinear  characteristic  of  the  travelling 
wave  tube  usually  used  for  economical  transmission. ) 


The  single  sequence  binary  stream  is  then  transmitted  to  k  re¬ 
ceivers,  which  are  k  filters  (or  correlators)  matched  to  the  k  mapped 
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basis  sequence  used  as  initial  carriers.  Hie  output  of  a  particular 
correlator  is  a  positive  or  negative  value,  depending  upon  vhether  a 
mapped  basis  sequence  or  its  complement  was  initially  transmitted,  in¬ 
dependent  of  the  information  state  of  the  other  carriers. 

It  will  be  shown  in  this  chapter  that  a  hard  limiting  syn¬ 
chronous  multiplex  system  can  be  designed  with  N-ary  sequences  with 
symbols  mapped  onto  the  roots  of  unity  or  phase  modulated  cosine  func¬ 
tions.  Hie  k  (mapped)  cyclically  N-orthogonal  basis  sequences,  S  , 

cL © 

S  ,  . . .  ,  S  ,  of  length  are  the  carriers*  to  be  multiplexed.  A 
*1  ak-i 

mapped  sequence  or  one  of  its  N-l  complements  is  transmitted.  (Hence 
logs  N  bits  of  information  are  transmitted  with  each  carrier  in  N^  time 
slots. )  As  in  the  Boolean  system,  the  optimum  multiplexing  is  achieved 
by  a  logic  device  operating  on  the  sum  of  the  components  of  the  carriers. 

Hie  output  of  the  logic  device,  f(  )  during  time  slot  g, 
g  =  0,  ...  ,  -1,  is 

a.  a  unit  amplitude  vector  with  a  phase  angle  equal  to  the 
angle  of  the  vector  which  is  the  linear  sum  of  the  com¬ 
ponents  of  k  carriers  (arbitrarily  complemented)  during 
the  particular  time  slot,  for  the  roots  of  unit  mapping, 

or 

b.  a  unit  amplitude  cosine  function  phase  modulated  by  the 
angle  defined  in  (a)  above. 

It  can  be  seen  that  the  logic  device  for  the  roots  of  unit  mapping  acts 
as  a  "hard- limiter",  provided  G.  C.  D.  [N,k]  =  1.**  Hie  interpretation  of 
the  optimum  logic  for  the  phase  modulated  cosine  mapping  will  be  discussed 


*  Throughout  the  remainder  of  the  chapter,  the  word  carrier  will  be  used 
to  indicate  a  mapped  sequence  and  complement  carrier  will  be  used  to 
indicate  a  mapped  complement  sequence. 

*#  This  provision  ensures  that  the  linear  sum  of  the  components  of  the 
carriers  is  not  zero.  The  output  of  the  logic  device  is  arbitrary 
if  the  sum  of  the  components  is  zero. 
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In  the  succeeding  sections. 

The  receiver  which  is  to  decide  on  the  information  contained 
in  the  ith  carrier  is 

a.  A  cross-correlator  between  f(Sj.)  and  mapped  onto  the 
roots  of  unity  The  phase  of  the  correlator  output 
determines  the  decision  on  which  complement  of  8^  was 
sent, 

or 

b.  A  set  of  N  cross-correlators  between  f ( S^)  and  8^  and 
its  N-l  complements  all  mapped  onto  phase  modulated 
cosine  functions.  The  correlator  which  yields  the 
largest  output  determines  the  decision  on  the  information 
contained  in  the  particular  carrier 

The  complete  multiplexed  system  is  shown  in  Figure  6  1  for 
roots  of  unity  mapped  sequences.  '  In  the  following  sections,  the  results 
summarized  in  thi6  introduction  are  derived  and  the  system  performance 
is  analyzed.  For  the  development  of  the  theory  it  will  be  assumed  that 
the  N-ary  sequences  are  mapped  onto  the  roots  of  unity  When  allusion 
is  made  to  real  sequences  the  theory  will  be  appropriately  qualified 
to  the  phase  modulated  cosine  mapping. 


6. 2  Transform  Analysis  of  N-ary  Sequences 

A  transform  pair  will  be  derived  which  will  be  used  to  evalu¬ 
ate  the  performance  of  the  multiplexed  N-ary  system.  The  transform 
pair  will  be  shown  to  relate  to  the  optimum  logic  function  and  the  mag¬ 
nitude  of  the  correlator  output. 

Define  a  function,  d>( s, x)  of  the  vectors 

s  =  (s^Sj,  ...  ,  s^j )  ,  s^  =  0, 1,  . . .  ,  N-l  , 
x  =  (x^,xj,  ...  ,  Xj^_j)  ,  xi  =  0,1,  . . .  ,  N-l  , 


ARY  MULTIPLEXING  SYSTEM 


ORRELATi 
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<Ks,x)  =  N"k/3  exp[j  y’  (fl0x0  +  +  8k-i  xk-i^  f6,1) 

It  is  apparent  that  there  are  functions  9(sc,x)  of 

x  for  a  given  vector  s  .  It  is  easily  verified  that  q>(s  ,x)  and 

c  c 

9(s.,x)  are  orthonormal,  that  is: 


l  V<?0,x>  9*(?d,x)  =  { 

all  1 

N*  x 


c  j-  d 

c  =  d 


(6.2) 


where  c,d  =  0,1,  ...  ,  N^-l,  and  V*(s*d,x)  is  the  complex  conjugate  of 
Similarly 


£  9(s,xe)  9*(s>xf)  =  | 

all  1 

I**  s 


e  /  f 


e  a  f 


(6.3) 


where  e,f  =  0,1,  . . .  ,  h  -1  . 

An  interesting  transform  pair  can  be  derived  with  function 
cf>(s,x)  as  the  kernel.  Define  a  function  F(s)  by: 

Ft?)  =  rk/3  \  f(x)  <p(?,?)  (6.4) 

La 

all. 

F*  x 


where  f(x)  is  a  unit  amplitude  vector  which  is  a  function  of  the  k 
components  of  x. 

Ihe  inverse  transform  f(x)  is  derived  as  follows.  Consider 
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l 

all 
N*  s 


<F*(s>xf) 


N"k/S 


1 

all 


)  f(xe)  v(s,xe)  ^s,xf) 
atl 


,  rk/J f(xe) 
all 


£  9(s,xe)  9*(s,xf) 
all 


Applying  the  orthonormal  property  (6.3)  it  is  seen  that 


)  F(s)  9*(s,xf)  = 


all 
^  8 


or 

f(x)  =  !^/3  ^  F(s)  q>*(8/x)  .  (6.5) 

all 
s 

The  equivalent  to  Parseval's  formula  is  derived  as 

^F(s)  F*( s )  =  y  f(x)  f*(x)  =  1  (6.6) 

all  all 

s  X 

6.3  Derivation  of  the  Optimum  Logic  and  Correlator  Output 
Consider  the  multiplexing  of  k  N-ary  sequences,  Sq, 

Sj,  ...  ,  S^-i  (at  present  arbitrary)  of  length  mapped  onto  the 

complex  roots  of  unity.  Then  if  the  N-ary  sequence  Sj  has  components, + 


+Note  the  change  in  notation  from  the  definition  of  a  sequence  given 
in  Chapter  2. 
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\  •••  then  the  component  of  the  mapping  of 

during  time  slot  g,  g  *  0,1,  ...  ,  N^-l,  is  =  exp  [J2n 
The  output  of  the  logic  network  during  time  slot  g  will  be  denoted  as 


'Tie) 


f(^vgy)*  It  is  clear  that  f(§^)  is  strictly  a  function  of 
Gi8)’  ••• .  4?D.  the  symbols  of  each  of  the  k  sequences  during  time 


o 

slot  g. 


A  set  of  sequences,  S^,  will  now  be  specified  and  the  result 


— •/  \ 

of  cross-correlating  f(vg  )  (for  all  g)  with  the  mapping  of  a  particu¬ 
lar  sequence  will  be  determined.  Choose  the  k  sequences  to  be  the  set  of 
k  cyclically  N-orthogonal  basis  sequences  of  the  N-orthogonal  sequences 

of  length  N*\  Thus,  in  the  notation  of  Chapter  4,  we  have  S  *  S  ,  , 

o  a 

o 

Sfc.i  =  S  ,  with  a  typical  sequence 
k-  x 

S.  *  S  =  (0,0,  ...  ,  0,  N— 1,  ...  ,  N-l,  ...  ,  1,  ...  1, )  , 

1  ai  — 


N1  N3 

repeated  times. 


r(g) 


At  the  receiver,  let  p_  be  the  cross-correlation  of 

fa. 


f(SVtt')  and  the  mopping  of  S  .  Thus,  p  becomes 

ai  fai 


(6.7) 


fa. 


~  y  f(C(e))  exp  C-j2n4e)/N] 
J'  g=o 


Noting  that  f(fc>^)  is  given  by 


f(46))  =  N^3  y  F(s)  9*(s,fc,[g)) 


ell 


(6.8) 


s 
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Y  exP  l-&*0  4g)  +  si  48^  +  •••  +  Sk-i  4?l)]  >  (6*9) 


all 

8 


p_  is  found  to  be 
fa. 


a*-: 


pfa1  '  J  1  £4'#Bo  4g>+  •••  +  (si+1>  4g>  + 


all  g=o 
— < 

8 


+  8k-i  4-i )  J 

(6.10) 


Now  define  a  vector,  v.,  in  the  s  space  as 

i 

v_L  e  (o,  ...  ,  N-l,  ...  ,  0) 

where  the  only  non- zero  term  occurs  in  the  ith  place.  Hence,  the  ex¬ 
pression  for  pfa  can  be  reduced  to 

N^-l 

pfa  'F(71)+^  l  F(?)  I  exph  + 


all  g=o 


s^i 


. .  +  (sj+1)  +  •••  +  sk-i  ^k?i)j 


(g)' 


(6. 11) 


where 


FCv^)  =  II"- /s  ^  f(x)  9(vt,x) 
all 
N^3? 


However,  if  the  sequences  are  chosen  as  the  k  cyclically  N-orthogonal 


133 


basic  sequences  with  components  given  by  (6.7)  then 
1^-1 

^  exp[-j  ■|1(so  ^  +  ...  +  (s±+l)  ^  +  ...  +  4?i) J  =  0 


since  for  all  s  /  v^,  the  symbols 

sq  ^g^  +  ...  +  ( s^+1)  ^g^  +  ...  +  sk-1  mod  N, 

for  g  *  0,  ...  ,  N*-!  are  the  symbols  of  one  of  the  N- orthogonal 

sequences. 

Hius,  the  value  of  the  correlator  output  for  the  k  cyclically 
N-orthogonal  sequences  is 


pfat  “ 


(6. 12) 


Hie  value  of  the  correlator  output  will  now  be  derived  if  the 

mapping  of  the  r  (r=l,  ...  ,  N-l)  complement  of  S  is  multiplexed 

ai 

instead  of  Sa  .  Hien,  ^g^  is  replaced  by  sj^-r  in  Equation  (6.7). 

th 

Also  under  the  condition,  the  r  complement  is  transmitted  the  expres¬ 
sion  for  f(t/g^),  Equation  (6.9)  is  changed  to  [f(^g^)l  , 


WS(s))J, 


=  £  F(?)  expj_-J  |ks0 


(g)  +  ...  ♦  s(  i,<g)  + 

0  1111 


all 


+  Sk_!  4?i)j 


(6.13) 


Hien  repeating  the  operations  defined  by  Equations  (6.10)  and  (6.11) 
with  the  substitution  [f(^g^)^,in  place  of  f(^8^),  the  output  of  the 
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1  correlator,  [p.  ]  ,  under  the  condition  that  the  r  complement  was 

rai  r 

transmitted  is 

tpfa^r  "  pfa  exPt"«3  2n  r/Nl  *  (6.14) 

It  is  noted  then  that  the  information  carried  by  the  i  carrier 

is  determined  uniquely  by  the  phase  of  the  correlator  output,  independent 

of  the  information  content  of  the  other  k-1  carriers.  The  contribution 
th 

of  the  i  sequence  to  the  total  output  is  then  maximized  by  maximizing 
the  value  of  over  all  possible  logic  functions,  f(x),  under  the 

condition  that  the  k  cyclically  N-orthogonal  sequences  are  used.  However, 
it  can  be  shown  that  for  any  set  of  input  sequences,  F(vi)  should  be  maxi- 
mized  to  maximize  the  contribution  of  the  i  input  independent  of  the 
condition  of  the  other  inputs.* 

Now  it  is  necessary  to  determine  the  logic  function  f(x)  which 
maximizes  F(v^),  where  f(x)  is  specified  as  a  unit  amplitude  complex 
vector.  If  we  assume  that  the  magnitude  o;f  the  outputs  from  all  cor¬ 
relators  is  to  be  identical,  then  F(v^)  =  F^,),  for  all  i,  i'.  We 
can  then  set 

k-1  _k  k-1 

HVj)  =  £  )  ^  )  f(x)  exp [-;}■—  x.,  ]  / 

i=o  i=o  aU 

x 

or  -k  _  kT1 

F( Vf)  =  ^  f(x)  ^  exp(-J  Y  xi^  (6.15) 

all  i=o 

x 

*The  proof  is  similar  to  the  proof  for  the  Boolean  multiplexed  system*^ 
and  hence  is  omitted. 
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Now  the  x^'s,  i  =  o,  ...  ,  k-1  are  the  symbols  of  each  of  the  k  multi- 

k-1  g  n 

plexed  sequences  during  a  particular  time  slot.  Hence  £  exp[-J  -rr  x  ] 

i=o  !t  A 

is  a  complex  vector  which  is  the  sum  of  the  mapped  symbols  of  the  k  se¬ 
quences  during  the  time  slot.  Thus,  F(Vj,)  is  maximized  if  for  all  ^ 

_  _  k-1  2n 

values  of  x  the  complex  vectors  f(x)  £  exp[-J  -rr  x  )  all  have  the  same 

i=o  w  1 

phase  angle.  If  the  phase  angle  of  F(v^)  is  to  be  zero  (this  is  an  ar¬ 
bitrary,  but  practical  choice),  then  F(v^)  is  maximized  if  f(x)  is  a 
vector  (unit  amplitude)  with  a  phase  angle  equal  to  the  negative  of  the 


k-1 


2n 


phase  angle  of  £  exp(-J  —  x.  ]•  Ihus,  the  optimum  logic  can  be  inter¬ 
im  N  l 

preted  as  a  hard  limiter  on  the  sum  of  the  complex  components  of  the  in¬ 
put  mapped  sequences,  so  that  the  output  vectors  of  the  logic  network 


k-1 


2n 


during  each  time  slot  are  all  of  unit  amplitude.  If  L  exp(j  x. )  =  0, 

i=o  w  1 

which  can  only  result  if  0. C.D.  [N,k]  j  1,  then  f(x)  can  be  chosen  to  be 
any  unit  amplitude  vector  since  there  is  no  contribution  to  the  maxi- 
mality  of  F(vi). 

It  is  of  interest  to  note  the  nature  of  the  optimum  logic  net¬ 
work  under  the  condition  that  phase  modulated  cosine  functions  are  used 
as  a  mapping  instead  of  the  Nth  roots  of  unity.  From  the  discussion  in 
Chapter  2  on  the  relationship  between  the  correlation  function  with  co¬ 
sine  function  and  roots  of  unity  signalling,  it  can  be  seen  that  the 
output  of  the  logic  network  f f ( x)  3C0S  vhich  maximizes  the  correlator 
outputs  for  cosine  function  signalling  is 


Cf(Z»cos  =  006  ^  +  * 


(6.16) 


k-1 


where  6 _  is  the  phase  angle  of  £  exp(,)  ~  x4).  However,  0_  can  be 

ii  i=0  W  1  la 


written  as 


k-1 

1 sln  (T  xi> 

®L  "  tan'1  Erf -  •  (6. 17) 

I cos  (T  xi> 

i=0 


Substituting  (6.17)  in  ( 6. 16) ,  and  expanding  cos(ujt  +  0^),  [ f*( x)  )cos 
becomes 


k-1 

yZn 

sin  —  x ^  sin  cut 

i=o 

k-1 

Y  cos  (cut  +  y  x±) 
i=o 

i=o  i=o  (6.18) 

Thus  the  optimum  output  of  the  logic  network,  under  the  cosine  func¬ 
tion  mapping,  is  proportional  to  the  linear  sum  of  the  mapped  components 
of  the  sequences  during  a  particular  time  slot.  Then  the  optimum  logic 
for  this  condition  can  also  be  interpreted  as  a  hard  limiter.  The  in¬ 
strumentation  of  the  logic  can  either  be  an  "infinitely"  fast  acting  AGC, 
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or  a  bandpass  limiter. 

For  either  mapping  it  has  been  shown  that  the  value  of  F(v^) 
(which  is  the  correlator  output  when  the  cyclically  N-orthogonal  se¬ 
quences  are  used),  specifies  the  system  performance.  It  is  seen  from 
(6.15)  that  F(vi)  is  a  function  of  N  and  k.  Values  of  F(v^)  for 
N  =  2,3,4  and  several  values  of  k  are  given  in  Table  6.1.  The  values 
for  N  =  2  are  taken  from  Titsworth's  papers. 


k 

OJ 

It 

S5 

Table  6. 1 

Values  of  F(vi) 

N  =  3 

N  =  4 

1 

1 

1 

1 

2 

0.5 

0.667 

0. 60U 

3 

0.5 

0.496 

0.684 

4 

0.375 

0.  455 

0.439 

5 

0.375 

0.405 

0. 404 

6 

0.  312 

0.358 

0.  362 

7 

0.312 

0.340 

It  has  been  noted  previously  that  there  is  no  crosstalk 
in  the  reception  of  the  information  in  the  i^  carrier  due  to  the 
complementing  of  the  other  signals  on  the  channel.  However,  from 
Table  6.1  it  is  noted  that  the  correlator  output  (normalized)  de¬ 
creases  as  the  number  of  users  increases,  resulting  in  an  effective 


decrease  in  the  signal  to  noise  ratio.  !Eie  loss  in  signal  to  noise 
ratio  due  to  the  hard  limiting  appears  to  be  lower  for  non- binary  sig¬ 
nalling.  An  interesting  comparison  can  be  made  between  the  performance 

21 

for  N  *  2  and  N  =  4-  It  has  been  shown  that  the  probability  of  error 
per  bit  of  information  is  lower  for  N  =  4  than  for  N  =  2  when  the  sym¬ 
bols  for  N  =  4  are  mapped  onto  phase  modulated  cosine  functions.  Since 
the  signal  to  noise  ratio  loss  due  to  hard  limiting  is  lower  for  N  =  4, 
it  is  apparent  that  significant  improvement  in  multiplexed  performance 
is  obtainable  hy  using  quaternary  encoding  instead  of  binary  encoding. 


I 
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7.  SUMMARY  AND  CONCLUSIONS 

7.1  Summary  of  Major  Results 

The  major  objectives  of  the  thesis  were 

a.  The  synthesis  of  new  classes  of  periodic  sequences,  with 
binary  and  non-binary  symbols,  which  exhibit  useful  autocorrelation  and 
cross-correlation  properties,  when  the  symbols  are  appropriately  mapped. 

b.  The  analysis  of  two  types  of  multiple-access  systems 
using,  os  carriers,  these  sequences  along  with  some  known  classes  of 
sequences. 

In  Chapter  3  the  method  of  interleaving  two  level  binary  sequences 
was  analyzed  as  a  method  for  synthesizing  lax^e  classes  of  sequences  with 
prescribed  correlation  properties.  A  large  class  of  almost  two- level 
sequences  with  low  cyclic  cross-correlation  function  was  synthesized. 

These  sequences  are  applicable  for  synchronization,  some  asynchronous 
carrier  systems  or  as  a  set  of  non-linear  error-correcting  codes  with  good 
data  rates.  A  class  of  sequences  exhibiting  autocorrelation  functions  with 
intermediate  minor  peaks  all  of  different  amplitude  and  cyclic  cross-corre¬ 
lation  functions  which  are  uniformly  low  was  synthesized,  lfcese  sequences 
are  applicable  for  synchronization  under  high  signal  to  noise  conditions 
and  also  as  carriers  for  the  asynchronous  multiplexing  system.  It  was 
shown  that  somt  classes  of  interleaved  sequences  are  easily  generated  by 
a  non-linear  filtering  technique. 

In  Chapter  4  classes  of  binary  and  N-ary  sequences  which  are  cyc¬ 
lically  orthogonal  (orthogonal  for  all  cyclic  shifts)  were  derived.  A  set 
of  binary  cyclically  orthogonal  sequences  are  the  k  basis  sequences,  (re- 


l4o 


ferred  to  as  the  Basic  Walsh  functions)  for  the  set  of  2^  Walsh  functions. 
All  of  the  classes  of  binary  cyclically  orthogonal  sequences  derived,  con¬ 
tained  sequences  all  of  different  least  period.  It  is  conjectured  that 
binary  cyclically  orthogonal  sequences  of  the  same  least  period  do  not 
exist.  This  conjecture  was  discussed  by  noting  the  Fourier  components  of 
mapped  binary  sequences. 

In  Chapter  5  the  problem  of  multiplexing,  on  a  linear  channel  with 
additive  white  Gaussian  noise,  the  carriers  from  k  different  transmitter- 
receiver  pairs,  when  each  transmitter  is  synchronous  with  its  corresponding 
receiver,  although  the  various  transmitter-receiver  pairs  are  asynchronous 
with  each  other,  was  analyzed.  Information  is  transmitted  by  sending  a 
carrier  or  its  negative.  It  was  shown  that  the  average  binary  error  prob¬ 
ability  when  the  k  carriers  are  a  set  of  mapped  randomly  chosen  binary 
sequences,  reduces  to  the  binary  PSK  formula  when  the  ratio  of  number  of 
sequences  to  sequence  length  is  zero  as  the  sequence  length  becomes  un¬ 
bounded.  The  average  probability  of  error  when  the  carriers  are  mapped 
Basic  Walsh  functions  was  derived,  and  it  was  shown  that  the  result  is  a 
function  of  the  particular  sequence  assigned  to  a  user,  and  the  number  of 
users  sharing  the  channel.  The  average  probability  of  error,  averaged  over 
all  of  the  cariers,  approaches  the  PSK  formula  as  the  number  of  users  (and 
of  course  the  sequence  length)  becomes  unbounded.  However,  only  loggL 
(with  L  “•  <*>)  users,  can  be  present  when  the  carriers  are  mapped  Basic 
Walsh  functions,  while  a  greater  number  can  be  present  if  mapped  random 
binary  sequences  are  used.  The  average  probability  of  error  was  computed 
when  the  carriers  are  a  set  of  sinusoids  of  periods  comparable  to  the 


periods  of  the  Basic  Walsh  functions.  Cie  average  probability  of  error 
was  higher  for  the  sinusoid  carriers. 

A  set  of  mapped  Bose-Chaudhuri  error-correcting  codes  was  consid¬ 
ered  as  carriers.  The  noiseless  channel  performance  of  these  carriers 
appeared  to  be  significantly  worse  than  the  performance  of  the  Basic  Walsh 
functions. 

Chapter  6  dealt  with  the  "hard- limiting"  multiplexing,  prior  to 
transmission,  of  a' set  of  mapped  N-ary  sequences,  with,  all  carriers 
synchronous.  Information  is  transmitted  by  sending  the  mapping  of  a 
sequence  or  one  of  the  N-l  complements.  The  optimum  set  of  sequences  were 
the  k  cyclically  N-orthogonal  sequences  of  length  N^  derived  in  Chepter  4. 

It  was  shown  that  the  signal  to  noise  ratio  degradation  due  to  "hard- 
limiting"  multiplexing  was  lower  jbr  K=3  and  4  than  for  the  binary  case. 

7.2  Problems  Y.c t  to  Be  Solved 

In  this  section  a  number  of  unsolved  problems  suggested  by 
the  research  described  in  this  thesis  are  listed. 

a.  The  synthesis  of  two- level  N-aiy  sequences  for  both 
mappings  presented  in  Chapter  2. 

b.  A  listing  of  all  of  the  autocorrelation  functions  which 
arbitrarily  interleaved  two- level  sequences  can  exhibit. 

c.  A  search  for  large  classes  of  sequences  (binary  and  N-ary) 
which  exhibit  good  cyclic  cross-correlation  properties,  besides  the  Bose- 
Chaudhuri  codes,  interleaved  sequences  and  cyclically  orthogonal  sequences 
considered  in  the  thesis.  It  appears  that  the  technique  of  choosing 
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sequences  with  distinct  Fourier  components  will  be  useful  for  the  synthe¬ 
sis  of  binary  sequences. 

d.  Investigation  of  cyclically  orthogonal  binary  sequences 
of  the  same  least  period. 

e.  Additional  searching  for  sequences  for  the  asynchronous 
linear  multiplexing  system  which  provide  performance  comparable  to  the  set 
of  "average"  randomly  chosen  sequences.  The  computer  derivation  of  the 
"flipped"  cross-correlation  function  of  mapped  Bose-Chaudhuri  codes,  and 
the  interleaved  sequences  considered  in  Chapter  3,  might  indicate  that 
these  carriers  are  applicable  for  this  system. 

f.  An  analysis  of  the  asynchronous  linear  multiplexing  of 
mapped  non-binary  sequences,  with  a  mapped  sequence  or  complement  sequence 
sent  as  information. 

g.  Derivation  of  sequences  for  an  asynchronous  "hard-limiting11 
multiplexing  system. 

7.3  Relationship  Between  Thesis  Problems  and  Prior  Work 

In  this  section  the  problems,  results  and  proofs  presented 
in  this  thesis  are  related  to  the  prior  work  which  hat  appeared  in  the 
literature. 

In  Chapter  2  the  concept  of  the  N-l  complement  sequences  of  an 

N-ary  sequence  was  introduced.  Although  the  existence  of  complement 
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sequences  has  been  "hinted"  at  in  the  literature,  the  technique  of  sig¬ 
nalling  with  mapped  complement  sequences  has  not  been  analyzed  to  the  ex¬ 
tent  that  binary  signalling  has.  The  proof  of  the  shift  and  add  property 


of  p-nary  maximal  length  sequences  and  the  derivation  of  the  autocorre¬ 
lation  function  of  the  iaximal  length  sequences  with  symbols  mapped  to 
-1,  0,  +1  haa  not  appeared  previously.  The  shift  and  subtract  property  is 
probably  known. 

The  analysis  of  the  correlation  functions  of  mapped  interleaved 

two-level  sequences,  presented  in  Chapter  4,  appears  to  be  a  problem 

which  has  not  been  previoulsy  solved.  Interleaving  has  been  used  pre- 

29 

viously  for  the  derivation  of  large  classes  of  error-correcting  codes 
and  the  product  "almost"  two- level  sequences  which  Titsworth  has  derived 
are  recognized  as  interleaved  complement  sequences.  Ihe  large  class  of 
"almost"  two-level  sequences  and  "different"  peak  sequences  derived  in 
Chapter  4  were  not  previously  known.  The  non-linear  filtering  technique 
as  a  method  for  generating  certain  classes  of  interleaved  sequences  is  an 

19 

extension  of  the  work  of  Raphael. 

It  is  probable  that  the  cyclically- orthogonal  property  of  the  Basic 
Walsh  functions  is  known.  However,  the  general  techniques  for  the  synthesis 
of  cyclically  N-orthogonal  sequences,  presented  in  Chapter  4,  have  not 
appeared  previously  in  the  literature.  No  work  has  appeared  on  the  Fourier 
analysis  of  binary  sequences  or  on  the  existence  of  cyclically  orthogonal 
binary  sequences  of  the  same  least  period.  The  results  presented  on  the 
generalized  Hadamard  matrices  are  probably  known. 

Although  a  significant  effort  has  been  directed  towards  the  solu- 
tion  of  the  multiple  access  problem, the  problem  of  the  binary  asynchro¬ 
nous  linear  multiplexing  of  carriers  has  not  been  previously  considered. 


The  technique  used  for  the  asymptotic  series  approximation  to  the  proba¬ 
bility  of  error  integral  for  the  Basic  Walsh  functions  was  suggested  by 
the  work  on  a  similar  integral. 35>36  ^e  derivation  of  the  asynchronous 
performance  of  the  random  sequences  is  similar  to  a  derivation  of  the 
performance  for  a  synchronous  system.  37 

The  transform  analysis  of  N-ary  sequences,  presented  in  Chapter  6, 
utilized  in  the  derivation  of  the  "hard-limiting''  performance  of  N-ary 
sequences  mapped  onto  the  roots  of  unity  is  an  extension  of  the  work  of 
Tit  sworth. 12,23 
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APPENDIX  A:  Proof  of  Theorem  2.2 


as 


The  normalized  autocorrelation,  p/^(  X)  was  defined  in  Section  2 

>a(x)  -t-  l  4m)<bi-x>*  •  <A-« 

i=0 


For  the  mapping,  0  -  0,  1-1,  ?.  -  -1  of  the  ternary  maximal  length  se¬ 
quence,  Equation  (A.l)  is  formed  by  the  sun  of  products  of  the  form, 

r-i  r-i 

0  *  0,  1  *  1,  (-1)  *  (-1) .  However,  since  there  are  3  2's  and  3 

l's  in  a  ternary  maximal  length  sequence  it  is  seen  that 


P,(  0)  = 


fr-i 


,r-i 


(3-1) 


From  -Theorem  2.1  it  is  noted  that  the  mod  3  sum  of  the  maximal 
length  sequence  ani  the  maximal  length  sequence  shifted  by  L/2  digits, 
yields  the  all  zero  sequence.  Hence,  for  X  =  L/2,  Equation  (A.l)  is 
formed  by  the  sum  of  products  of  a  mapped  symbol  and  its  mapped  additive 
inverse .  The  products  appearing  are  1  *  ( -l) ,  ( -l)  *  1,  0  •  0  resulting 
in 


PA 


L 


The  result  for  the  remaining  integral  shifts  will  now  be  proven. 
It  is  noted  if  any  two  ternary  sequences  are  summed,  (mod  3)>  there  are 
three  sub-sums  which  generate  each  of  the  three  symbols  in  the  resultant 
sequence.  These  sub-sums  are 


0  +  1N;  1  +  1)  2  +  l\ 

1  +  0)  =  1  0  +  2  =2  1  +  2  =  0  . 

2  +  2)  2  +  0)  0+0) 
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It  will  now  be  established  that  if  a  ternary  maximal  length  se¬ 
quence  is  summed  with  a  shifted  version  of  itself  (except  for  shifts  of 
0,  L/2) ,  each  of  these  sub-sums  will  occur  3r’2  times,  except  the  last 
sum,  (0+-0)  which  will  occur  3r'2  -  1  times.  Once  this  is  established, 
it  will  be  clear  that  p  ( X)  =0,  X  /  0,  L/2. 

Write  the  original  maximal  length  sequence,  and  its  first  (r-l) 
cyclic  shifts  as  the  row  vectors  of  an  r  x  L  matrix  M. 


*0 

ai  a2 

1 

H 

1 

$ 

• 

• 

• 

aL-i 

• 

ao  ai 

•  • 

...  e-jj-2 

• 

• 

j^-r+i 

•  # 

^-r+a  aL-r+3 

• 

The  row  space  of  M  is  a  basis  for  the  L  cyclic  shifts  of  the  original 
sequence,  where  the  vector  generated  by  twice  the  first  row  is  the  se¬ 
quence  shifted  by  L/2  digits. 

From  the  characteristic  of  the  maximal  length  sequence,  the 

j*  y 

L  =  3  "1  columns  of  M  are  the  3  -  1  nonzero  r- tuples  over  GF(3). 

Hence,  if  we  consider  the  corresponding  elements  of  row  1  and  any  other 
row  as  forming  L  2-tuples,  it  is  clear  that  in  this  set  of  L  2-tuples, 
each  of  the  9  possible  2-tuples  over  GF(3)  will  occur  3r_2  times  except 
(0,0)  which  will  occur  3r~2  -  1  times.  Thus  the  result  is  established 
for  X  =  1,  ...  y  r— 1. 

Next,  form  a  new  matrix  M’,  of  the  form 


I 


U7 


M' 


(L+D/3 


(I /H)/3 


0  0 


...  0 


11  ...  1 


(L+l)/3 


2  2  ...  2 


000  1  ...  1  2  ...  2 


0  0 


2  2  0  0 


2  2 


(L+D/9  (L+l)/9 


0  12 


...  2 


0  1  2  ...  2  0  1  2  ...  2 


by  permutation  of  the  rows  and  columns  of  M,  and  the  addition  of  a  column 
containing  all  zeros.  It  is  apparent  that  M1  can  be  formed  from  the 
characteristic  of  the  maximal  length  sequence.  Replace  any  row  of  M', 
except  the  first  row,  by  a  vector,  V,  which  is  any  linear  combination  of 
the  r  rows  except  twice  the  first  row.  If  the  vector  Is  viewed  as  divided 
into  3  segments  -  the  first  (L+l)/3  symbols,  the  second  (L+l)/3  symbols 
and  the  last  (L+l)/3  symbols  -  then  each  of  these  segments  will  contain 
(L+l)/9  ones,  (L+l)/9  twos,  and  (L+l)/9  zeros.  Thus  in  the  set  of  (L+l) 
2-tuples  which  are  formed  by  taking  corresponding  elements  of  row  1  and 
V,  each  of  the  9  possible  2-tuples  over  GF(3)  occurs  3r-2  times.  Discard¬ 
ing  one  of  the  (0,0)  2-tuples,  which  results  from  the  addition  of  the  all 
zero  column  in  M’,  it  is  seen  that  the  9  possible  sub-sums,  resulting 
from  the  summing  of  a  maximal  length  ternary  with  any  shifted  version  of 
the  sequence,  (except  L/2),  occur  3r‘2  times  except  the  0+0  sub-sum 
which  occurs  3r'2  "  1  times.  Hence  Equation  (A.l),  for  the  mapped  se¬ 
quence,  reduces  to 


[l*  1  +  (-D  •  (-D  +  (-D  *1  +  1*  (-1)]  =  0  . 


« 


PA(X) 

X/0,L/2 


L 


IkB 


l 


APPENDIX  B:  Proof  that  Evejy  2-Fold.  Interleaved  Sequence  Can 
Be  Generated  by  Multiple  Nonlinear  Filtering 

For  n  =  2,  the  output  sequence  of  the  multiple  nonlinear  filter, 
in  polynomial  representation,  reduces  to 

soT(x)  3  k270^^*2)  ©x2^1+1"y°^a1(x2)  ] 

,  ,  T  ( B.l) 

+  xi2^0+1[a0(x2)  ®  x2^1-?,0^a1(x2)  ]  ,  mod  (xsL-l)  . 

In  order  to  simplify  (B.l),  let  ■/[  =  7\  -  7Q.  Then  SoT(x)  becomes 

SqT(x)  3  x27o[a0(x2)  ©x2^7l+1^a  (x2)] 

,  S  T  ( B.2) 

+  x27o+1[aQ(x2)  0  x2^1^a1(x2)  ]  ,  mod  (x~L-l) 

It  will  now  be  shown  that  each  of  the  states  of  the  filter  gen¬ 
erates  a  different  output  interleaved  sequence.  Since  there  are  as  many 
filter  states  as  interleaved  sequences,  the  proof  will  be  complete.  The 
following  cases  will  be  considered  for  70>  7X>  a0(x),  ai(x)» 

Case  a.  aQ(x)  =  0(x)  ,  ax(x)  =  0(x) .  Then  S0rj.(x)  =  0(x2)  +  xO(x2),  in¬ 
dicating  that  the  output  sequence  is  the  all  zero  sequence. 

Case  b.  aQ(x)  =  0(x),  ai(x)  =  hj^x) .  Then  SQj(x)  =  x2bh  (x2) 

+  x2^”1^hm(x2) ,  mod  (x2L-l),  where  b  =  70  +  7X  +  1,  where  b  =  0,1,  ..., 
L-l.  The  interleaved  sequences  for  different  values  of  7 i  are  all  differ¬ 
ent,*  and  for  no  value  of  7i  can  the  all  zero  sequence  of  Case  (a)  result. 


*The  result  is  independent  of  70  since  sw0  is  open  (b  =  70+7 {+1  =  7i+l) . 
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Case  c.  a0(x)  =  h^x),  ax(x)  =  0(x).  Then 

SoT(x)  s  x2Chm(x2)  +  x2c+1hm(x2)  ,  mod  (x21-!)  , 

where  c  =  70,  c  =  0,1,  . ..,  L-l.  The  interleaved  sequences  for  different 
values  of  70  are  different,  and  for  no  values  of  y0  can  the  sequences 
of  Cases  (a)  or  (b)  result. 

Cased.  aQ(x)  =hm(x),  ax(x)  =  h^x) .  Then 

soT(x) 

+  x2[ro*i(rl)>lhn(j[2)  (  mod  B'} 

It  is  clear  that  thi3  polynomial  cannot  reduce  to  the  polynomial  repre¬ 
sented  in  Case  (a).  In  order  to  establish  that  Case  (b)  cannot  result. 

It  is  necessary  to  demonstrate  that 

yQ  +  1(7^1)  t  70  +  i(y[)  +  1  , 


or  equivalently  that, 

i(7i+D  h  i (y[)  +  1 

The  properties  of  the  function  I(X)  defined  by  (2.8)  will  now  be  used. 
From  (2.8)  we  know  that 

xI(7i+1)hm(x)  s  hjx)  ©x^(x)  ,  mod  (x^l)  , 

and 

x1^)^*)  =  xhntx)  ©x^hnfx)  ,  mod  (xL-l)  . 


However,  from  the  uniqueness  of  the  shift  and  add  property  it  is  clear 
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that  x^^1+1^hn(x)  f  xI(?'1)+1hm(x)  f  m0<i  (x^-l)  indicating  that  I(y^+1) 

/  I(7>1). 

In  order  to  show  that  Case  (c)  cannot  result  from  Equation  (B.3) 
it  is  only  necessary  to  indicate,  because  of  uniqueness,  that  I(7i+l) 

/  1(7;). 

As  a  final  step  it  vill  be  shown  that  the  L2  different  combina¬ 
tions  of  7o  and  7i>  in  (B.3)  each  generate  different  interleaved  sequences. 

For  y[  =  L-l,  S 0t(x)  reduces  to 

Sot(x)  h  0(x2)  +  x2^o+I(L"1)]+i  hjjjCx2)  ,  mod  (x^-i)  , 

which,  for  the  L  values  of  yQ,  represents  all  cyclic  shifts  of  the  maxi¬ 
mal  length  sequence  in  the  "second"  position  and  the  all  zero  sequence 
in  the  "first"  position. 

For  y[  =  0,  SoT(x)  reduces  to 

S0«p(x)  h  a^fro+Ki)]^/^)  +  xOCx2)  ,  mod  (x^-l)  , 

which  represents,  for  the  L  values  of  70,  all  cyclic  shifts  of  the  maxi¬ 
mal  length  sequence  in  the  "first"  position  and  the  all  zero  sequence  in 
the  "second  position. 

For  7^  ^  0,  L-l  the  output  sequence  S0j  is  an  interleaving  of 
shifted  maximal  length  sequences.  In  order  to  show  that  these  interleaved 
sequences  are  all  different,  it  is  necessary  to  establish  that,  if  70  /  6C, 
and  y l  £  6^,  then  both  of  the  following  equalities  cannot  be  satisfied 
s  inrultaneously . 
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yQ  +  I(7i+1)  =  50  +  I(6i+1)  ) 
7q  +  I (7[)  =  60  +  I(6i)  ) 


(B.4) 


This  is  accomplished  by  showing  that  only  d  =  0  satisfies  the  equality 


I(ri+1)  -  i (j[)  =  Kri+1+d)  -  Kri+d) 


(B.5) 


Equation  (B.5)  can  be  transformed  to 


I(7{+1)  +  I(7i+d)  -  1(7 1*  1+d)  +  Kri)  •  (B.6) 

From  the  definition  of  I(X),  (B.6)  implies 

xK'1+  ^^(x)  ©x7l+dbia(jc)]  =  xI^7l^[hm(x)  ©x7l+_+dhm(x)  ]  ,  mod  (xL-l)  , 


or 


:I(W«  •<  W\<*)  -  *i(y\(*W  (xl-i). 


or 


xd^xI^1+  ^+Ylhw(x)  0xI^^+Yl+1l^(j;)l  s  xI^^bW©  xl(Y4m(4  mod  (xL-l) 

(B.  7) 

However,  (B.  7)  can  be  transformed  to 

xd  {x7i  [hj x)  ©:cVl+1hm(x)"l©xYl+;L^hm(x)  ©  xYlhm(x)Jj 


5  [hm(x)  ©xYl+1hra(x)  ©hm(;c)  ©xrihm(x)  j  ,  mod  (xL-l) 


m  w  m 


or 

xdi  h  (x)  ©::Yl+1h  (x)l  =  jx^h  (x)  QxYl+1] h  (x)  1  ,  mod  (x^-l) 
l  m  m  j  L  m  n  <  ' 


(B.8) 


From  (B.3)  it  is  seen  that  only  d  =  0  satisfies  (B.5),  establishing 
the  proof. 
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APPENDIX  C:  Derivation  of  "Flipped"  Cross -correlation  Function  for 
Basic  V7alsh  Functions 


( f) 

In  this  appendix,  (X^.)  will  be  derived  for  i  >  i  and 
=  n  2in  +  \r  where  n  =  0,1,  ...  ,  2^ 1-1  -l,  and  w  =  0,1,  ...  ,  '21. 
The  derivation  for  the  values  of  not  included  above  is  quite  similar 
and  only  the  final  result  will  be  indicated. 

Let  Sa_  -  [(ai)c,  ...  ,  (a.)b,  ...  ,  (a*)^},  =  [(aJQ, 

•••  $  **•  >  ( a^,) L-i^  anc*  =  £go>  •••  >  •••  >  Gl-i  J  •  Tlien 

the  first  ,\  .  symbols  of  S  ,  [(a.).,  b  =  0,  . . .  ,  kc  -  1],  are 

~  a^  id  *> 


b  =  0,1,  ...,  2^1;  2i+1,  ...,  3«2X-1;  ...  ; 
n2i+1,  n2i+1+l,  . ..,  n2i+1  +  w-1 

b  =  21,  2i+l,  ...,  2i+1-lj  3-21,  . U^-l*  ...; 
(2n-l)21,  (2n-l)2i+l,  . ..,  n2i+1-l  . 


\s  (f) 

The  first  X.  symbols  of  x  S'  , 
**  aI 

are  (a  )T  ,  =  1,  c  =  0, 1,  ... 

S  =  S  +  ^  S^  are 

g  at  •  aiX, 


,+c  C  “  °*  '  Xi  " 

Xj 

X/.-1.  The  first  X^  symbols  of 


d  =  0,1,  •  • ,  21-!;  2i+l ,  3’  2^1;  ...  ; 

n2i+1,  ...,  n2^+1  +  w-1 

d  =  21,  2i+l,  ...,  2i+1-l;  Z’21,  ...,  V2X-1;  ...; 
( 2n-l)21,  . ..,  n2i+l-l  . 


I 
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From  the  symbol  structure  indicated  above,  it  is  seen  that  there  ore  1  y' 
more  ones  than  zeros  in  the  first  \  =  nl’i+1  +  u  digits  of  S  .  Thus 


from  ( 5.7 ),  ;.<f^(n2i+1  +  w)  =  -  2y/L  (L  =  2k)  for  n  =  0,  ...  ,  Z-1'1 -1, 
and  v  =  0,1,  ...  ,  2*. 

For  the  same  constraint  on  n,  but  yith  u  =  2 i,  2*  +  1,  . ..,  2i+1 , 
the  "flipped"  cross-correlation  is  found  to  be,  -^(n2i+-  +  y)  = 

I  a/ 


rr 

o 


-  2(2i+:-v)/L. 


-  - - i-1  pi-i-i  J, 


For  n  =  2  1-1 ,  2 


a.  1  2^” 1  -  1 

JL  f  •  •  •  f  w  Xf 


p(f)[n2i+1'™)  = 
1" 


+  2v/L 


2(2i+1-w)/L 


w  =  0, 1,  ...  ,  2 


y  =  21,  ...  ,  2i+ 


Finally  for  n  =  2a(2A'-*-1  ),  2o(2'v~*"’- )  +  1,  ...  , 
2o(2^ i”1  )  +  a*’"1”1  -  1,  (o  =  0,1,  ...  ,  2k-"-:-  1), 


p(^(n2i+1+y)  =  ) 

ix  \ 


f  -  2w/L 


2(21+1-u)/L 


v  =  0,3,  ...  ,  2 


w  =  21,  ...  ,  2t+1, 


and  vith  n  =  (2o+l)2Ui~l,  ...  ,  ( 2a:l)  2~_i-x  +  2X'i-1-  1, 


p<f)(n2i+1+v)  = 

1  A 


2y/L 


^  +  2(2in-y)/L 


w  =  0, 1,  ...  ,  2 

v  =  21,  ...  ,  2i+1 
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APPENDIX  D:  Derivation  of  Probability  of  Error  With  Basic  Walsh 
Functions  For  k=2 

In  this  appendix,  the  average  probability  of  error,  when  corre¬ 
lating  with  the  carriers  associated  with  the  Basic  Walsh  function 
sequences  Sg^  or  Sa^,  with  two  carriers  on  the  channel,  (k=2)  is  deter¬ 
mined.  We  are  given  the  density  functions,  for  i=0,  jt=l, 


P\Pon) 


(D.l) 


and 


pCp^)  =  p(Poi)  =  \i  6(Poi)  +  — 

'  2S 

0 


-  S  <  Pol  <  S  (D.2) 
2  “  ”2 
elsewhere 


Then  p(Z  |  do=+l)  is  given  by  the  convolution  of  p(pon)  and 
p(poi)>  or 


P(Z 


|  di=+D  = 


p(Poi  )Pon(z"Poi  )dPoi 


u 

-oo 


ps/2 

5 '6(  poi )+.-.  ] 

-S/2  '  SJ 


exp 


(D-3) 


Thus  the  expression  for  the  probability  of  error,  P^)  =»  P(ej)  becomes 


P(ei)  =  P(eo)  =  p(Z  |  d^+lJdZ  = 

-JCO 


■(.^R0i/:.s£~|dpold2 

%S/T  J 


(D-10 


1_ 

hS 


S/2 


-S/2 


A,s/t 


where,  erf(v),  the  error  function  of  v  is  defined  as 


U  T 

2  r  -x3 

erf(y)  =  -a  e  dx  . 
«/tt 


(D-6) 


The  integration  of  the  error  function  in  (D-5)  is  performed  by  parts. 
Ihe  final  result  for  the  average  probability  of  error  is  then 
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APPENDIX  E:  Asymptotic  Series  Approximation  to  P(e0) 


An  asymptotic  series  approximation  to  the  integral  expression 
(526)  for  i  =  o,  which  yields  accurate  results  for  k  >  4  will  now  be 
derived.  For  i  =  o  (5.26)  can  be  transformed  to 


co  .  P 

1  1  sin  2k-1y  f  sin  y  l"^"1  -y22sk-4/G 

2  "  n  y  L  2y  2J  6 

o 


(E.l) 


The  procedure  to  be  used  for  the  above  expression  is  sim''  to  the  pro- 

35.36 

cedure  previously  used  for  the  asymptotic  series  approximation  to 


CD 

T  f  sin  x  'Nn 
J  V  x  J 


o 


dx 


An  approximation 


to 


fsin  y 
L  2y 


will  first  be  derived. 


is  noted  that 


It 


d  T a  /'sin  y  .  P]  y  cos  y  -  sin  y 
2y  zA  =  y  sin  y  +  y- 


(E.2) 


Writing  the  series  approximation  to  the  right  hand  side  of  (5.2)  we 
find 


Performing  the  long  division  and  integration,  and  retai  ling  the  first 
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three  terms 


tn 


i  -s.i.n.  y 

V  2y 


(E.4) 


Thus: 

f sin  y  1  >ik-1  -(v*-l)y°/l2  k"1[‘]Jfo  +  iB^lU h] 

\  2y  2  J  “  e  e 


(E.5) 


which  can  be  approximated  as 


-(k-l)r/l2  r  (k-l)y«  ,  (k-l)y‘,-| 

V  2y  2  J  **  L  1440  18,144  J 


Substituting  (E.6)  into  (E.  1)  the  expression  for  P(eQ)  is  approximated 
as 


1 

n 


oo 


P  sin  2k~1y 


exp 


i 

o 


00 

r  [sin  2k-y 
o 


00 


-  \  j [sln  2k'lyJ[i57I^-]  exi>  [-  y3  C.rl 


k-  -** 


(E.7) 


The  contribution  of  the  third  integral  can  be  shown  to  be  negligible 
for  0  <  G2  <  20.  Also  for  G2  in  this  range  and  k  >  4  :.e  car.  use  the 


approximation 


15° 


e:cp 


r 


r— 

V12 


Jk- 


G= 


exp 


r  -jlI 

1  G= 


2  k- f 


n 


fi 


12 


(E.8) 


Incorporating  these  approximations  into  (E.  j)  and  performing  the  ir.te- 

34 

gration  with  the  aid  of  integral  357- 2b  in  Grobner,  we  get  the  approxi¬ 
mation 


(jt-1) 

lMiOit 


G& 

2'^k-tr 


+ 


(k-1)3  G7 

(T^SoJn  FKOT  l  y 


(E.9) 


where  F[x,y,  z]  is  a  confluent  hypergeometric  function. 

It  is  seen  that  the  resultant  approximation  for  P(e^)  is 
the  binary  PSK  error  probability  plus  some  error  terms  which  approach 


zero  for  k  -  co. 
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APPENDIX  F:  Probability  of  Error  Averaged  Over  k  Basic  Walsh 
Function  Carriers 


In  this  appendix  the  probability  of  error  averaged  over  all  k 

carriers  on  the  channel,  for  k  -•  ®,  will  be  derived.  When  correlating 

with  the  carrier  associated  with  S_  =SQ  ,  the  output  of  the  ith  corre- 

ai  ak-E 

lator  becomes  Gaussian  distributed  with  mean  zero  and  variance,  , 

^k-E 

E=l,  2,  .  • . ,  k  -  ao  , 


=  M  ♦  £ 

\.l  2T  3 


(F-l) 


Let 


cl  =  cl  -M  .  (F-2) 

E  ^-E  2T 

Then  the  probability  of  error  P(e  )  when  correlating  with  the  i=k-E 
carrier  is 


(r-3) 


where 


G  = 
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We  are  interested  in  Prrl  the  error  probability  averaged  over  all  of  the 

n 

carriers.  Thus 


P  =  \  P(e  ) 

w  k-®  k  2,  1  k-E;  9 

E=1 


(F.k) 


or 


’■ ' "  { £<*  - 1  "fcsj  * 


E 

*  +  ~F 


(F.5) 


Expand  the  error  function  in  the  series 


«  2c -1 

2  c  .  .c+1  xt 


erf(^)  =  1  '-1’ 


c=l 


(2c-i)(c-i): 


(F.  6) 


where 


= 


E  7i  2^3 


L?+  ■#■] 2  t + 


2aE2G2  -ii 


S3 


•  * 


( F.  7) 


2a/ G2 

If - <  i  then  x_  can  be  expanded  in  a  power  series  in 

^  2c r/G3  2a/ G2 

ascending  powers  c  — -g —  .  Now  „  it-.. 

o  S3 


achieves  its  maximum  value  for 


161 


E=1  or  E=2, 
aler.t  to 


with 


1 


9 


Hence  for  all  G3  <  9*,  x^,  is  equiv- 


x 


E 


(F.8) 


Then 


_  lim  1 
oo  k 


£{** i 


E=1 


2_  ?  (-l)0^20-1 

v  c=l 


(F.9) 


Since  all  of  these  infinite  series  are  uniformly  convergent  we  can  inter¬ 
change  the  summation  and  limiting  operations  in  an  arbitrary  manner.  Thus 
the  expression  for  Pw  c*m  be  reduced  to 


p„  *  {i  -  i  ■  2-  £ 

^  ./tT  ,  (2c- 


c+1^2c-l 


00 


cti  (2c-i)(c-i): 


lim  1 
k-« 


E=1  d=0  *  « 


(F.  10) 


*The  average  probability  of  error  for  binary  PSK  signalling  with  G  =9 
is  1.1  •  10-fe. 
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It  will  now  be  shown  that 


00 


11m  1 
k- 


2ovzG2  .d-Cc-l 


=  1 


The  first  term  of  the  above  summation  is  1, 

A  typical  remaining  term  is  of  the  form 


lim  c 
k-“  k 


P 


> 


where  0  is  a  positive  integer  and  C\  is  a  finite  constant. 
But, 


Q/2-^[£(@)(e-i)V] 

m=0 


It  is  then  verified  that 


)  Pw#  averaged  over  the  k  Basic 


v 


(F.  11) 


(F.  12) 


(F.13) 


Walsh 


Hence  the  probability  of  error 
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functions,  k  -*  ®,  is 


00 


£.  iS i=iL 

ttA  (?c_ 


c+1  2c-l 


^  (2c-l)(c-l): 


=  £  -  ?  erf  (G)  , 


for  G2  <  9. 
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